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Introduction 



Cosmology and Quantum Gravity are the two main areas of Physics the research 
collected in this dissertation is about. 

Cosmology is the study of the universe as a dynamical system. It is a rather 
peculiar chapter of Physics for more than one reason. The main one, I believe, is 
that we have but one universe (someone may not agree on this) and there is no 
way to do experiments on the universe, in the proper sense of the word. We can 
at best observe it, and that's what astronomers do very well. Another reason, 
at least for myself, is on purely intellectual basis: Cosmology is one of the areas 
(if not the area) of Science that most touches in deep the questions Science was 
born for. One more argument, which is somehow mixed with the preceding one, 
is that Cosmology is in a certain sense a rather young science: Einstein's General 
Relativity (1916) is the first theoretical framework that has allowed a scientific 
and mathematical approach to the problem of the study of the universe as a 
single system. Before that, questions like 'where does the universe come from?' 
and 'what is the fate of the universe?' where mere philosophy, with no scientific 
attempt to be answerecQ 

Since its birth. Cosmology has taken gigantic steps, both on the phenomeno- 
logical and theoretical viewpoints. The Friedmann models (see section [T]) are the 
first theoretical models comprehending the dynamical expansion of the universe, 
together with the Big Bang initial singularity, both compatible with observations 
and accepted from the scientific communitjj^ 

I dare say that expansion and initial moment of the universe have somehow 
given birth (for what it is possible) to two macro lines of research. Observations 
of the Cosmic Microwave Background have indeed fostered the research on the 
nature of the universe in its very beginning. We talk then about "early-time 
Cosmology", referring to the study of primordial perturbations amplified by some 

^ Maybe it is superfluous to remark that they are not answered even nowadays. The point are 
not the answers, which will probably never come, but the scientific and rigorous attempt 
to give them. 

^The Big Bang, being a singularity of the theory, cannot be observed, strictly speaking. What 
is accepted, is the existence, some 15 billion years ago, of a phase of the universe in which 
everything (spacetime itself) was contracted in a tiny and extremely hot bubble. But on 
the nature of that bubble, on its real size and so on, the debate is very far from being 
closed. 



1 



Introduction 

(quantum) mechanism (e.g. inflation), giving birth to the large-scale structures 
we observe nowadays, such as galaxies, clusters and so on. On the other hand 
we have the discovery of the cosmic acceleration [6l [7], that has guided many 
theoreticians in the study of cosmological models able to reproduce this (and 
others to it related) experimental evidence. 

The first part of this collection is focused on the second of this problem, the 
one regarding the entire evolution of the universe, its expansion, its acceleration 
and so on; while I will not touch the cosmological perturbations or amplification 
mechanisms such as inflation. 

Quantum Gravity is not a self-consistent theory, yet. It is a work in progress 
attempt to find a coherent set of tools in order to formulate some kind of predic- 
tions at physical scales where both Quantum Mechanics and General Relativity 
should hold. And we know that taken as they are. Quantum Mechanics and 
General Relativity cannot be both valid. Thus, we cannot properly speak of a 
theory, since a theory is something self-contained (at least at some degree) that 
has the possibility of making scientific predictions. 

Loop Quantum Gravity is one of these attempts, so one usually says that it 
is a candidate for Quantum Gravity. The keystone of Loop Quantum Grav- 
ity, namely the feature that identifies it among other candidates, is to quantize 
(canonically) the metric tensor (or, better, an appropriate manipulation of the 
metric tensor) without assuming the existence of a somehow fixed classical back- 
ground, around which fluctuates the "quantum part" of the metric. This is known 
as background independence. 

Since its birth in 1986, LQG has reached amazing and important results, the 
main of which - I dare say - is the suggestive prediction of the 'granular' nature 
of space: LQG indeed predicts the existence of a minimum length in space. 
However, much is still to be done; and the main drawback, which is actually 
common to all the candidates Quantum Gravity theories, is the lack (more or 
less total) of testable predictions. 

A decade ago a new "spin-off" theory was born in the framework of LQG: 
spinfoam theory. It is an attempt to define quantum gravitational amplitudes 
between space-geometries (given by LQG) in terms of a sum-over-histories. The 
histories are called 'spinfoams', a kind of bubble-like representation of spacetime 
at fundamental level. Many problems are still open, and it is a very stimulating 
field of research, particularly for young researchers, since there are entire lines 
still to be explored and thousands of bridges with other theoretical frameworks 
still to be built. 
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The main results of the research here presented are the following: 

• we analyzed in depth two-field cosmological models, with one scalar and 
one phantom scalar field. The procedure of reconstruction of two-field po- 
tentials reproducing a dynamical evolution (specifically an evolution with a 
crossing of the phantom divide line) has been worked out: we have discov- 
ered that two-field models have a huge freedom with respect to single-field 
models, namely there is an infinite number of potentials which give, with 
a specific choice of initial conditions, the wanted dynamics for the model. 
Moreover, a thorough analysis of the phase space of two specific two-field 
models is carried out, showing that inside a single model, varying the ini- 
tial conditions on the fields, qualitatively different families of evolutions 
are present, with different cosmological singularities as well. 

• In line with the preceding result, we have coupled our two-field models with 
cosmic magnetic fields, which are experimentally observed quantities. We 
have shown the sensitivity of the amplification of the magnetic fields to the 
change of underlying cosmological model, but keeping the same background 
evolution of the universe. This, in principle, is a way to experimentally 
discriminate between cosmological models with different potentials. 

• The greatest problem of phantom scalar fields is quantum instability. We 
have shown that in the framework of PT symmetric Quantum Theory 
(somehow adapted to Cosmology) this instability can be cured: we have 
an effective phantom scalar field, with stable quantum fluctuations. 

• We analyzed one-ficid cosmological models comprehending a peculiar ver- 
sion of the Big Bang and Big Crunch singularity, namely singularity with 
finite and non-zero radius. We performed a detailed analysis of the phase 
space and reveal the presence of different classes of evolutions of the uni- 
verse. 

• In the framework of spinfoam theory for Quantum Gravity, we proposed 

a fixing of the face amplitude of the spinfoam sum, somewhat neglected 
until now, motivated essentially by a form of "unitarity" of gravitational 
evolution. 

I list here the publications in which the above said research has been collected 
(chronologically from the latest) : 

• E. Bianchi, D. Regoli and C. Rovelli, "Face amplitude of spinfoam quantum 
gravity" Classical and Quantum Gravity CQG 27, 185009 (2010) 
(arXiv: 1005.0764 [gr-qc]). 
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• A. A. Andrianov, F. Cannata, A. Y. Kamenshchik and D. Regoli, "Phan- 
tom Cosmology based on PT -symmetry", International Journal of Modern 
Physics D IJMPD 19, 97 (2010). 

A. A. Andrianov, F. Cannata, A. Y. Kamenshchik and D. Regoli, "Cosmol- 
ogy of non-Hermitian {C)PT -invariant scalar matter", J. Phys. Conf. 
Ser. 171, 012043 (2009). 

• F. Cannata, A. Y. Kamenshchik and D. Regoli, "Scalar field cosmologi- 
cal models with finite scale factor singularities", Physics Letters B 670 
(2009) 241-245 (arXiv:0801.2348vl [gr-qc]). 

• A. A. Andrianov, F. Cannata, A. Y. Kamenshchik and D. Regoli, "Two- 
field cosmological models and large-scale cosmic magnetic fields". Journal 
of Cosmology and Astroparticle Physics (JCAP) 10 (2008) 019 
(arXiv:0806.1844vl [hep-th]). 

• A. A. Andrianov, F. Cannata, A. Y. Kamenshchik and D. Regoli, "Recon- 
struction of scalar potentials in two-field cosmological models". Journal of 
Cosmology and Astroparticle Physics (JCAP) 02 (2008) 015 
(arXiv:0711.4300 [gr-qc]). 

This work is thus intended to pursue two targets: 

• to provide, for the fields of interest, a sufficient framework of notions and 
results that are the current background on which research is performed 
and are necessary to understand (or closely related to) the subject of the 
following goal 

• to resume and explain the results of the research I have been working on 
during this three-year PhD course, under the supervision of Dr. Alexan- 
der Kamenshchik. 

Gravity is the main character of the play. Matter is the other one. 'Matter' 
of course in a cosmological sense, i.e. everything that is embedded in spacetime. 
Talking in field theory jargon, we have nothing but fields: the gravitational field 
(the metric) and all the others: the matter fields (including those whose (quan- 
tum) excitations are the well-known elementary particles, and possibly others, 
which should model/explain such things as the dark matter, dark energy and, 
generically, other types of exotic matter/energy). All living on the spacetime 
manifold. 

We know that if we stay away from the quantum regime. General Relativity 
is the stage on which our characters plajj^ GR, in the intuitive rephrasing by 

■^There are of course many proposals to modify somehow GR, e.g. the Hof ava-Lifshitz pro- 
posal [9 , the f{R) models jTD]|Tl|, and others, but I will not consider them here at all. 
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Wheeler [8j, is simply encoded in the following aphorism-like sentence: matter 
teaches spacetime how to curve, spacetime teaches matter how to move. The 
"how", is encoded in Einstein's equation and in general in the Einstein-Hilbert 
action for GR. 

Classical Cosmology is precisely the study of this interplay between spacetime 
geometry and matter, on the scale of the entire universe. The idea is to build 
models with various kinds of matter, possibly encoded in terms of fields, that are 
able to reproduce the observed data on the dynamical evolution of the universe. 
More specifically for what concerns this thesis, the key empirical observation is 
the cosmic acceleration j6l[7], and the possibility of a kind of "super-acceleration", 
that has to be somehow explained theoretically, which we tried to, for our (little) 
part. 

When one enters the quantum regime, attention must be paid. We know that 
there is a typical scale, the Planck scale, under which we expect both GR and 
QM to be valid. However there is no accepted framework for this scale. No 
complete Quantum Gravity theory exists, yet. Actually, it is the Holy Grail of 
contemporary theoretical physics. 

The second (shorter) part of the present thesis is entirely focused on one 
specific candidate of Quantum Gravity: loop quantum gravity. I will have a 
more appropriate occasion to discuss this thoroughly |TT| but for the moment 
let me stress that I believe it fundamental to study what happens to spacetime 
under the Planck scale in order to properly catch the nature of the relation 
between geometry and matter. I think Quantum Gravity is going to revolutionize 
our concept of the interplay between gravity and matter as much as Einstein's 
theory has drastically changed the idea of matter fields top of a passive, fixed 
background. My personal contribution in this respect, is in the framework of 
spinfoam models for quantum gravity: a kind of path-integral formulation of 
quantum gravity, based and founded on the results of canonical Loop Quantum 
Gravity. 

What follows is divided in: 

• Part |T], dedicated to classical Cosmology: it starts with an introductory 
chapter intended to provide the necessary concepts of classical Cosmology; 
this is followed by four 'research'-chapters, each dedicated to one of the four 
publications about Cosmology listed before (or, equivalently, see yy,^,||2], 
[4J). The second of these four chapters, namely the one dedicated to "PT- 
phantom Cosmology", includes also a short review of PT-symmetric Quan- 
tum Mechanics, whose ideas play a key role in the research explained in 
(what remains of) that chapter. 

• Part [TT| focused on loop quantum gravity and spinfoam models: here I ded- 
icate more time to properly introduce the basics of these theories. Namely, 
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four entire chapters are devoted to a review of fundamental concepts of 
this theoretical framework. Indeed loop quantum gravity stands on much 
more advanced and sophisticated basis than classical Cosmology, and, since 
the work I have done is about a rather technical aspect of spinfoam mod- 
els, I thought it absolutely necessary to properly review all the underlying 
framework, trying to be as self contained as possible. 

The last chapter of this second part resumes the results and explains the 
content of the paper [5j, about a proposal for fixing the face amplitudes of 
spinfoam models. 



A brief conclusion 9.4 is intended to give a global view ex post on all the issues 
discussed. 
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Part I. 

Cosmology and Dark Energy 

models 



1. Cosmology: the basics 



Cosmology is the study of the universe as a dynamical system. At a classical 
level the tools to handle this kind of study are given by the Einstein's General 
Relativity theory, in nuce by equations 

1 p SvrG 
R^lu - -^gfiuR = —^T^y , (1.1) 

where R^y is the Ricci tensor, g^^^ is the metric tensor, R = g^'^R^y is the scalar 
curvature and T^^^ is the stress-energy tensor. As for the constants, G is Newton's 
gravitational constant and c is the speed of lighlj^ 

Of course, to (hope to) resolve these equations - in order to catch the dynamics 
of the universe, encoded in the solution g^v{x) - one has to know the source 
term T^^ point by point in spacetime. This can be achieved only by drastic 
exemplifications on the nature of the source and of the spacetime. Namely one 
assumes a rather strict - but also very reasonable at this large scale - symmetry: 
roughly speaking the assumption is homogeneity and isotropy of space. More 
precisely one assumes that the spacetime 4d manifold admits a foliation is 3d 
spaces that are homogeneous and isotropic, i.e. one admits the existence of a 
family of observers that see uniform space-like surfaces. This can be rephrased 
in a more suggestive and simple way by saying that we exclude the existence 
of special regions in space. This is known as cosmological principle and, as we 



will see in a moment, it drastically exemplifies equations (1.1). Of course this 
principle is just a device through which one tempts to extract physics from a 
general model which would be by far too complicated. One can relax in many 
ways this principle, for example in studying cosmological perturbations or by 
admitting anisotropics in some ways [98j. However, for our purposes, it won't 
be necessary to abandon this principle. 

The assumption of such a principle results in the following form of the line 
element in spacetime (the Friedmann-Lemaitre- Robertson- Walker line element, 
or simply FLRW): 

ds^ = dt^ - a\t)dl\ (1.2) 



^The reader can find detailed informations in every good textbook on General Relativity. Let 
me just say that for general GR framework I really like the book by Hawking and Ellis [T^] 
while for classical Cosmology (and for much more) I find very useful the book by Landau 
and Lifshitz [13]. 
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where 

df = - r\d9^ - sin^ 9d^') . (1.3) 

1 — Kr^ 

Recall that k is the curvature parameter and equals 1, or -1 for a closed, open- 
fiat and open-hyperbolic universe respectively. In all what follows we shall take 
k = 0, i.e. we assume a (spatially-)fiat universe. This is in accordance with all 
the observational datsEl 



In equation (1.2) the only dynamical variable is the scale factor a{t), the only 



dynamical degree of freedom in classical (uniform) Cosmology. As is clear from 



(1.2), it represents the 'size' of the universe at a certain moment, and thus it 
gives information on the expansion or contraction of the universe itself. 

Recall that for a non-dissipative fluid one has T^j, = {p + s)u^j_Ui, — pg^v with 
p and e the pressure and the energy density of the fluid, respectively, and = 
dXfj_{s)/ds is the 4- velocity of the fluid. With this in mind Einstein's equations 



become much more simple than (1.1), namely 



3 - +3^ = SnGe , (1.4) 
2ad + a"^ + k 



a2 



-SttGp . (1.5) 



These are the well-known Friedmann equations. Actually one usually replaces 
the second of these equations with the simpler requirement of conservation of 
energy, i.e. 

e = 3-{p + e), (1.6) 
a 



and thus uses equations (1.5) and (1.6) as the set of two independent equations 
that governs the dynamics of the univers^ 



Solving the system (1.5), (1.6) for a dust-like fluid (a fluid with zero pressure, 
which is a good approximation for visible matter like galaxies, clusters of galaxies 
and so on) one obtains the well-known Friedmann models, summarized in figure 



Remark. These results, I want to stress it, are striking. Doing research in 
this field one is often 'obliged' to investigate technical (and often exotic) details 
somewhat losing the importance of the Friedmann models: with a couple of 
amazingly simple equations one is able to obtain a qualitative (past and future) 



^Actually it is true that fc = is coherent with all the data, but some argue this is only 
a proof that the space is locally flat, and still it may have a positive or negative global 
curvature. I am not going to discuss this (important) issue here. 

^Intuitively, one for determining the source (say e) in terms of the scale factor, and then the 
other to determine the scale factor itself. 
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a[t] 




evolution of the universe. Of course these results must be corrected in order to 
fit observations introducing different kind of matter /energy, must be generalized 
to embrace anisotropics or inhomogeneities and must be supported by some 
quantum gravity theory to investigate what physics is inside the initial (and 
possibly final) singularity. But much of classical Cosmology is already inside 
these simple models. 

In what follows we shall make a huge usage of the Hubble parameter 

Ht) = ^ , (1.7) 
a[t) 

in terms of which the Friedmann equations (1.5), (1.6) becom^ 



h' = e- ^ 
a 



2 ' 



and, for the fiat case. 



e = —3h{e + p) 



h'^ = e , (1.8a) 
i = -Sh{e + p). (1.8b) 



and thus can be written in the compact form 



h = ~l{e + p). (1.9) 



*I use here, and everywhere in what follows about Cosmology, the convention SttG/S = 1. 
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For sake of completeness, let me recall the form of the curvature in terms of 
the Hubble parameter: 

R=-6(^h + h^ + ^^ . (1.10) 
which is important to understand how h and a determines curvature singularities. 



1.1. Accelerated expansion 

Up till here we have talked about a universe filled with galaxies, and we have 
seen its large scale evolution, which is an expansion starting from the Big Bang 
singularity and ending eventually either in a future singularity after a contraction 
(the closed case) or in an endless cold expansion. 

However it is now well-known that this is not the case. The cosmic expansion 
is accelerated [HI E]. This is the empirical observation that acted like a spark 
for the explosion of a large amount of cosmological models trying to incorporate 
this effect |14]-|19]. [20t [21] . |22]-[32] in the decade after it was first discovered 
(1998). The first part of this thesis lays in this line. 

Now I shall introduce the general ideas and the framework for these models 
as well as the technical jargon necessary to understand them. 

Take a generic cosmic fluid, and write its equation of state as follows 

p = we , (1-11) 

The w parameter is crucial to discriminate between different kind of expansions, 
as well as of fluids, of course. Indeed, Friedmann equations tell us that in order 
to have acceleration (a > 0) one has to have 

w<-^. (1.12) 

Such a fluid is generically known as Dark Energy (DE) [ll]-[T8]. The w = —1 
(constant) case is the ACDM model of Cosmologjj^ and is actually the best 
candidate model to fit observational data, as well as the simplest DE model. 
As is well-known, at a sufficiently large time scale (when the baryonic and dark 
matter contribution becomes negligible) the evolution of ACDM is de Sitter-like. 

There are some motivations to investigate DE models different from the simple 
ACDM. Some more practical and some more "philosophical". Let me spend some 
words on this issue. The not interested reader of course can skip the entire 
following subsection. 



^Obviously we are here talking of the DE side of the story. In every decent cosmological 
model there must be also a baryonic sector as well as a dark matter sector, which, when 
talking specifically about DE, are often understood. 
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1.1.1. ACDM: ups and downs 

Let me remark one thing: ACDM is in accordance with all observational data. 
This, as contemporary theoretical physics teaches us, does not prevent theoreti- 
cians to investigate alternatives, if there are some motivations supporting this 
research. 

A practical (and I think more common-sense) motivation is the following: it is 
true that observations do not exclude ACDM, but there are other possibilities not 
excluded as well. Moreover, some observations indicate a best fit for w which is 
less then -l|^(we shall see that this is not a painless difference), and contemplate 
the possibility of "crossing" this benchmark w = —1, thus with a non-constant 
value ofw. In my humble opinion I think this is enough for researchers to explore 
alternatives, even if exotic; obviously keeping clear in mind that there is a good 
candidate model in accordance with data. 

Theoretical, or "philosophical", motivations have been put out as well (I draw 
arguments mainly by |56]). They are generically of two types. The first is the 
'coincidence problem': observations tell us that the dark energy density is about 
2.5 times bigger than baryonic + dark matter(dama) density (about 74% for 
DE and 26% for the rest). Since the baryonic and dark matter contributions 
dilute in time (as a~^^^) while the cosmological constant does not, it will even- 
tually dominate the evolution of the universe and its density should be much 
greater than the observed + dama sector, starting from some cosmic era on. 
The argument is that it is not likely that we are just in that era in which the 
DE density is comparable to the observed+dama one, it would be too much of 
a coincidence. Think of it as a kind of cosmological principle in time. Thus, 
the argument claims, ACDM must be wrong. This is a rather poor argument, 
in my point of view. First, because probability arguments are to be taken very 
carefully, always. Second, because it is quite clear that in a universe dominated 
by the cosmological constant, matter and thus us, wouldn't exist. This is a kind 
of anthropic (counter) argument, but - I think - of a very light and reasonable 
kind. 

The second type is based on QFT reasonings: the cosmological constant (i.e. 
the DE in ACDM) is a kind of vacuum energy; QFT 'predicts' the existence of 
a energy of the vacuum, but if we compare the QFT calculation with the (tiny) 
observed value for A they differ by 120 orders of magnitude. So ACDM must be 
wrong. For this argument there are counter-arguments as well: it is true that 
QFT predicts a vacuum energy, but it is actually a huge amount of energy, and 
it is not at all observed. If QFT vacuum energy was really there, any region of 
space with a quantum field would have a huge mass, and it would certainly be 
observed. The Casimir effect reveals only the effect of a difference (or, better, a 

^See I15]-|51] and |53]-IS5]. 
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change) in vacuum energy. Thus it is hkely that there would be some unknown 
mechanism in QFT that prevents the vacuum energy to gravitate, or protect it 
from huge radiative corrections, rather than ACDM is 'wrong'. 

Thus, I do think it is worthwhile exploring alternatives of ACDM, but mainly 
for practical and observational reasons. And still I consider the ACDM model 
as a good standpoint for classical Cosmology. 



1.1.2. The case -1 < u; < -1/3 

Let us briefly sketch the behavior of this class, i.e. the class of models with state 
parameter w of DE type but greater than -1. The calculations starting from the 
evolution equation (1.9) are straightforward 

hit) = ^ . (1.13) 

3{w + l)(t - to) 

with a Big Bang initial singularity, in the sense that 

a{t) — -> , (1.14) 

indeed 

a{t) = (t - toY/'^-^'K a = -^(t-to)-^. (1.15) 

But differently from the Friedmann models - where a is increasingly high going 
towards the initial instant - here a goes to zero as well. 

One could as well tempt to analyze cases where w varies with time. One 
remarkable example is the so-called Chaplygin gas |26], characterized by the 
following state equation 

P = -^, (1.16) 

with A a positive constant. This fluid gives a sort of interpolation between a 
dust-like era and a de Sitter era. 



1.1.3. The case w < —1: the Big Rip type singularity 

I usually call the evolution of the models with w constant and less than -1 
"super-acceleration". Indeed they are characterized by 

h>0, (1.17) 

that is 

a , , 

- > ^ . 1.18 
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The dynamical evolution reads 



a{t) = {t- to)-^/^l'"+^l , a = :^r^At - to)-^ ■ (1.19) 

Notice that now the exponents are both negative, thus 

a > oo , a > —CO . (1-20) 

These features - scale factor and its velocity both infinite at a finite instant - 
define the Big Rip singularity, which is typical of the super-accelerated mod- 
els |57t [58] . Sometimes, the fluid responsible of such an exotic evolution is 
dubbed phantom energy |36|-|45|. 

I hope it is now clear to the reader that the benchmark w = —1 (somehow 
represented by the ACDM model) discriminates between dynamical evolutions 
qualitatively different: cosmic acceleration with a Big Bang and a cold infinity 
expansion versus a "super-acceleration" with a different cosmological singularity. 
The line w = —1 has deserved a name, the phantom divide line. 

Remark {the phantom stability problem). Let me anticipate here the big draw- 
back of the phantom energy. A kind of energy with state parameter w < —1 
violates all the energy conditions prescribed by General Relativity (see e.g. |12|). 



Indeed, when a phantom energy is modeled by a scalar field (see section 2.1) it 
requires a negative kinetic term. Its hamiltonian becomes, obviously, unbounded 
and such a field is plagued by manifest quantum instability: its quantum fluctua- 
tions grow exponentially with time |59 t l60 | . This is actually the main critic done 
to the use of such things as phantom fields. And, I have to admit, it is a quite 
strong and reasonable critic. However, ghost fields like these, were studied long 
before the phantom field was first introduce; and in chapter |4] (or equivalently 
in [3j) I will describe a mechanism that could fix this plague, paying the price of 
enlarging the framework of Quantum Theory to its PT symmetric version. 
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2. Two-field cosmological models 



This chapter is devoted to reviewing and explaining the content of /Ij. The 
idea is to work with two scalar fields, one of which phantom, and try with this 
model to reproduce an evolution from a Big Bang to a Big Rip singularity. The 
procedure of reconstruction of the model starting by a specif evolution is carried 
out in detail. Then, the dynamical phase space of two of such models is studied, 
analyzing the dynamical classes included in it by varying the initial conditions 
on the fields. 

2.1. Introduction 

As we have seen in the previous section, the discovery of cosmic acceleration 
j6l [7] has stimulated the construction of a class of dark energy models |14] - 
[19] . [20],[21j, [22j-[32j describing this effect. In what follows we shall deal with 
cosmological models based on scalar fields, i.e. the matter content of the universe 
is modeled by means of scalar fields. Notice that the cosmological models based 
on scalar fields were considered long before the observational discovery of cosmic 
acceleration [33]-|35]. 

According to some authors, the analysis of observations permits the existence 
of the moment when the universe changes the value of the parameter w from 
w > —1 to w < —1 [46j-[52]. This transition is called "crossing of the phantom 
divide line". The most recent investigations have shown that the phantom divide 
line crossing is still not excluded by the data |53|-[55|. 

It is easy to see that the standard minimally coupled scalar field cannot give 
rise to the phantom dark energy, because in this model the absolute value of 
energy density is always greater than that of pressure, i.e. \w\ < 1. A possible 
way out of this situation is the consideration of the scalar field models with the 
negative kinetic term, i.e. phantom field models, as I will briefly review in a 
moment. Thus, the important problem arising in connection with the phantom 
energy is the crossing of the phantom divide line. The general belief is that while 
this crossing is not admissible in simple minimally coupled models its explanation 
requires more complicated models such as "multifield" ones or models with non- 
minimal coupling between scalar field and gravity (see e.g. |61)-|66]). 

Some authors of [T] described the phenomenon of the change of sign of the 
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2. Two-field cosmological models 



kinetic term of the scalar field implied by the Einstein equations [67J EH]- It 
was shown that such a change is possible only when the potential of the scalar 
field possesses some cusps and, moreover, for some very special initial conditions 
on the time derivatives and values of the considered scalar field approaching 
to the phantom divide line. At the same time, two-field models including one 
standard scalar field and a phantom field can describe the phenomenon of the (de- 
)phantomization under very general conditions and using rather simple potentials 

PJ-IZZ!. 

In the paper under consideration fl] we have focused to the drastic difference 
between two- and one-field models. 

The reconstruction procedure of the (single) scalar field potential models is 
well-known Pj,t74ti75j,t76j,t77j,t78,i79j,t80jand IMJ-IMJ- Let me recapitulate 
it briefly. 

If the matter is represented by a spatially homogeneous minimally coupled 
scalar field, then the energy density and the pressure are given by the formulae 

e = l02^V(0), (2.1) 

P=l^'- ' (2-2) 
where V{(l)) is a scalar field potential. Friedmann equations thus read 

h' = e='^^' + Vi<P) , (2.3) 
• dV 

+ 3/i0 + — = , (2.4) 
d0 

Where the second is clearly the Klein-Gordon equation on a FRLW background. 
Combining these equations we have 

V = l + h\ (2.5) 

and 

^' = -lh. (2.6) 



Equation (2.5) represents the potential as a function of time t. Integrating 



equation (2.6) one can find the scalar field as a function of time. Inverting this 



dependence we can obtain the time parameter as a function of and substituting 



the corresponding formula into equation (2.5) one arrives to the uniquely recon- 
structed potential V{(j)). It is necessary to stress that this potential reproduces 
a given cosmological evolution only for some special choice of initial conditions 
on the scalar field and its time derivative |67| 176] . Below I shall show that in 
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2.1. Introduction 



the case of two scalar fields one has an enormous freedom in the choice of the 
two-field potential providing the same cosmological evolution. This freedom is 
connected with the fact that the kinetic term has now two contributions. 



Notice that equation (2.6) immediately implies one thing: in order to have (at 
least a phase of) super- acceleration (which, I recall, is defined by /i > 0) with an 
energy density produced by a scalar field, one has to have a negative sign of the 
kinetic term. This is even clearer if one take the specific Hubble parameter 



h{t) 



3{w + l)t 



(2.7) 



i.e. the Hubble parameter for a fluid with w = constant, and tries to reconstruct 



the scalar field model. Integrating (2.6) one has 



0(t) = ±, 



9{w + 1 



Int 



(2. 



which clearly implies w > —1. Had we chose a scalar field with a negative kinetic 
term, we would have been able to reproduce the w < —1 case. 

In order to examine the problem of phantom divide line crossing we shall be 
interested in the case of one standard scalar field and one phantom field ^ , 
whose kinetic term has a negative sign. In this case the total energy density and 
pressure will be given by 



1 • 



1 • 



P 



^''-le-v{<i>,o- 



(2.9) 
(2.10) 



The relation (2.5) expressing the potential as a function of t does not change in 



form, but instead of equation (2.6) we have 



3 



(2.11) 



Now, one has rather a wide freedom in the choice of the time dependence of one of 
two fields. After that the time dependence of the second field can be found from 



equation (2.11). However, the freedom is not yet exhausted. Indeed, having two 
representations for the time parameter t as a function of or ^, one can construct 



an infinite number of potentials V(0, ^) using the formula (2.5) and some rather 
loose consistency conditions. It is rather difficult to characterize all the family of 
possible two-field potentials, reproducing given cosmological evolution h{t). In 
the following, I describe some general principles of construction of such potentials 
and then consider some concrete examples. 
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2. Two-field cosmological models 



2.2. The system of equations for two-field 
cosmological models 



The system of equations, which we study contains (2.5) and (2.11) and two 
Klein-Gordon equations 



(2.12) 
(2.13) 



From equations (2.12) and (2.13) we can find the partial derivatives ^^j^'''-' and 
^^^'^'^^ as functions of time t. The consistency relation 



(2.14) 



is respected. 

Before starting the construction of potentials for particular cosmological evo- 
lutions, it is useful to consider some mathematical aspects of the problem of 
reconstruction of a function of two variables in general terms. 



2.2.1. Reconstruction of the function of two variables, 
which in turn depends on a third parameter 

Let us consider the function of two variables F{x, y) defined on a curve, parame- 
terized by t. Suppose that we know the function F{t) and its partial derivatives 
as functions of t: 

F{x{t),y{t)) = F{t) 

dF{x{t),y{t)) _ OF 
dx 
_ dF , 
dy 



dx 

dF{x{t),y{t)) 



dy 



it) , 
it) . 



These three functions should satisfy the consistency relation 

m 

As a simple example we can consider the curve 

x{t) = t, y{t) = e , 



dF,^. dF,^. 



(2.15) 
(2.16) 

(2.17) 



(2.18) 



(2.19) 



20 



2.2. The system of equations for two-field cosmological models 



while 



F{t) = , 
dF 
dx 
dF 
dy 



t 



(2.20) 
(2.21) 

(2.22) 



and equation (2.18) is satisfied. 



Thus, we would like to reconstruct the function F{x, y) having explicit expres- 



sions in right-hand side of equations (2.15)- (2.17). This reconstruction is not 



unique. We shall begin the reconstruction process taking such simple ansatzes 
as 



F^{x,y) = G^{x) + HM 
F,{x,y) = G2{x)H2{y) , 
Fs{x,y) = {Gs{x) + H3{y)y 



The assumption (2.23) immediately implies 

dFi _ dCi 
dx dx ' 

dFi _ dHi 
dy dy 



Therefrom one obtains 



Gi{x) 



dFi 



dx 



y(t{x'))dx' 



dHi 



dy 



Hence 



Fi{x,y) 



dF, r^H, , 

y{t{x))dx + / -—{t[y))dy 



dx 



dy' 



For an example given by equations (2.19)-(2.22) the function Fi{x,y) is 
Fi{x,y)= / t^{x')dx' + / t{y')dy' = / x' W ± 



y'dy' 



Explicitly 



Fi{x,y) 



- {x^ + 2y3/2) , X > , 
\ (x^ - 2y3/2) , X < . 



(2.23) 
(2.24) 
(2.25) 

(2.26) 
(2.27) 

(2.28) 
(2.29) 

(2.30) 
(2.31) 

(2.32) 
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Similar reasonings give for the assumptions (2.24), and (2.25) correspondingly 

1 dF^ 



F2{x,y) =exp|y" 



F3{x,y) 



dx 
a 



F dx 

L.idF 
dx 



{t{x))dx + 
it{x)) + 



1 OF 
F~dy 



{t{y))dy 



dy 



a 



dy 



(2.33) 



(2.34) 



For our simple example (2.19)-(2.22) the functions F2{x,y) and F-,,{x^y) have 
the form 

F2{x,y) = xy , 



^3(3;, y) 



(a;3/a(_)3/a2^3/2a) 



a; > , 
X < . 



(2.35) 
(2.36) 



Thus, we have seen that the same input of "time" functions (2.20)-(2.22) on the 
curve (2.19) produces quite different functions of variables x and y. 

Naturally, one can introduce many other assumptions for reconstruction of 
F[x,y). For example, one can consider linear combinations of x and y as func- 
tions of the parameter t and decompose the presumed function F as a sum or a 
product of the functions of these new variables. 

Now I present a way of constructing the whole family of solutions starting 
from a given one. Let us suppose that we have a function FQ{x,y) satisfying all 
the necessary conditions. Let us take an arbitrary function 



f{x,y) = f 



t{y) 



which depends only on the ratio t{x)/t{y). We require also 

/(i) = i, r(i) = o. 



(2.37) 



(2.38) 



i.e. the function reduces to unity and its derivative vanishes on the curve 
{x(t),y{t)). Then it is obvious that the function 



F{x,y) = Fo{x,y)f 



t{y) 



(2.39) 



is also a solution. This permits us to generate a whole family of solutions, 
depending on a choice of the function /. Moreover, one can construct other 
solutions, adding to the function F{x, y) a term proportional to {t{x) — t{y)Y. 
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2.2. The system of equations for two-field cosmological models 



2.2.2. Cosmological applications, an evolution "Bang to 

RFI 
ip 

To show how this procedure works in Cosmology, we consider a relatively simple 
cosmological evolution, which nevertheless is of particular interest, because it 
describes the phantom divide line crossing. Let us suppose that the Hubble 
variable for this evolution behaves as 

where A is a positive constant. At the beginning of the cosmological evolution, 
when t — )■ the universe is born from the standard Big Bang type cosmological 
singularity, because h{t) ~ 1/t. Then, when t — )■ tji, the universe is super- 
accelerated, approaching the Big Rip singularity h{t) ~ — t). Substituting 
the function (2.40) and its time derivative into equations. ( |2.11[ ) and (2.5) we 
come to 

A(2t -tn + SA) , , 

For convenience let me choose also the parameter A as 

tE 

3 



A = ^ , (2.43) 



Then 



hit) = T^TTT^ , (2.44) 



3t{tR - 1) 

and 

Let us consider now a special choice of functions (j){t) and ^{t) used alreadjj^ 
in [69]-[72j, 

0(t) = -- arctanh ~ ^ , (2.46) 
3 V ^-R 

e(t) = ^ln-^. (2.47) 

6 tfi — t 



^Notice that the origin of two scalar fields has been associated in [M'-f75] with a non-hermitian 
complex scalar field theory and there a classical solution was found as a saddle point in 
"double" complexification. I will discuss in more detail the link between non-hermiticity 
and phantom fields in chapter |4] 
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2. Two-field cosmological models 



The derivatives of the potential with respect to the fields (f) and ^ could be found 
from the Klein-Gordon equations ( 2.12[ ) and (2.13): 



We can obtain also the time parameter as a function of or ^: 

*<^' = cosh-(-3^/4) • P'^") 

" exp(-3?/y2) + 1 ■ '^'^^^ 

Now we can make a hypothesis about the structure of the potential V{^, (p): 

V2{iA)=G{0H{^) ■ (2.52) 
Applying the technique described in the subsection A, we can get G{^): 

lnG(e) = J (^^) (2.53) 

2tR 3 tmn-my ^ 

and 

G(0 = exp(3y20 . (2.54) 

To find H{(f)) one can use the analogous direct integration, but we preferred to 
implement a formula 

„,.., Vim ,„ 



which gives 



and hence, 



Finally, 



Of 

Hm) = ^ , (2.56) 



i/(0) = -^cosh6(-3(/)/4). (2.57) 
^2(e, 0) = ^ cosh6(-30/4) exp(3y20 . (2.58) 
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2.2. The system of equations for two-field cosmological models 



Here we have reproduced the potential studied in [69j-|72|. 
Making the choice 

^i(e,0) = G(o + i^(0), 

we derive 



(2.59) 



3tl 



Vi(e,0 

sinh^(-30/4) + smh\-34>/4) - 1 



3 3 



+ lnsinh^(-30/4) 



(2.60) 



sinh^(-30/4) 

Now, we can make another choice of the field functions (f){t) and C,{t), satisfying 

V2 



the condition (2.41): 



,/N v2, t 

3 V ^-R 

The time parameter t is a function of fields is 

tR 



m 



exp (-30/^2) + 1 
tanh2(3e/4) ' 



(2.61) 
(2.62) 

(2.63) 
(2.64) 



Looking for the potential as a sum of functions of two fields as in equation 



(|2.59|) after lengthy but straightforward calculations one comes to the following 

1 + I sinh^ (3^/4) + 3 sinh^ (3^/4) 



potential: 



3ti 



3sinh^(35/4) 

+21nsinh2(3^/4) + | exp(-30/v^) - 3^20 

-2exp(30/y2) - |exp(0/y2)] . (2.65) 

Similarly for the potential designed as a product of functions of two fields 



(2.52) we obtain 



sinh2(3^/4) cosh2(3^/4) exp(-3v^0) 



(2.66) 



One can make also other choices of functions (j){t) and ^(t) generating other 
potentials, but I shall not do it here, concentrating instead on the qualitative 
and numerical analysis of two toy cosmological models described by potentials 
and ([2:58|. 
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2.3. Analysis of cosmological models 

It is well known |85] that for the qualitative analysis of the system of cosmological 
equations it is convenient to present it as a dynamical system, i.e. a system of 
first-order differential equations. Introducing the new variables x and y we can 
write 

^6 = X , 



X 



y 



—3 sign(/i) X 
^ y = -3 sign(/i) y 
Notice that the reflection 



i-4 + v{^,<p 



dV 

dcj> 

dV 
95 



(2.67) 



X — )■ —X , 

y ^ -y , 
t -t 



(2.68) 



transforms the system into one describing the cosmological evolution with the 
opposite sign of the Hubble parameter. The stationary points of the system 



(2.67) are given by 



X 



0, y = 



dV 

90 



dV 



(2.69) 



2.3.1. Model I 



In this subsection I shall analyze the cosmological model with two fields - stan- 
dard scalar and phantom - described by the potential (2.65). For this potential 
the system of equations (2.67) reads 



= X , 

i= y , 

x= -3 sign(/i) x^f + sinh2(3e/4) cosh'^(3e/4)e-3v^'^ 

+ |^sinh'(3e/4)cosh6(3e/4)e-3^^ (^.70) 
y= -3 sign(/i) y^f ~ ? + 9^ sinh2(3e/4) cosh%3^/A)e-^^<t> 

^ .inh(3€/4)cosh^(3€/4) ^1 ^03^=^(3^/4) + smh\3^/A)] C'^v^-^ . 
k R 

It is easy to see that there are stationary points 

^ = cj)o,^ = 0,x = 0,y = 0, (2.71) 
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where (pQ is arbitrary. For these points the potential and hence the Hubble 
variable vanish. Thus, we have a static cosmological solution. We should study 



the behavior of our system in the neighborhood of the point (2.71) in linear 
approximation: 



= 


X , 


i = 


y , 


X = 


, 


y = 





(2.72) 



-3\/2(/)o 



One sees that the dynamics of in this approximation is frozen and hence we 
can focus on the study of the dynamics of the variables ^, y. The eigenvalues of 
the corresponding subsystem of two equations are 



Ai,: 



-300 



2tR 



(2.73) 



These eigenvalues are real and have opposite signs, so one has a saddle point in 



the plane {^,y) and this means that the points (2.71) are unstable. 

One can make another qualitative observation. Freezing the dynamics of ^ 
independently of 0, namely choosing y = 0,C, =0, which implies also ^ = ^ = 0, 
one has the following equation of motion for (p: 



(f) + 3h(j) = . 



(2.74) 



Equation (2.74) is nothing but the Klein-Gordon equation for a massless scalar 



field on the Friedmann background, whose solution is 

V2. 



Int , 



and which gives a Hubble variable 



h{t) 



1 

3t 



(2.75) 



(2.76) 



This is an evolution of the fiat Friedmann universe, filled with stiff matter with 
the equation of state p = e. It describes a universe born from the Big Bang 
singularity and infinitely expanding. Naturally, for the opposite sign of the 
Hubble parameter, one has the contracting universe ending in the Big Crunch 
cosmological singularity. 

Now, I describe some results of numerical calculations to have an idea about 
the structure of the set of possible cosmological evolutions coexisting in the 
model under consideration. We have carried out two kinds of simulations. First, 
we have considered neighborhood of the plane y, ^ with the initial conditions on 
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I III -4 I- 1 

Figure 2.1.: An example of section of the 4d phase space obtained with numerical 
calculations. We can see a series of trajectories corresponding to 
different choices of the initial conditions. Every initial condition in 
the graphic is chosen on the dashed "ellipse" centered in the origin 
and is emphasized by a colored dot. A "saddle point like" structure 
of the set of the trajectories clearly emerges. 



the field such that 0(0) = 0, 0(0) = (see figure 2.1 ). The initial conditions for 
the phantom field were chosen in such a way that the sum of absolute values of 
the kinetic and potential energies were fixed. Then, running the time back and 
forward we have seen that the absolute majority of the cosmological evolutions 
began at the singularity of the "anti-Big Rip" type (figure 2.2). Namely, the 
initial cosmological singularity were characterized by an infinite value of the 
cosmological radius and an infinite negative value of its time derivative (and 
also of the Hubble variable). Then the universe squeezes, being dominated by 
the phantom scalar field C,- At some moment the universe passes the phantom 
divide line and the universe continues squeezing but with h < 0. Then it achieves 
the minimal value of the cosmological radius and an expansion begins. At some 
moment the universe undergoes the second phantom divide line crossing and its 
expansion becomes super-accelerated culminating in an encounter with a Big 
Rip singularity. Apparently this scenario is very different from the standard 
cosmological scenario and from its phantom version Bang-to-Rip, which has 
played a role of an input in the construction of our potentials. The second 
procedure, which we have used is the consideration of trajectories close to our 
initial trajectory of the Bang-to-Rip type. The numerical analysis shows that 
this trajectory is unstable and the neighboring trajectories again have anti-Big 
Rip - double crossing of the phantom line - Big Rip behavior described above. 
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1.0 



0.5 



TOO 



-80 



-60 



-40 



-20 



-0.5 



-1.0 



Figure 2.2.: Typical behavior of the Hubble parameter for the model I trajecto- 
ries. The presence of two stationary points (namely a maximum fol- 
lowed by a minimum) indicates the double crossing of the phantom 
divide line. Both the initial and final singularities are characterized 
by a Big Rip behavior, the first is a contraction and the second is 
an expansion. 
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However, it is necessary to emphasize that a small subset of the trajectories of 
the Bang-to-Rip type exist, being not in the vicinity of our initial trajectory. 



2.3.2. Model II 

In this subsection we shall study the cosmological model with the potential 



(2.58). Now the system of equations (2.67) looks like 



X 



X , 

y , 

-3 sign(/i) - ^ + ^ 



2 



cosh'^(30/4)exp{3y2^} 
-Ticosh^(30/4)sinh(30/4)e3^^ , 



(2.77) 



-3sign{h)y 



+ 



2V2 



5^ _ _L ^ 
2 2 9*2 



cosh^(30/4) exp{3\/2C} 



cosh^(30/4)e3v^« . 



Notice that the potential (2.58) has an additional reflection symmetry 

l^2(e,0)= ^2(^,-0) . 



(2.78) 



This provides the symmetry with respect to the origin in the plane {4>,x). The 



system (2.77) has no stationary points. However, there is an interesting point 

= 0, X = , (2.79) 

which freezes the dynamics of and hence, permits to consider independently 
the dynamics of C, and y, described by the subsystem 



e= y, 

y= -3 sign(/i) y 



2 



I 2 3v^C I 2v|p3V2C 



(2.80) 



9t 



3tf 



Apparently, the evolution of the universe is driven now by the phantom field and 
is subject to super-acceleration. 

In this case the qualitative analysis of the differential equations for ^ and ?/, 
confirmed by the numerical simulations gives a predictable result: being de- 
termined by the only phantom scalar field the cosmological evolution is char- 
acterized by the growing positive value of h. Namely, the universe begins its 
evolution from the anti-Big Rip singularity {h = —00) then h is growing passing 
at some moment of time the value h = (the point of minimal contraction of 
the cosmological radius a{t)) and then expands ending its evolution in the Big 
Rip cosmological singularity {h = +00). 
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Another numerical simulation can be done by fixing initial conditions for the 



phantom field as ^(0) = 0,y{0) = (see figure 2.3). Choosing various values of 
the initial conditions for the scalar field 0(O),a:(O) around the point of freezing 
= 0, X = we found two types of cosmological trajectories: 

1. The trajectories starting from the anti-Big Rip singularity and ending in the 
Big Rip after the double crossing of the phantom divide line. These trajectories 
are similar to those discussed in the preceding subsection for the model I. 

2. The evolutions of the type Bang-to-Rip. 

Then we have carried out the numerical simulations of cosmological evolutions, 
choosing the initial conditions around the point of the phantomization point with 
the coordinates 

m = , 
m = f . 

4 1 
^(0) = -arctanh^^ , 
^ 3 72 ' 

1/(0) = ^ . (2.81) 

This analysis shows that in contrast to the model I, here the standard phan- 
tomization trajectory is stable and the trajectories of the type Bang-to-Rip are 
not exceptional, though less probable then those of the type anti-Big Rip to Big 
Rip. 



2.4. Conclusions 

In the paper [T] we have considered the problem of reconstruction of the potential 
in a theory with two scalar fields (one standard and one phantom) starting with 
a given cosmological evolution. It is known ( see e.g. |73| 176] ) that in the case 
of the only scalar field this potential is determined uniquely as well as the initial 
conditions for the scalar field, reproducing the given cosmological evolution. 
Changing the initial conditions, one can find a variety of cosmological evolutions, 
sometimes qualitatively different from the "input" one (see e.g. [76j). In the case 
of two fields the procedure of reconstruction becomes much more involved. As 
we have shown here, there is a huge variety of different potentials reproducing 
the given cosmological evolution (a very simple one in the case, which we have 
explicitly studied here). Every potential entails different cosmological evolutions, 
depending on the initial conditions. 

It is interesting that the existence of different dynamics of scalar fields corre- 
sponding to the same evolution of the Hubble parameter h{t) can imply some 
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X 

20 - 
15 - 




-20 I- 

Figure 2.3.: An example of section of the 4D phase space for the model II ob- 
tained with numerical calculations. We can see a family of different 
trajectories corresponding to different choices of the initial condi- 
tions. Every initial condition in the graphic is chosen on the dashed 
"ellipse". The origin represents the point of freezing. 

observable consequences connected with the possible interactions of the scalar 
fields with other fields. In this case, the time dependence of the scalar fields con- 
sidered above can directly affect physically observable quantities. Indeed, some 
results in this direction are presented in chapter [3} devoted to the explanation 
of the paper ||2j. 
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3. Two-field models and cosmic 
magnetic fields 



The present chapter is intended to present the content of the paper J^. Starting 
from the models studied in the preceding chapter, we couple one of the fields to a 
cosmic electromagnetic field, whose existence is experimentally confirmed. The 
idea is to try and find a discrimination between cosmological models presenting 
the same background evolution, by means of the effects of the coupling with other 
(observable) fields. 



3.1. Introduction to the content of the paper p] 

In the preceding chapter I have shown the procedure of reconstruction of the 
potential in two-field models. It was shown that there exists a huge variety of 
potentials and time dependences of the fields realizing the same cosmological 
evolution. Some concrete examples were considered, corresponding to the evo- 
lution beginning with the standard Big Bang singularity and ending in the Big 
Rip singularity [571 EH] • 

One can ask oneself: what is the sense of studying different potentials and 
scalar field dynamics if they imply the same cosmological evolution? The point 
is that the scalar and phantom field can interact with other fields and influence 
not only the global cosmological evolution but also other observable quantities. 

One of the possible effects of the presence of normal and phantom fields could 
be their influence on the dynamics of cosmic magnetic fields. The problem 
of the origin and of possible amplification of cosmic magnetic fields is widely 
discussed in the literature [861 ISTf 155]. In particular, the origin of such fields can 
be attributed to primordial quantum fluctuations |89l l90l [91] and their further 
evolution can be influenced by hypothetic interaction with pseudo-scalar fields 
breaking the conformal invariance of the electromagnetic field [92] -[96], |97| . 
In the paper under consideration we analyzed the evolution of magnetic fields 
created as a result of quantum fluctuations, undergoing the inflationary period 
with unbroken conformal invariance and beginning the interaction with pseudo- 
scalar or pseudo-phantom fields after exiting the inflation and entering the Big 
Bang expansion stage, which is a part of the Bang-to-Rip scenario described in 
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the preceding chapter. We used different field reahzations of this scenario and 
we shall see how the dynamics of the field with negative parity influences the 
dynamics of cosmic magnetic fields. 

What follows is so organized: in section 3.2 I very briefly recall the Bang-to- 
Rip scenario and the two models that will be considered (the same of chapter |2| 
in section 3.3 I introduce the interaction of the fields (phantom or normal) with 



an electromagnetic field and write down the corresponding equations of motion; 
in section 3.4 I describe the numerical simulations of the evolution of magnetic 
fields and present the results of these simulations; the last section 3.5 is devoted 
to concluding remarks. 



3.2. Cosmological evolution and (pseudo)-scalar 
fields 



I shall consider a spatially flat Friedmann universe with the FLRW metric (1.2) 



Notice that the physical distance is obtained by multiplying dl by the cosmo- 
logical radius a{t). We would like to consider the cosmological evolution char- 
acterized by the following time dependence of the Hubble variable h{t), where, 
I recall, "dot" denotes the differentiation with respect to the cosmic time ^ 



hit) 



tR 



3t{tn - 1) 



(3.1) 



I called this kind of scenario Bang-to-Rip in section [2.2.2[ at small values of t the 
universe expands according to power law: a ~ t^^^ while at t — )■ t^j the Hubble 
variable explodes and one encounters the typical Big Rip type singularity. (The 



factor one third in (3.1) was chosen for calculation simplicity). 

In the paper we are considering [2j we kind of continued the work presented in 
the preceding chapter, or equivalently in [T] . Specifically we considered precisely 
the two models of the preceding chapter (see sections 2.3.1, 2.3.2) and coupled 
them with a magnetic field. 



We have seen in chapter [2] that, analyzing the Friedmann equation (1.8a 



for two-field models there is huge variety of potentials V^(0, realizing a given 



^Cfr. equation <\2M ]. 

use here, in accordance with [TJ|2], the following system of units h = l,c 



1 and 8nG = 3. 

In this system the Planck mass m p , the Planck length I p and the Planck time t p are equal 
to 1. Then when we need to make the transition to the "normal", say, cgs units, we should 
simply express the Planck units in terms of the cgs units. In all that follows we tacitly 
assume that all our units are normalized by the proper Planck units. Thus, the scalar field 
entering as an argument into the dimensionless expressions should be divided by the factor 

y/m,p/tp. 
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3.2. Cosmological evolution and (pseudo)-scalar fields 



evolution, in contrast to models with one scalar field. Moreover, besides the 
freedom in the choice of the potential, one can choose different dynamics of the 
fields and ^(t) realizing the given evolution. For simplicity I repeat here 



the main ingredients of the two models of 2.3.1 , 2.3.2: The first potential is (cfr. 
with (|2]66|) 

^^(^' = ^ cosh6(-30/4) exp(3y20 . (3.2) 



and the fields 



(hit) = — arctanh^ / -^^ 



(3.3) 



(3.4) 



(The expression for the potential should be multiplied by the factor mp/{clp), 
while the expressions for the fields (j){t) and ^{t) should be multiplied by ^Jmp/tp. 
For the relation between Planck units and cgs ones see e.g. |98]). If we would 
like to substitute one of these two fields by the pseudo-scalar field, conserving 
the correct parity of the potential, we can choose only the field because the 
potential V/ is even with respect to 0, but not with respect to ^. In what follows 



I shall call the model with the potential (3.2), the pseudo-scalar field (3.3) and 



the scalar phantom (3.4) model / 



Consider another potential (cfr. with (2.58)) 



with the fields 



9^1 



sinh^(3^/4) cosh2(3^/4) exp(-3y20) . 



(3.5) 



v^, t 
In 

3 tR-t 



(3.6) 



m = iarctanh,/! 
3 V ^-R 



(3.7) 



This potential is even with respect to the field ^. Hence the model II is based 



on the potential (3.5), the pseudo-phantom field (3.7) and the scalar field (3.6) 



They will be the fields with the negative parity which couple to the magnetic 
field. 
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3.3. Post-inflationary evolution of a magnetic 
field interacting with a pseudo-scalar or 
pseudo-phantom fields 

The action of an electromagnetic field interacting with a pseudo-scalar or pseudo- 
phantom field 4> is 

S = -\j d'x^iF^^F^" + a(j)F,,Fn , (3.8) 

where a is an interaction constant and the dual electromagnetic tensor F^'^ is 
defined as 

F^^ = ^E'"'P''Fp^ , (3.9) 

where 

with the standard Levi-Civita symbol 



e 



^pupa], e0123 = +1 • (3-11) 



Variating the action (3.8) with respect to the field we obtain the field 
equations 

Vf.F^" = -ad^cPF^'' , (3.12) 

Vf.F"" = . (3.13) 
The Klein-Gordon equation for the pseudo-scalar field is 

V^V^^ = -aFpuF^^ . (3.14) 

The Klein-Gordon equation for the pseudo-phantom field (which is the one that 



couples with the magnetic field in the model //) differs from equation (3.14) by 
change of sign in front of the kinetic term. In what follows we shall neglect the 
influence of magnetic fields on the cosmological evolution, i.e. we will discard 



the electromagnetic coupling in equation (3.14). 

If one wants to rewrite these formulae in terms of the three-dimensional quan- 
tities (i.e. the electric and magnetic fields) one can find the expression of the 
electromagnetic tensor in a generic curved background, starting from a locally 
flat reference frame — in which it is well known the relation between electro- 
magnetic fields and F — and using a coordinate transformation. It is easy to 
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3.3. Interaction between magnetic and pseudo-scalar /phantom field 



see that we have, for the metric ((1.2)): 










-aEi 


-aE2 


-aEs \ 


1 




aEi 





B3 


-B2 






aE2 







Bi 




v 


aEi 


B2 




/ 



(3.15) 



The field equations (3.12), (3.13) and (3.14) rewritten in terms of E and B 
become 

V ■ E = -aVcp- B . (3.16a) 



^0(0^^) - V X (aB) = -a[do(p{a^B) - V0 x (aE)] , 
do{a^B) - V X (aE) = , 
V-5 = . 



(3.16b) 
(3.16c) 
(3.16d) 



For a spatially homogeneous pseudo-scalar field equations (3.16a) and (3.16b) 
look like 

V-^ = 0, (3.17a) 



doia^E) - V X (aB) = -ado^ia^B) . 



(3.17b) 



Taking the curl of (3.17b) and substituting into it the value of E from (3.16c) 
we obtain 



d^ia^B) + h{t)do{a^B) - ^^'^^f^^ - -9o0V x {a^B) = 



(3.18) 



where A*^^) stands for the three-dimensional Euclidean Laplacian operator. 
Let me introduce 

F{x,t) = a^{t)B{x,t) (3.19) 



and its Fourier transform 

Fik,t) 



(27r)3/2 



(3.20) 



Here the field B is an observable magnetic field entering into the expression for 
the Lorentz force. The field equation for F{k,t) is 



F{k,t) + h{t)F(k,t) + 



%a 



(j)kx 



F{k,t)=0, 



(3.21) 
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where "dot" means the time derivative. This last equation can be further sim- 
plified: assuming k = (fc, 0, 0) and defining the functions F± = {F2 ± i-F3)/\/2 
one arrives to 

F± = , (3.22) 

where I have omitted the arguments k and t. 

Assuming that the electromagnetic field has a quantum origin (as all the fields 
in the cosmology of the early universe [99]) the modes of this field are represented 
by harmonic oscillators. Considering their vacuum fluctuations responsible for 
their birth we can neglect the small breakdown of the conformal symmetry and 
treat them as free. In conformal coordinates {1], x) such that the Friedmann 
metric has the form 

ds^ = a^{ri){dr]'^ - dl^) , (3.23) 

the electromagnetic potential Ai with the gauge choice Aq = 0, djA^ = satisfies 
the standard harmonic oscillator equation of motion 

Ai + k'^Ai = . (3.24) 

Hence the initial amplitude of the field A^ behaves as Ai = l/\/2k, while the 
initial amplitude of the functions F is ^Jkjl. The evolution of the field F 
during the inflationary period was described in [;97] , where it was shown that the 
growing solution at the end of inflation is amplified by some factor depending 
on the intensity of the interaction between the pseudo-scalar field and magnetic 
field. 

Here we are interested in the evolution of the magnetic field interacting with 
the pseudo-scalar field after inflation, where our hypothetic Bang-to-Rip sce- 
nario takes place. More precisely, I would like to see how different types of 
scalar-pseudo-scalar potentials and field dynamics providing the same cosmo- 
logical evolution could be distinguished by their influence on the evolution of 
the magnetic field. We do not take into account the breaking of the conformal 
invariance during the inflationary stage and all the effects connected with this 
breakdown will be revealed only after the end of inflation and the beginning of 
the Bang-to-Rip evolution. This beginning is such that the value of the Hubble 
parameter, characterizing this evolution is equal to that of the inflation, i.e. 

^(^O) = il"" , . = /^inflation ^ lO'^^S"^ . (3.25) 

In turn, this implies that we begin evolution at the time moment of the order of 



F± + hF± + 



± a- 



38 



3.4. Generation of magnetic fields: numerical results 



I shall consider both the components and F_ 



scenarios / and // described in the preceding section 3.2 



and we shall dwell on the 
Anyway, our assump- 



tion regarding the initial conditions for equation (3.22) can be easily modified in 
order to account for the previous possible amplification of primordial magnetic 
fields as was discussed by |97]. Thus, all estimates for the numerical values of the 
magnetic fields in today's universe should be multiplied by some factor corre- 
sponding to the amplification of the magnetic field during the infiationary stage. 
Hence, our results refer more to differences between various models of a post 
infiationary evolution than of the real present values of magnetic fields, whose 
amplification might be also combined effect of different mechanism [86l |87| |88]. 

j-i96j, uniEniET], m 



3.4. Generation of magnetic fields: numerical 
results 



In this section I present the results of numerical simulations for the two models 



/ and //. In these models, the equations of motion for the modes F±{k, t) (3.22) 



reai 



tR-t 
t 



2/3 



± a 



tn-t 



1/3 



k<l> 



(3.26) 



where $ stands for the scalar field 
the model //, so that 



in the model / and for the phantom in 



tR 



II 



^R 



3 ^t{tR - t) ■ 



(3.27) 



Equation (3.26) is solved for different values of the wave number k and the 
coupling parameter a. (The parameter a has the dimensionality inverse with 
respect to that of the scalar field; the wave number k has the dimensionality of 
inverse length; the time = lO^^s). Qualitatively let me remark that in (3.22) 



the coupling term infiuence becomes negligible after some critical period. After 
that the magnetic fields in our different scenarios evolve as if the parameter 



a in (3.22) had been put equal to zero. Indeed, it can be easily seen that the 



interaction term vanishes with the growth of the cosmological radius a. Then the 



^The reader can easily verify that this equation is obtained imposing the normaHzation 
« I today = Ij where today-time is taken to be near the crossing of the phantom divide 
line, i.e. at t ~ This implies in turn that at the beginning of the "Bang-to- Rip" 

evolution the cosmological radius is a{to) ~ 10^"'^^. 
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distinction between the two models is to be searched in the early time behavior 
of the field evolution. 

Noting that in both our models / and // the time derivative $ is positive, by 



inspection of the linear term in equation (3.22) we expect the amplification to be 
mainly given for the mode F_ provided the positive sign for a is chosen; so we 
will restrict our attention on F_. We can also argue that the relative strength 



of the last two terms in the left-hand side of (3.26) is crucial for determining 
the behavior of the solution: when the coupling term prevails (I remark that 
we are talking about F_ so this term is negative in our models) then we expect 
an amplification, while when the first term dominates we expect an oscillatory 
behavior. For future reference it is convenient to define 

k k ftn-t^^'^ 



A{t-<^)^—^ = ^[^\ , (3.28) 
which is just the ratio between the last two terms in the left-hand side of equa- 



tion ( |3.26[ ). 

Indeed our numerical simulations confirm these predictions. Let us consider 
the model I with a = l(tp/mp)^/^ and k = 10~^^/p^, where Ip is the Planck 
length. Such a value of the wave number k corresponds to the wave length of 
Ikpc at the present moment. We obtain an early-time amplification of about 
2 orders of magnitude, with the subsequent oscillatory decay. Notice that the 
parameter A in this model at the beginning is very small: this corresponds to 
the dominance of the term proportional to $ and, hence, to the amplification 
of the field F_. At the time scale of the order of lO^^tp, where tp is the Planck 
time, this regime turns to that with big values of A where the infiuence of the 
term proportional to $ is negligible. 

For the same choice of the parameters a and k in the model // the amplifica- 
tion is absent. 



In figure 3.1 you see the time dependence of the function A for the mod el / 
for the values of a = l(tp/mpY^'^ and k = 10~^^/p^ chosen above. The figure 3.2 
manifests the amplification of the magnetic field in the model /. 

Naturally the effect of amplification of the magnetic field grows with the cou- 



pling constant a and diminishes when the wave number k increases. In figure [373 
are displayed the results for the case of a = 100(tp/mp)^/^, which is admittedly 
extreme and possibly non realistic, but good for illustrative purposes. Here the 
amplification is more evident and extends for a longer time period. 

In [2], we also tried to make some estimates of the cosmic magnetic fields in 
the universe today, using the correlation functions. The correlation function for 
the variable F is defined as the quantum vacuum average 

G,,{t,x-y) = {0\Fi{t,x)F,{t,y)\0) (3.29) 
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3.4. Generation of magnetic fields: numerical results 




Figure 3.1.: Plot of the ratio A for the model /, with the parameter choice 
a = l^tp/mpY^"^, k = 10~^^lp^. It can be easily seen that at a time 
scale of order lO^^tp the ratio becomes greater than 1. 




Figure 3.2.: Plot of the time evolution of the absolute value of the complex field 
F (given in Planck units) in model / with the parameter choice 
a = l{tp/mpY^'^, k = 10^^^/p^. The behavior, as said above in the 
text, consists in an amplification till a time of order lO^^tp, after 
which the oscillations begin. 
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3. Two-field models and cosmic magnetic fields 



|F| 




Figure 3.3.: Plot of the time evolution of the absolute value of the complex field 
F (given in Planck units) in model / with the parameter choice 
a = 100(tp/mp)^/^, k = 10~^^lp^. The behavior, consists in an 
amplification till a time of order lO^Hp, after which the oscillations 
begin. 



and can be rewritten as 

Gij{t,x-y) 



(27r)3 

Integrating over the angles, we come to 

k'^dk sm{k\x — y\) 



e-iH£-y)F^^t,k)FUt,k). 



(3.30) 



G,j{t,x-y) 



27?^ k\x — y\ 



F,{t,k)F;{t, k). 



(3.31) 



To estimate the integral (|3.3l|) we notice that the main contribution to it comes 

(3.32) 



from the region where k ^ l/\x — y\ (see e.g.jHS]) and it is of order 

1 



L3 



where L = \x — y\. In this estimation the amplification factor is 

F(l/L) 



(3.33) 



where the subscript i is not present since we have taken the trace over polariza- 
tions. 
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3.5. Conclusions 



Now we are in a position to give numerical values for the magnetic fields at 
different scales in the model / for different values of the coupling parameter a. 
These values (see Table 3.1 ) correspond to three values of the coupling parameter 
a (1,10 and 100 {tp/mpY^'^) |^ and to two spatial scales L determined by the 
values of the wave number k. We did not impose some physical restrictions on 
the value of a. It is easy to see that the increase of a implies the growth of the 
value of the magnetic field B. 

Let me stress once again that we ignored the effects of amplification of the 
magnetic fields during inflation to focus on seizable effects during evolution. 



Table 3.1.: The Table displays the values of the magnetic field B 

corresponding to the chosen values of a and k. The length L 
refers to the present moment when a = 1. 





a = ly^tp/mp 


a = lO^ytp/mp 


a = lOO^ytp/nip 


k = 10-^%^ {L = 


Ikpc) 


B ~ lO-'^^G 


B ~ 10-*^°^ 


B ~ lOOG 


k = 10-5°/p^ (L = 


lO^V) 


B ~ lO-^^G 


B ~ lO-^^G 


B ~ W^G 



Finally notice that our quantum "initial" conditions correspond to physical 
magnetic fields which for presented values of k are Bin ~ k'^/a^ is equal to 
10-34G for k = and Bin ~ l^'^^G for k = lQ-^>%\ 



3.5. Conclusions 

We have seen that the evolution of the cosmic magnetic fields interacting with 
a pseudo-scalar (pseudo-phantom) field is quite sensitive to the concrete form of 
the dynamics of this field in two-field models where different scalar field dynamics 
and potentials realize the same cosmological evolution. 

The sensitivity of the evolution of the magnetic field with respect to its he- 
licity is confirmed, given the sign of the coupling constant a and that the $ is 
a monotonic function of time (as it is really so in our models). We gave also 
some numerical estimates of the actual magnetic fields up to the factor of am- 
plification of such fields during the inflationary period. The toy model of the 
Bang-to- Rip evolution studied in this paper [2], cannot be regarded as the only 

'^It is useful to remark that at this length scales values of a less than 1 make the coupling 
with the cosmological evolution negligible. 
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responsible for the amplification of cosmic magnetic fields implying their present 
observable values. It rather complements some other mechanisms acting before. 
However, the difference between cosmic magnetic fields arising in various models 
(giving the same expansion law after the inflation) is essential. It may provide 
a discriminating test for such models. 
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4. Phantom without phantom in 
a PT-symmetric background 



This chapter is devoted to presenting the content of the paper I3j. In nuce, the 
idea is that if we put ourselves in the framework of PT symmetric Quantum 
Theory we can have an effective phantom field which is stable to quantum fluc- 
tuations. "Effective" in the usual sense, i.e. the real field is a standard scalar 
field, hut becomes a phantom field once we 'sit' on a (specific) classical solution. 

Before going into the details of the model presented in J3], let me briefly in- 
troduce the key ideas and concepts of PT symmetric Quantum Mechanics. 

4.1. PT-symmetric Quantum Mechanics: a brief 
introduction 

PT symmetric Quantum Mechanics was born in 1998 in the seminal paper |1UU| 
by Carl Bender and Stefan Boettcher. The motivating idea is very simple: how 
to make sense of non-hermitian hamiltonians. Indeed it had been known for 
long (since late 50 's) that some non-hermitian hamiltonians come out naturally 
in some systems [101J - [1Q5J . but everybody where rather skeptical about their 
physical sense. 

However, since the 80 's, there where hints in the following direction: some 
non-hermitian hamiltonians do have a real and bounded spectrum |1U6[ I1U7| . 
You guess the point: Bender and collaborators rigorously proved since the very 
beginning of the theory in 1998 that it is not necessary to require hermiticity 
in order to have a real and bounded spectrum. We can weaken that axiom 
of QM with the more physical requirement: the hamiltonian operator must be 
PT symmetric. 

A PT transformation is a combined transformation made of parity 




(4.1) 
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and time reversal 




(4.2) 



thus 

PT <j p^p (4.3) 
i 

we can think of it as a spacetime reflection. 

PT symmetric hamiltonians do have real and bounded spectrum, thus they 
have the usual physical interpretation as operators whose eigenvalues are the 
energy levels of the system under investigatioiij^ 

As Gell-Mann said (and as QM teaches us) "everything that is not forbidden is 
compulsory". This is the sort of principle that guided - together with the need of 
dealing with non-hermitian hamiltonians - the flourishing of the PT symmetric 
approach in the last decade or so. 

Remarkably Bender and collaborators proved also another (even more) funda- 
mental result [llOj: in PT symmetric QM it is possible to have unitary evolution. 
This is not at all a trivial result. Indeed the self-adjointness of the hamiltonian 
is 'responsible' not only for the reality of the energy levels, but also - and I dare 
say most importantly - for the unitarity of time evolution, i.e. the conservation 
of probability in time. Indeed in standard QM unitarity is completely due to the 
hermeticity of the hamiltoniai]|^ the evolution operator is U{t) = e~*^* and is 
unitary as long as H = W . I am not going in details in this respect, I only say 
that it has been proved that it can be defined a proper inner product with respect 
to which PT symmetric QM does have a unitary evolution |110^ IllH 1112] . 

Remark. One key point must be stressed: non-hermitian hamiltonians classi- 
cally mean complex forces. This implies that we have to consider complex dy- 
namical variables x and p. This is completely acceptable in the quantum world, 
once you understand that this only means that the position and momentum op- 
erators are not observables in PT symmetric quantum theory. Observables must 
have real spectrum, this has to be true indeed in any theory. 

However, this is a problem in the corresponding classical limit. One has to 
include somehow complex valued trajectories. Here, I must admit, I see a draw- 
back of this approach. Frankly, I am not aware of any reasonable interpretation 



^Actually it is only a sub-class of PT symmetric operators that have this property. This 
sub-class is said to have an un-broken PT symmetry. See I109| for more details, 
am using here hermeticity and self adjointness as synonimous. Strictly speaking, they are 
not, because of domain subtleties. 
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about it, and it seems to me that in the hterature this point is quite often 
underestimated. 

One more thing to say is that, at present, no empirical clear and definitive 
evidence of the existence of systems with complex hamiltonian has been found. 

It is far beyond the scope of this thesis to discuss the details of this approach, 
let alone the many insights that have been reached in these years in many re- 
spects. However, in the following section I present a class of very well studied 
and archetypal complex PT symmetric hamiltonians, which I find particularly 
instructive and also very useful to understand the mechanism that plays a central 
role in the model studied in |3]. 



4.2. PT-symmetric oscillator 

Take the hamiltonian 

H, = p'^ + x\ixY , (4.4) 

this is complex for e 7^ and PT symmetric]^ This is a generalization of the stan- 
dard harmonic oscillator. We shall call this class the PT symmetric oscillator, 
as is sometimes done in the literature. 



I will review here mainly the classical analysis of the hamiltonians (4.4), i.e. 
the study of the trajectories of a point particle subject to the complex potential 
x'^{ixy. Notice that the richness of the model arises precisely from considering 
it in the complex x plane, as it is natural to do since we have complex forces. 
The equation of motion is 

^=2z(2 + 6)(^x)i+% (4.5) 
that can be integrated, e.g. using the energy first integral H = E, to give 

^^ = ±v/eTM^, (4.6) 

where E is the constant value of the energy of the particle. We set for simplicity 
E = 1] obviously this does not spoil the generality of the following considerations. 
Let us analyze separately the cases e = 0, 1, 2. 



^For e < this hamiltonian is said to have a broken PT symmetry, and the reahty of the 
spectrum is not guaranteed anymore |1081I109] . Thus we restrict to positive (or null) value 
of e. 
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Figure 4.1.: Qualitative plot of the trajectories in the complex x plane of the 
system with hamiltonian (4.4) with e = 0. 



e = 

This is the standard harmonic oscillator. It has two turning points x = —1, +1, 
both on the real axis. The trajectory depends on the choice of a;(0). Here we 
can choose an initial condition wherever we want on the complex x plane. The 
solutions are all periodic of period 27r, and are all oscillating. If we start on the 
real axis, then - as we know from the standard oscillator - we remain on the 



real axis, otherwise we oscillate on ellipses on the complex plane (see figure 4.1). 
Notice that the trajectories are P invariant (reflection with respect to the origin) 
T invariant (reflection with respect to the real axis) and PT invariant (reflection 
about the imaginary axis). 



e = 1 

In this case there are three turning points, namely the solution of the equation 

ix^ = 1 , (4.7) 



i.e. X- = e~^*'^/^, x+ = e~*'^/^ and xq = i (see flgure 4.2). In this case we 
have two classes of trajectories: periodic and non-periodic. Actually the non- 
periodic trajectories are a null-measure subset of all the trajectories. Namely, 
trajectories starting in a;(0) = ix with x real and greater than or equal to 1, goes 
up to to inflnity on the imaginary axis. In other words, trajectories starting at 
the Xq turning point or somewhere upper but always on the imaginary axis, are 
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Figure 4.2.: Qualitative plot of the trajectories for the case e = 1. 

not periodic. All the other are periodic and have the same period, by virtue of 
Cauchy's theorem. Notice that the trajectories are manifestly PT-invariant, but 
P and T symmetries taken separately are lost. 



This is the most interesting case. Notice that the hamiltonian (4.4) now reads 

H2=p^-x\ (4.8) 

which is actually a real but unbounded hamiltonian. Being it real we can analyze 
it on real trajectories, and we know that the system is completely unstable: the 
particle will roll down the maximum in x = going away at infinity. The energy 
is doomed to fall to — oo. However, we can also expand our horizon and see 
what happens in the complex plane. There are four (complex) turning points, 
solutions of = — 1 



± 



^±3j7r/4 



x^ = e— , (4.9) 

two above and two under the real axis (see figure 4.3). The trajectories are 
here divided in three families: one encircling the top turning points xf2j one 
around the bottom turning points Xi2 stuck on the real axis - again of 

null measure on the whole set. The latter is of course the one we talked about 
previously. Indeed it is unstable, in the sense that it is not periodic and goes to 



infinity. So, here you are: while on standard basis you look to hamiltonian (4.8) 
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and think 'unstable, good for nothing', in the complex plain you'd better bet on 
the periodicity/stability of the same system described by (4.8). 

This last example shows the key point: a manifest instability can be 'absorbed' 
by the 'complexity'. Or better (but less suggestively), what is manifestly unsta- 
ble in a "real" world, might not be so in a "complex" world. 



With these brief examples I end this introduction to PT symmetric QM. 
Actually, this was more a presentation of the ideas behind PT symmetry rather 
than of the theory itself. Indeed there is a huge amount of work, both in the 
foundations of the theory and in more border-line issues. The interested reader 
is encouraged to look to the very good reviews by Bender |108l 1109] . 

4.3. Introduction and resume of the paper p] 

As I pointed out at the beginning of the chapter, complex (non-hermitian) hamil- 
tonians with PT symmetry have been vigorously investigated in Quantum Me- 
chanics and Quantum Field Theory |106l 1107) . A possibility of applications to 
quantum cosmology has been pointed out in |113j . In |^ we mainly focused 
attention on complex field theory. We explored the use of a particular complex 
scalar field lagrangian, whose solutions of the classical equations of motion pro- 
vide us with real physical observables and well-defined geometric characteristics. 
In the said paper we proposed a cosmological model inspired by PT sym- 
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metric Quantum Theory, choosing potentials so that the equations of motion 
have classical phantom solutions for homogeneous and isotropic universe. Mean- 
while quantum fluctuations have positive energy density and thus ensuring the 
stability around a classical background configuration. 

We considered the complex extension of matter lagrangians requiring the re- 
ality of all the physically measurable quantities and the well-definiteness of ge- 
ometrical characteristics. It is worthwhile to underline here that we considered 
only real space-time manifolds. Attempts to use complex manifolds for studying 
the problem of dark energy in cosmology can be found here |114t I115[ I116| . 

In [2], we start with the model of two scalar fields with positive kinetic terms. 
The potential of the model is additive. One term of the potential is real, while 
the other is complex and PT symmetric. We find a classical complex solution 
of the system of the two Klein-Gordon equations together with the Friedmann 
equation. The solution for one field is real (the "normal" field) while the solution 
for the other field is purely imaginary, realizing classically the phantom behavior. 
Moreover, the effective lagrangian for the linear perturbations has the correct 
potential signs for both the fields, so that the problem of stability does not arise. 
However, the background (homogeneous Friedmann) dynamics is determined by 
an effective action including two real fields one normal and one phantom. As a 
byproduct, we notice that the phantom phase in the cosmological evolution is 
inevitably transient. The number of phantom divide line crossings, (i.e. events 
such that the ratio w between pressure and energy density passes through the 
value —1) can be only even and the Big Rip never occurs. The avoidance of Big 
Rip singularity constitutes an essential difference between our model and well- 
known quintom models, including one normal and one phantom fields [TT7| 1118] . 
The other differences will be discussed in more detail later. 

What follows is so organized: in sections 4A and 4^ the cosmological model 
we considered is described, together with a brief explanation of the "interplay" 



between PT symmetric Quantum Mechanics and two-field models; in section 4.6 



I present the results of the qualitative analysis and of the numerical simulations 
for the dynamical system under consideration; conclusions and perspectives are 



presented in the last section 4.7 



4.4. Phantom and stability 



As I pointed out at the end of section 1.1.3 the great problem of phantom fields 
is the quantum instability [5^ lUU] . Their hamiltonian is unbounded from below, 
thus its quantum fluctuations grow exponentially. This should remind the reader 



of the PT oscillator case |4.2[ which indeed I will again discuss, in a slightly more 
specific context. The idea is to exploit the PT symmetric framework to "go 
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around" (in the complex plane) this kind of instability. 

We shall study the fiat Friedmann cosmological model described by the (usual) 



FLRW metric ( |1.2[ ). 

Let us consider the matter represented by scalar fields with complex poten- 
tials. Namely, we shall try to find a complex potential possessing the solutions 
of classical equations of motion which guarantee the reality of all observables. 
Such an approach is of course inspired by the quantum theory of the PT sym- 
metric non-hermitian hamiltonians, whose spectrum is real and bounded from 
below. Thus, it is natural for us to look for lagrangians which have consistent 
counterparts in the quantum theory. 

Let me elucidate how the phantom-like classical dynamics arises in such la- 
grangians. For this purpose, at risk of being repetitive, choose the one-dimensional 
PT symmetric potential of an-harmonic oscillator V^'^^'^\q) = Xq^{iqy, < e < 



2 which has been rapidly analyzed above (4.2). For illustration, let me take here 
the most interesting case of e = 2,V^^\q) = — Ag^. As said above, the classi- 
cal dynamics for real coordinates q{t) offers the infinite motion with increasing 
speed and energy or, in the quantum mechanical language, indicates the absence 
of bound states and unboundness of energy from below. However, just there is a 
more consistent solution which, at the quantum level, provides the real discrete 
energy spectrum, certainly, bounded from below. It has been proven, first, by 
means of path integral |106| and, further on, by means of the theory of ordi- 
nary differential equations jll9] . In fact, this classically "crazy" potential on a 
curve in the complex coordinate plane generates the same energy spectrum as 
a two-dimensional quantum an-harmonic oscillator V^^\q) = X{qf + 5'|)/4 with 
real coordinates gi 2 in the sector of zero angular momentum |106) . Although 
the superficially unstable an-harmonic oscillator is well defined on the essen- 
tially complex coordinate contour any calculation in the style of perturbation 
theory (among them the semi-classical expansion) proceeds along the contour 
with a fixed complex part (corresponding to a "classical" solution) and varying 
unboundly in real direction. In particular, the classical trajectory for V^'^\q) 
with keeping real the kinetic, (g)^ and potential, — Ag^(t) energies (as required 
by its incorporation into a cosmological scenario) can be chosen imaginary. 



q = i^, e = -2Ae^ e = c-xe, c>o 

which obviously represents a bounded, finite motion with |^| < (C/A)^/^. Such 
a motion supports the quasi-classical treatment of bound states with the help of 
Bohr's quantization. Evidently, the leading, second variation of the lagrangian 
around this solution, q{t) = i^{t) + Sq{t) gives a positive definite energy, 

£(2) = prnit) -H, H = ^/(t) + uemsqit))' 
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realizing the perturbative stability of this an-harmonic oscillator in the vicinity of 
imaginary classical trajectory. It again reflects the existence of positive discrete 
spectrum for this type of an-harmonicity. However the classical kinetic energy 
— is negative, i.e. it is phantom-like. 

In a more general, Quantum Field Theory setting let us consider a non- 
hermitian (complex) lagrangian of a scalar field 

L=h,^d'^<P*-V{<P,^*) , (4.10) 
with the corresponding action, 

S[<P,g] = J d''xV^(^L + ^Rig)^ , (4.11) 

where g stands for the determinant of a metric g^'^ and R{g) is the scalar curva- 
ture term and the Newton gravitational constant is as usual normalized to S/Stt 
to simplify the Friedmann equations. 

We employ potentials V($, $*) satisfying the invariance condition 

(l^($, $*))* = 1/($*,$) , (4.12) 

while the condition 

(\/($, $*))* = ■!/($, $*) , (4.13) 

is not satisfied. This condition represents a generalized requirement of {C)PT sym- 
metry. 

Let's define two real fields, 

0^ + !($_$*). (4.14) 

Then, for example, such a potential can have a form 

^($, $*) = V^($ + $*,$- $*) = V{(j),ix) , (4.15) 

where V{x, y) is a real function of its arguments. In the last equation one can 
recognize the link to the so called {C)PT symmetric potentials if to supply the 
field X with a discrete charge or negative parity. When keeping in mind the 
perturbative stability we impose also the requirement for the second variation 
of the potential to be a positive definite matrix which, in general, leads to its 
PT symmetry |69| [70] . 

Here, the functions and x appear as the real and the imaginary parts of the 
complex scalar field $, however, in what follows, we shall treat them as inde- 
pendent spatially homogeneous variables depending only on the time parameter 
t and, when necessary, admitting the continuation to complex values. 
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4.5. Cosmological solution with classical 
phantom field 

It appears that among known PT symmetric hamiltonians (lagrangians) pos- 
sessing the real spectrum one which is most suitable for our purposes is that 
with the exponential potential. It is connected with the fact that the properties 
of scalar field based cosmological models with exponential potentials are well 
studied[74| [761 11201 I121j . In particular, the corresponding models have some 
exact solutions providing a universe expanding according to some power law 
a{t) = aof. We shall study the model with two scalar fields and the additive po- 
tential. Usually, in cosmology the consideration of models with two scalar fields 
(one normal and one phantom, i.e. with the negative kinetic term) is motivated 
by the desire of describing the phenomenon of the so called phantom divide line 
crossing. At the moment of the phantom divide line crossing the equation of 
state parameter w = p/e crosses the value w = —1 and (equivalently) the Hub- 
ble variable h has an extremum. Usually the models using two fields are called 
"quintom models" |69]-[72]. As a matter of fact the phantom divide line crossing 
phenomenon can be described in the models with one scalar field, provided some 
particular potentials are chosen[67t i68j or in the models with non-minimal cou- 
pling between the scalar field and gravity |122) . However, the use of two fields 
make all the considerations more simple and natural. In the framework we are 
considering here, the necessity of using two scalar fields follows from other re- 
quirements. We would like to implement a scalar field with a complex potential 
to provide the effective phantom behavior of this field on some classical solutions 
of equations of motion. Simultaneously, we would like to have the standard form 
of the effective Hamiltonian for linear perturbations of this field. The combi- 
nation of these two conditions results in the fact the background contribution 
of both the kinetic and potential term in the energy density, coming from this 
field are negative. To provide the positivity of the total energy density which 



is required by the Friedmann equation (1.8a) we need the other normal scalar 



field. Thus, we shall consider the two-field scalar lagrangian with the complex 
potential 

; 2 • 2 

L = y + ^ - Ae"^ + Be'^'' , (4.16) 

where A and B are real, positive constants. This lagrangian is the sum of 
two terms. The term representing the scalar field is a standard one, and 
it can generate a power-law cosmological expansion [7^ [7S| 1120^ I121| . The 
kinetic term of the scalar field x is also standard, but its potential is complex. 
Notice that the exponential potential has been analyzed in great detail in |123l 
1124 ]. The most important feature of this potential is that the spectrum of the 
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corresponding Hamiltonian is real and bounded from below, provided correct 
boundary conditions are assigned. 

Inspired by this fact we are looking for a classical complex solution of the 
system, including two Klein-Gordon equations for the fields and x'- 



4> + 3h^ + Aae'^* = , 

X + 3hx- iB(3e'^^ = , 
and the Friedmann equation 



/,2 ^2 



(4.17) 
(4.18) 

(4.19) 



2 2 

The classical solution which we are looking for should provide the reality and 



positivity of the right-hand side of the Friedmann equation (4.19). The solution 



where the scalar field is real, while the scalar field x is purely imaginary 



X = i^, ^ real , 



(4.20) 



uniquely satisfies this condition. Moreover, the lagrangian (4.16) evaluated on 



this solution is real as well. This is remarkable because on homogeneous solutions 
the lagrangian coincides with the pressure, which indeed should be real. 



Substituting the equation (4.20) into the Friedmann equation (4.19) we shall 
have 



2 2 



(4.21) 



Hence, effectively we have the Friedmann equation with two fields: one (0) is 
a standard scalar field, the other (^) has a phantom behavior, as we pointed 
out above. In the next section we shall study the cosmological dynamics of the 



(effective) system, including (4.21), (4.17) and 



^ + 3h^- B(3e 







(4.22) 



The distinguishing feature of such an approach rather than the direct con- 
struction of phantom lagrangians becomes clear when one calculates the linear 
perturbations around the classical solutions. Indeed the second variation of the 
action for the field x gives the quadratic part of the effective lagrangian of per- 
turbations: 

L,ff = \5x -Bp\^^^-{5x)\ (4.23) 

where xo is a homogeneous purely imaginary solution of the dynamical sys- 
tem under consideration. It is easy to see that on this solution, the effective 



lagrangian (4.23) will be real and its potential term has a sign providing the 
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stability of the background solution with respect to linear perturbations as the 
related Hamiltonian is positive, 

H!,% = l6^ + Bp'e-^^"{6xf, Sn ^ 6x. (4.24) 

Let us list the main differences between our model and quintom models, using 
two fields (normal scalar and phantom) and exponential potentials |117[ 1118] . 
First, we begin with two normal (non-phantom) scalar fields, with normal ki- 
netic terms, but one of these fields is associated to a complex (PT symmetric) 
exponential potential. Second, the (real) coefficient multiplying this exponential 
potential is negative. Third, the background classical solution of the dynamical 
system, including two Klein-Gordon equations and the Friedmann equation, is 
such that the second field is purely imaginary, while all the geometric charac- 
teristics are well-defined. Fourth, the interplay between transition to the purely 
imaginary solution of the equation for the field x the negative sign of the 
corresponding potential provides us with the effective lagrangian for the linear 
perturbations of this field which have correct sign for both the kinetic and po- 
tential terms: in such a way the problem of stability of the our effective phantom 
field is resolved. Fifth, the qualitative analysis of the corresponding differential 
equations, shows that in contrast to the quintom models in our model the Big 
Rip never occurs. The numerical calculations confirm this statement. 

In the next section we shall describe the cosmological solutions for our system 
of equations. 

4.6. Cosmological evolution 

First of all notice that our dynamical system permits the existence of cosmolog- 
ical trajectories which cross the phantom divide line. Indeed, the crossing point 
is such that the time derivative of the Hubble parameter 

h = -li^'-i') (4.25) 

is equal to zero. We always can choose = at t = tpoi provided the 
values of the fields (j){tpDL) and ^{tpDi) are chosen in such a way, that the 
general potential energy Ae°'^ — Be~^^ is non-negative. Obviously, tpDL is the 
moment of the phantom divide line crossing. However, the event of the phantom 
divide line crossing cannot happen only once. Indeed, the fact that the universe 
has crossed phantom divide line means that it was in effectively phantom state 
before or after such an event, i.e. the effective phantom field ^ dominated over 
the normal field 0. However, if this dominance lasts for a long time it implies 
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that not only the kinetic term — 1^^/2 dominates over the kinetic term 0^/2 but 
also the potential term —Bexp{—P^) should dominate over Aexp{a(f)); but it is 
impossible, because contradicts to the Friedmann equation (4.21). Hence, the 



period of the phantom dominance should finish and one shall have another point 
of phantom divide line crossing. Generally speaking, only the regimes with even 
number of phantom divide line crossing events are possible. Numerically, we 
have found only the cosmological trajectories with the double phantom divide 
line crossing. Naturally, the trajectories which do not experience the crossing at 
all also exist and correspond to the permanent domination of the normal scalar 
field. Thus, in this picture, there is no place for the Big Rip singularity as well, 
because such a singularity is connected with the drastically dominant behavior 
of the effective phantom field, which is impossible as was explained above. The 
impossibility of approaching the Big Rip singularity can be argued in a more 
rigorous way as follows. Approaching the Big Rip, one has a growing behavior 
of the scale factor a{t) of the type a{t) ~ {tsR — t)^^-, where g > 0. Then the 
Hubble parameter is 

h{t) = — — (4.26) 



and its time derivative 



hit) 



(tBR - tf 



Then, according to equation (4.25) 



if - 

2^ 2^ 



3(tBi? - ty ■ 



(4.27) 



Substituting equations (4.27), (4.26) into the Friedmann equation (4.21), we 
come to 



q 



itBR-tf 'OtBR-tf 



In order for this to be satisfied and consistent, the potential of the scalar field 
should behave as l/{tBR. — t)"^. Hence the field should be 



bo \n{tBR - t) 

a 



(4.28) 



where 0o is an arbitrary constant. Now substituting equations (4.26) and (4.28) 



into the Klein-Gordon equation for the scalar field (4.17), the condition of the 



cancellation of the most singular terms in this equation which are proportional 

to l/{tBR — ty reads 



2 + 6g + Aa^ exp(a</)o) = . 



(4.29) 
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This condition cannot be satisfied because all the terms in the left-hand side of 
equation (4.29) are positive. This contradiction demonstrates that it is impossi- 
ble to reach the Big Rip. 

Now I would like to describe briefly some examples of cosmologies contained in 
our model, deduced by numerical analysis of the system of equations of motion. 
In figure 4.4 a double crossing of the phantom divide line is present. The evolu- 
tion starts from a Big Bang-type singularity and goes through a transient phase 
of super-accelerated expansion ("phantom era"), which lies between two crossings 
of the phantom divide line. Then the universe undergoes an endless expansion. 
In the right plot I present the time evolution of the total energy density and of 
its partial contributions due to the two fields, given by the equations 

e = + , 

which clarify the roles of the two fields in driving the cosmological evolution. 
The evolution presented in figure 4^ starts with a contraction in the infinitely 
remote past. Then the contraction becomes super-decelerated and turns later in 
a super-accelerated expansion . With the second phantom divide line crossing 
the "phantom era" ends; the decelerated expansion continues till the universe 
begins contracting. After a finite time a Big Crunch-type singularity is encoun- 
tered. From the right plot we can clearly see that the "phantom era" is indeed 
characterized by a bump in the (negative) energy density of the phantom field. 
In figure 4.6 the cosmological evolution again begins with a contraction in the 
infinitely remote past. Then the universe crosses the phantom line: the contrac- 
tion becomes super-decelerated until the universe stops and starts expanding. 
Then the "phantom era" ends and the expansion is endless. 
In figure 4/7 the evolution from a Big Bang-type singularity to an eternal ex- 
pansion is shown. The phantom phase is absent. Indeed the phantom energy 
density is almost zero everywhere. 



4.7. Conclusions 

As was already said many times, the data are compatible with the presence of 
the phantom energy, which, in turn, can be in a most natural way realized by 
the phantom scalar field with a negative kinetic term. However, such a field 
suffers from the instability problem, which makes it vulnerable. Inspired by the 
development of PT symmetric Quantum Theory we introduced the PT symmet- 
ric two-field cosmological model where both the kinetic terms are positive, but 
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h 




Figure 4.4.: (left) Plot of the Hubble parameter representing the cosmological 
evolution. The evolution starts from a Big Bang-type singularity 
and goes through a transient phase of super-accelerated expansion 
("phantom era"), which lies between two crossings of the phantom 
divide line (when the derivative of h crosses zero). Then the universe 
expands infinitely, (right) Plots of the total energy density (blue), 
and of the energy density of the normal field (purple) and of the 
phantom one (green). 




Figure 4.5.: (left) The evolution starts with a contraction in the infinitely remote 
past. At the first phantom divide line crossing the contraction be- 
comes super-decelerated and turns in a super-accelerated expansion 
when h crosses zero. The second crossing ends the "phantom era"; 
the decelerated expansion continues till the universe begins contract- 
ing. In a finite time a Big Crunch-type singularity is reached, (right) 
Plots of the total energy density (blue), and of the energy density 
of the normal field (purple) and of the phantom one (green). 
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Figure 4.6.: (left) The cosmological evolution begins with a contraction in the 
infinitely remote past. With the first phantom divide line crossing 
the contraction becomes super-decelerated until the universe stops 
[h = 0) and starts expanding. With the second crossing the "phan- 
tom era" ends and the expansion continues infinitely, (right) Plots 
of the total energy density (blue), and of the energy density of the 
normal field (purple) and of the phantom one (green). 
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Figure 4.7.: (left) Evolution from a Big Bang-type singularity to an infinite ex- 
pansion, without any crossing of the phantom divide line. This 
evolution is thus guided by the "normal" field 0. (right) Plots of the 
total energy density (blue), and of the energy density of the normal 
field (purple) and of the phantom one (green). Notice that the en- 
ergy density of the phantom field (green) is very close to zero, thus 
the total energy density is mainly due to the standard field. 
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the potential of one of the fields is complex. We studied a classical background 
solution of two Klein-Gordon equations together with the Friedmann equation, 
when one of this fields (normal) is real while the other is purely imaginary. The 
scale factor in this case is real and positive just like the energy density and the 
pressure. The background dynamics of the universe is determined by two effec- 
tive fields - one normal and one phantom - while the lagrangian of the linear 
perturbations has the correct sign of the mass term. Thus, so to speak, the 
quantum normal theory is compatible with the classical phantom dynamics and 
the problem of instability is absent. 

As a byproduct of the structure of the model, the phantom dominance era is 
transient, the number of the phantom divide line crossings is even and the Big 
Rip singularity is avoided. 
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The present chapter is intended to present the results of ^j. This is kind of 
off-topic, with respect to the preceding work, but not at all by far. Indeed the idea 
is to study the reconstruction of a single scalar field model, able to reproduce a 
specific dynamical evolution, which seemed to us particularly interesting, since 
it presents some kind of mild version of the Big Bang singularity. Dynamical 
analysis is performed on the model and its phase space is divided into classes of 
qualitative different cosmic evolutions. Here there is no crossing of the phantom 
divide line, nor phantom fields. 



5.1. Introduction and review of paper [4| 

We have repeatedly said that the discover of cosmic acceleration [6, 7J has stim- 
ulated the research of new cosmological models. This development of model- 
designing art has revealed cosmological evolutions possessing various types of 
singularities, sometimes very different from the traditional Big Bang and Big 
Crunch. The most popular between them is, perhaps, the Big Rip cosmolog- 
ical singularity [57, 58j arising in super-accelerating models driven by some 
kind of phantom matter. Other types of singularities are sudden singulari- 
ties [Ml M Dsn, Big Brake [7S|, and so on |I2S]-|IM]. Here we would like 
to study the singularities which are close to the known Big Bang, Big Crunch 
and Big Rip singularities, but arising at finite values of the cosmological fac- 
tor (different from zero and infinity as well). Similar singularities were recently 
considered in |l32l [1331 EMI. 

In the paper under consideration we construct potentials which can drive the 
cosmological evolution towards (or from) such singularities. Combining qualita- 
tive and numerical methods we study the set of possible cosmological histories 
in the suggested models to show that the presence of such singularities in a 
cosmological model under consideration depends essentially on initial conditions 
and that the same model can accommodate qualitatively different cosmological 
scenarios. 



The structure of what follows is: in section 5.2 we construct some potentials 
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corresponding to evolutions with "mild" singularities. In the third section 5.3 



we 



analyze their dynamics. The conclusion 5.4 is devoted to an interpretation of 
the obtained results. 



5.2. Construction of scalar field potentials 

Let me repeat here very briefly the idea of the reconstruction of potentials for 



cosmological models (cfr. section 2.1). 



We shall consider fiat Friedmann models with the metric (1.2) The Hubble 



parameter h(t) = a/a satisfies the Friedmann equation (1.8a) and also its 'ma- 



nipulation' (1.9) 



If the matter is represented by a spatially homogeneous minimally coupled 
scalar field, then the energy density and the pressure are given by the formulae 



(2.1) and (2.2), respectively; while equations (2.5) and (2.6) give the e xpre ssion 



of the potential and of as functions of time. Integrating equation (2.6) one 
can find the scalar field as a function of time. Inverting this dependence we can 
obtain the time parameter as a function of (j) and substituting the corresponding 



formula into equation (2.5) one arrives to the uniquely reconstructed potential 



V{(j)). It is necessary to stress that this potential reproduces a given cosmological 
evolution only for some special choice of initial conditions on the scalar field and 
its time derivative. 

It is known that the power-law cosmological evolution is given by the Hubble 
parameter h{t) ~ |. We shall look for a "softer" version of the cosmological 
evolution given by the law 

Kt) = ^ , (5.1) 

where 5 is a positive constant and 0<a<l. At t = a singularity is present, 
but it is different from the traditional Big Bang singularity. Indeed, integrating 
we obtain 



In 



a(0) 



1 — a 



If t > the right-hand side of equation (5.2) is finite and hence one cannot 
have a(0) 



in the left-hand side of this equation, because it would imply a 

a(0) 



contradiction, making In divergent. Hence a(0) > 0, while 



«(0) — exp 



S 



-t 



l-Q 



1 — a 



-7- OO 



(5.3) 



This type of singularity can be called "mild" Bing Bang singularity because the 
cosmological radius is finite (and non-zero) while its time derivative, the Hubble 
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variable and the scalar curvature are singular. It is interesting to note that 
when t — )■ oo both a(t) and d(t) tend to infinity, but they do not encounter any 
cosmological singularity because the Hubble variable and its derivatives tend to 
zero. 

Let us reconstruct the potential of the scalar field model, producing the cos- 



mological evolution (5.1) using the technique described above. Equation (2.6) 
gives 



,2 ^ g+l 

±\l-aSt- — 



(5.4) 



We shall choose the positive sign, without loosing generality. Integrating, we get 

m- 



-ab - 
3 1 



a 



(5.5) 



up to an arbitrary constant. Inverting the last relation we find 



2aS 



1/2 



a 



(5.6) 



Hence, using equation (2.5) we obtain 

V{<P) = 



aS 



3 1-a 
2aS 2 * 



4q 

1-Qf 



3 l-g 
2aS 2 



2{a + l) ' 



(5.7) 



This potential provides the cosmological evolution (5.1 ) if initial conditions com- 



patible with equations (5.4) and (5.5) are chosen. Naturally, there are also other 



cosmological evolutions, generated by other initial conditions, which will be stud- 
ied in the next section. 



5.3. The dynamics of the cosmological model 
with = I 

In order to achieve some simplification of calculations we shall consider a par- 
ticular model, namely the one with the choice a = \. In this case 

a{t) = a(0)e2^^ , (5.8) 

and 

^ ' (730/2)4 3(v^0/2)6 ■ ^ ' ' 
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The Klein-Gordon equation reads 



+ 30 sign /i J -02 + —— + = . 

^^^V^ (v^0/2r 3(v^0/2)6 (v^0/2)5 (730/2)7 

(5.10) 

This equation is equivalent to the dynamical system 



X 



c2 I 165^ 3254 

(v^</./2)4 3(V3«!./2)6 

32v^S-* _ 32^35-* 
(v^0/2)5 (V30/2)7 • 



-3x sign(/^)^/^ + - (5.11) 

+ 



The qualitative analysis of dynamical systems in cosmology was presented in 
detail in [85j. 

First of all, let us notice that the system has two critical points: = ±;^, x = 0. 
We consider the linearized system around the point with the positive value of 0: 

(p= X, 

x= -3x sign(/i)^52 + 32-35V, 



....5..,,., + 32 ■35V, ^^"^^^ 

where (f = (j) — 2/ y/S. 

The Lyapunov indices for this system are (for h > 0) 

Ai = -Sv^^^ , (5.13) 

A2 = 4^35^ (5.14) 

For negative h, corresponding to the cosmological contraction, the signs of Ai 
and A2 are changed. The eigenvalues are real and have opposite signs, hence both 
the critical points are saddles. The universe being in one of these two saddle 
points means that it undergoes a de Sitter expansion or contraction, according 
to the sign of h, with the value of h = ho given by 

452 

ho=^. (5.15) 

For each saddle point there are four separatrices which separate four classes of 
trajectories in the phase plane x, corresponding to four types of cosmological 
evolutions. 

In order to simplify the study of the dynamics let us note that the potential 
is an even function of the scalar field and that the saddle points are also 
symmetrical with respect to the x axis. Thus, it is sufficient to consider only 
one of this saddle points. We shall carry out our qualitative analysis taking 
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/// 




Figure 5.1.: Plot of the potential V^(0) as given by equation (5.9). Here only the 
positive 0-axis is shown, which is enough to understand the behav- 
ior thanks to the parity of the function. We have also delineated 
the four different dynamical behaviors of the scalar field, clearer to 
understand taking into account the phase portrait (see figure 5.2). 




()-axis phase portrait for the dynamical system (5.11) with 
The four separatrices of the saddle point 



Figure 5.2.: Positive 

/i > 0. The four separatrices of the saddle point (a,/?, 7, 5) in- 
dividuate four regions {I, II, III, IV) with different behaviors of the 
trajectories, as explained in the text. 
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into account both figure |5.1[ giving the form of the potential, and figure |5.2 



representing the phase portrait in the plane 0, x. 

First let us consider trajectories which begin at the moment t = 0, when 
the initial value of the scalar field is infinite, its potential is equal to zero and 
the time derivative of the scalar field is infinite and negative. In terms of the 
figure 5.1 it means that we consider the motion of the point beginning at the 
far right on the slope of the potential hill and moving towards the left (i.e. 
towards the top of the hill) with an infinite initial velocity. Such a motion for 
h > describes a universe born from the standard Big Bang singularity. Further 
details of this evolution depend on the asymptotic ratio between absolute values 
of <p and at t — 0. If this ratio is smaller than some critical value then the 
scalar field does not reach the top of the hill and at some moment it begins 
to roll down back to the right. During this process of rolling down the scalar 
field increases, the potential is decreasing and the velocity becomes positive 
and increasing. However the universe expansion works as a friction and at some 
moment its influence becomes dominant causing an asymptotic damping to zero 
of the velocity. The universe expands infinitely with h{t) — )■ 0. In the phase 
portrait (figure 5.2) such trajectories populate the region //. This region is 
limited by the separatrices /3 and 7. The first one corresponds to the positive 
(for h > 0) eigenvalue A2, while 7 corresponds to Ai. If the ratio <p/(j) has the 
critical value, then the scalar field reaches asymptotically the top of the hill of 
the potential, meaning that the universe becomes asymptotically de Sitter: in 
the phase portrait it is nothing but the curve 7. 

When the ratio introduced above is larger than the critical one we encounter a 
different regime. In this case the scalar field passes with non-vanishing velocity 
the top of the hill and begins to roll down in the abyss on the left. The velocity 
is growing, but the potential becomes negative and at some moment the total 
energy density of the scalar field vanishes together with the Hubble parameter h: 
this means that the universe starts contracting. This contraction provides the 
growing of the absolute value of the velocity of the scalar field and the kinetic 
term again becomes larger then the potential one. Moreover, both the terms in 
the Klein-Gordon equation increase the absolute value of (p. One can easily show 
that the regime in which the time derivative becomes equal to —00 at some 
finite value of is impossible, because it implies a contradiction between the 
asymptotic behavior of different terms in equation (5.10). Thus, the universe 
tends to the singularity squeezing to the state with the value of equal to zero 
and an infinite time derivative 0. To understand which kind of singularity the 
universe encounters, we need some detail about the behavior of the scalar field. 
Let us suppose that, approaching the singularity at some moment to, the scalar 
field behaves as 

(l){t) = Mto-tr , (5.16) 
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2 



where < /i < 1. Then the first and second time derivative are 

m = -/i0o(to - tr~\ m = - 1)00(^0 - tf^. (5.17) 

The potential behaves as 



To have the Hubble variable well defined we require that the kinetic term is larger 



than the absolute value of the negative potential term (5.18), i.e. 2fi — 2< — 6/i, 
or < |- Now two opposite cases may hold: (i) the friction term in the Klein- 
Gordon equation could dominate the potential term or (ii) the opposite situation. 
For (i) to be the right case, one should require 2/i — 2 < —7^ or /i < |. In this 
situation the asymptotic behavior of the second time derivative of should be 
equal to that of the friction term, or, in other words /i — 2 = 2/i — 2, that is = 0, 
which obviously is not relevant. Thus we have to consider the range | < < |. 
In this case the potential term should be equal to the second time derivative of 
(j), which implies: 

f^=l (5.19) 

and 

f • (5.20) 
Substituting the values of /i and 0o into the expression for h{t), we obtain 

h{t) = -^=1= . (5.21) 

Thus, we see that the singularity we are approaching is of the "mild" Big Crunch 
type. 



In the phase portrait (figure 5.2) these trajectories occupy the region /// 



limited by the separatrices 7 and 6. The cosmological evolutions run from the 



Big Bang singularity to the mild Big Crunch one. However, the figure 5.2 is not 
sufficient to describe the complete behavior of the universe under consideration, 
because at some moment the Hubble variable changes sign and we should turn 



to the figure 5^, giving the phase portrait for the contracting universe h < 0. 
Note that increasing the velocity with which the scalar field overcomes the top 
of the hill, implies increasing the moment of time when the point of maximal 
expansion of the universe is reached. The limiting case in which this moment 
tends to infinity corresponds to the separatrix 6. 

The third regime begins from the mild Big Bang singularity, when the scalar 



field is equal to zero and its time derivative is infinite and positive. In figure 5.1 
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Figure 5.3.: Phase portrait for the dynamical system under discussion (equa- 
tion (5.11)) with h < 0, i.e. describing evolutions of the universe 
characterized by contraction. As for the /i-positive case, we can see 
the four separatrices of the saddle point and the corresponding four 
regions of different dynamical behaviors of the trajectories. 



that situation is represented by the point climbing from the abyss to the top 
of the hill. If the velocity term is not high enough, at some moment the field 
stops climbing and rolls back down. During this fall the Hubble variable changes 
sign and the universe ends its evolution in the mild Big Crunch singularity. The 
corresponding trajectories belong to the region IV of our phase plane, bounded 
by the separatrices 6 and a. The universe has its finite life time between the 
mild Big Bang and the mild Big Crunch singularities. The situation when the 
scalar field arrives exactly to the top of the hill and stops corresponds to the 
separatrix a. 



The fourth set of cosmological trajectories is generated by the scalar field 
climbing from the abyss and overcoming the top of the hill with the subsequent 
infinite expansion: the scalar field is infinitely growing and the Hubble parameter 
tends to zero. These trajectories occupy the region / and our original cosmolog- 
ical evolution (5.1) belongs to this family. Naturally there are other four classes 
of cosmological evolutions, which can be easily obtained by inverting the time 
direction. 
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5.4. Conclusions 

Let me sum up the results. Wishing to describe a cosmological evolution be- 
ginning from the singularity characterized by a finite and non-zero initial (or 
final) cosmological radius and an infinite value of the scalar curvature due to 
the infinite value of the Hubble parameter, we have constructed a scalar field 
potential, providing such an evolution. Then, using the methods of qualitative 
analysis of the differential equations, we have shown that the proposed model 
accommodates four different classes of cosmological evolutions, depending on 
initial conditions. Numerical simulations have confirmed our predictions. 

The main results of this work are: (i) the realization of a concrete cosmolog- 
ical model with a scalar field, where finite cosmological radius singularities are 
present and (ii) the complete description of all the possible evolutions of this 
model depending on the initial conditions. It is important to remark that, given 
the fixed scalar field potential, one has different types of evolutions encountering 
different kinds of singularities. 

Let us note, that there are some studies [HSl [M IHZ], [I2H]-[ini], UM 
11331 1134] devoted to the general analysis of various kinds of new cosmological 
singularities. In this work we were not looking for an exhaustive classification of 
different possible cosmological models possessing some kind of singularities, but 
rather we wanted to study in a complete way a particular cosmological model, 
having some interesting properties. 
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Part II. 

Loop Quantum Gravity and 
Spinfoam models 



Introduction 



The title of the present thesis talks about the "relation between Geometry and 
Matter". Just to be clear from the very beginning: I will not talk about the 
coupling of matter with gravity in the Loop Quantum Gravity framework. Even 
though some work has been done in this respeclj^ , this is not quite the sense we 
had in mind when we thought about that title. Let me spend some words about 
it. 

Gravity, or the gravitational field, is indeed a rather peculiar object in physics 
and - I dare say - in Nature. This is true even at the classical level. This 
has become apparent when Einstein discovered the particular "position" of the 
gravitational field with respect to all the other (matter) fields. Indeed the gravi- 
tational field is a field like all the others, with its dynamical equations and so on. 
But it is the background on which all the other fields live as well. All the mat- 
ter fields are defined using the metric g in their equations, they live and "play" 
on a certain metric manifold - the "game rules" - which "is" the gravitational 
field itself. Moreover, there is interplay between them, since Einstein's equations 
tell us that the game rules depend on how the matter plays and viceversa, in a 
recursive/non-linear interplay. This duality of being both a dynamical field and 
the background spacetime on which the fields live is at the heart of the peculiar 
nature of Gravity. 

This is a classical picture. When we quantize the matter fields in the Quantum 
Field Theory framework, we somewhat forget this relation, and assume a fixed 
flat background, without gravity. We can try at best to do Quantum Field 
Theory on curved background, with all its subtleties, still relying on that critical 
and peculiar relation I was talking about; and, most importantly, we still need 
a background spacetime metric on top of which our matter fields live. 

But now a question urges: what if we quantize the gravitational field? We 
can argue that this quantization is doomed to revolutionize once again the way 
in which we think about matter fields living on a spacetime, just as General 
Relativity revolutionized the concept of fixed spacetime. Indeed Loop Quantum 
Gravity gives a picture of space, and of spacetime, which is purely combinatorial 
and relational, where the very concept of "event" looses its meaning: totally 

^See e.g. |135l 11361 1137] for matter coupling in the canonical LQG approach, and |139| - [H5] 
for the spinfoam approach. 
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Introduction 

different from tlie smootfi metric manifold to wliicli we are accustomed. Tlie 
loss of a background is crucial in the way we think matter fields. 

In this precise sense I think that the study of Quantum Gravity is fundamental 
to understand the "relation between Geometry and Matter". 

The following is a brief review of the main aspects of Loop Quantum Gravity 
and of spinfoam theory, the latter being an attempt to understand the quantum 
dynamics of the gravitational field. After this review I will present the research 
work I have done together with Eugenio Bianchi and Carlo Rovelli, during my 
period of research in the Marseille Quantum Gravity group. It is about a rather 
technical aspect of the spinfoam approach to quantum gravitational dynamics. 

Non-perturbative Quantum Gravity: some good 
reasons to consider this way 

Loop Quantum Gravity is a form of canonical quantization of General Relativity. 
We'd better say an attempt in this direction, since the target is still far from being 
reached. However, there are quite amazing and interesting results that should not 
be under-estimated. I will briefiy review the standard approach to the canonical 
quantization of gravity and then pass to a review of the fundamental concepts 
of LQG. But first let us stress the two main ingredients that has permitted LQG 
to become what it is: 

• background independence, 

• focus on connection (rather than metric). 

A few words on each of these concepts are due, before diving into the theory. 
Quantization of gravity means quantization of the gravitational field, the metric 
tensor field g{x) = g^iy{x)dx^ <S> dx'-' , i.e. the geometry of spacetime. One should 
naively expect that such a quantization would bring to some form of "quan- 
tized geometry", meaning a quantization of distances, volumes and such. This 
should happen at a very small distance scale, of course. Dimensionally speaking, 
we know that this scale should be the Planck length, which is of the order of 
10~^^ cm. Moreover, Quantum Mechanics has taught us that linear superposi- 
tions of (quantum) states must be taken into account, in order to give a proper 
(correct) picture of nature. These two (quite speculative) reasonings are here 
used just to introduce the concept of background independence. Background 
independence is a technical way of saying nothing more than this: if one seeks 
for a complete quantization of the metric, one should quantize the metric itself 
and not its small perturbations around a fixed background geometry. In fact, the 
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latter is what is done in the perturbative (or background-dependent) approach 
to Quantum Gravity: to assume the following splitting 

9,u{x) = g%{x) + V(^) , (5.22) 

where (7°^, is treated as a classical field, a background metric, and /i^j, are its 
quantum fluctuations, i.e. it is what one tries to quantize. This is obviously a 
legitimate procedure if one wants to analyze the behavior of a (self-interacting) 
spin-2 particle on a (generally curved) g^^^ background. The problem is the fol- 
lowing: this splitting breaks down at high energies. Namely this kind of theory 
has been proved to be UV-unrenormalizable [144J. Background-dependence sup- 
porters believe that this unrenormalizability is a hint of a more fundamental 
theory still to be found. This is a chance, indeed. Background-independence 
supporters instead believe that this is just a hint of a non-proper way to face the 



problem, namely to believe that the splitting (5.22) should be valid at all energy 
scales. Indeed in the regime we are interested in (namely the Planck scale) it 
is reasonable to argue that the fluctuations won't be small; thus it would be 
a nonsense (read 'wrong') to assume that the geometry (along with the causal 
structure) is determined by g^^ alone. We do not know what is the fabric of 
spacetime at the Planck scale, but - as argued above - it is likely to be some 
quantum superposition, some grain-like structure, no more capable (in general) 
to be described by a smooth fleld + small fluctuations. In this respect back- 
ground independence is both more conservative - because after all makes use of 
the "old" covariant approach - and more radical - since it starts from the belief 
that 'pure' quantum spacetime should be deeply different from that of QFT. 

Let us come to the second issue: the focus on connection. In GR one usually 
takes the metric as the "main character" of the play. It is seen as the fundamental 
fleld of the theory. In LQG instead, the metric tensor is more like a secondary 
player. Actually it is no news at all that GR can be recast in the language 
of differential forms and tetrad fleld (we will soon see in what sense), and the 
language of tetrads - as opposed to that of metric tensor - leads naturally to 
put the focus on connections. At a classical level, everything is equivalent]^ but 
for what concerns quantization, the choice of variables will prove to be crucial, 
and this is one of the key points of the LQG approach. 



^Actually there is one (to my knowledge) difference: the tetrad formulation is "an extension" 



of standard GR. We shall see in section 6.2 in what sense 
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6. Canonical Quantum Gravity: 
from ADM formalism to 
Ashtekar variables 

The present chapter is dedicated to review the basics of canonical quantization of 
gravity. Everything is well-known and established, so I will skip tedious calcula- 
tions and try to present only the key passages of this subject, especially the ones 
which turn to be important in the LQG approad^ 

6.1. Hamiltonian formulation of General 
Relativity 

The starting point is the well known Einstein-Hilbert action: 

S[g,u] = Y^I d'x^R . (6.1) 

In the following I shall always set 167rG = 1 for simplicity, and I will recover the 
constants when to stress the importance of physical scales. 

In order to perform the canonical (i.e. hamiltonian) analysis, one needs to 
know which are the conjugated variables to the metric and then perform the 
Legendre transform. Indeed the hamiltonian formulation looses the manifest co- 
variance, separating time from space - so to say - in order to find the conjugate 
momenta. To this aim, one assumes that it is possible to foliate the spacetime 
manifold (A^ from now on) into spacelike three dimensional hypersurfaces. This 
amounts to say that Ai is diffeomorphic to E x M. This assumptions is actually 
not that restrictive: a theorem by Geroch |145[ |146] proves that it has to be so if 
Ai is globally hyperbolic, namely if it has no causally disconnected regions. Ob- 
viously we are far from saying that we are defining an absolute time coordinate! 
One can choose each timelike direction as "time". This amounts to say that there 
will be an infinity of diffeomorphisms : — t- S x M, each inducing a time 

refer the reader to the good textbook by Baez and Muniain |147| for a thorough and 
somewhat LQG-oriented presentation of this subject. 
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Figure 6.1.: Splitting of dr into normal and tangential components with respect 
to S. 

coordinate r = (j)^t on Ai (namely the pullback of the t coordinate in M). This 
splitting procedure is known as the Arnowitt-Deser-Misner decomposition |148] . 
Now we decompose the timelike (coordinate) vector field dr into its compo- 



nents tangential and normal to S (figure 6.1): 

dr = Nn + N ; (6.2) 

the shift vector (field) belongs to the tangent bundle TS of S; n is the unit 
normal to S (i.e. g{n,v) = Vi; G TS , g{n,n) = —1) and its coefficient is 
called the lapse function. Using the following identified 

^700 = 9{dr, dr) = -N^ + gabN'^N' ; goaN'' = g{dr, N) = NaN'' ^ goa = Na , 

(6.3) 

the expression of the spacetime line element in terms of these variables is: 

ds^ = g^^dx'^dx'' = -{N^- NaN'')dt'^ + 2A'„dMx" + gatdx^dx^ . (6.4) 

Notice that the induced gr-metric on S is given by 

QtMiy = gfiu - n^n^ , (6.5) 

(which has, correctly, goo = 0) and thus we will talk of qab as the metric on S. 
The object = 5'^ — n^n^ is indeed the projector on S 



0; V normal to S , 
v; V tangent to S . 

We want to use these variables qab,Na,N, to rewrite the Einstein-Hilbert ac- 



tion (6.1). Namely, defining the extrinsic curvature 



Kab = Cnqab = ^N-^qab V^A^, VfeA^J , (6.6) 



^From now on latin indices from the beginning of the alphabet (a, b,c. . .) will denote spacelike 
indices ranging from 1 to 3. 
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where Cn stands for the Lie derivative along the n direction (thus the extrinsic 
curvature contains informations on how S is embedded in Ai) and is the 
covariant derivative compatible with qab, it is a matter of algebra to find 



L = = q^/^NR = q^/^N{^R + i:i{K^) - [TiKf] 

The momenta conjugate to the 3-metric are 

dL 



P 



ah 



.ab\ 



(6.7) 



(6.8) 



and there are no momenta for and N°- (they are identically zero). We are now 



able to perform the Legendre transform and write down the action (6.1) as 



S[Tr''\qab,Na,N] 



xE 



-1/2 ab 



+N {q'/'[R^'^ ~ q-'n^bn^' + y-'n'])] , (6.9) 



where ^R is the scalar curvature of qab and vr = 7i"-^qab- Notice that the absence 
of momenta for the lapse function and the shift vector means they carry no 
dynamics, namely that their equation of motion are the following constraint 
equations: 

V'iqab^Tr"') = -2Vf (g-i/V-^) = (6.10) 

and 

Siqab, vr"^) = -q'^'[R^'^ - q-'i^abT^"' + \q^'A = , (6.11) 
usually referred to as the vector constraint and the scalar constraint., respectively. 



With this symbol convention the action (6.9) becomes simply 



S[K''\qab,Na,N]= [ dtd'x [TT'^'qab-NbV'-NS] . (6.12) 

We can also readily identify the hamiltonian density 

nin^', qab, Na, N) = NbV\n^\ qab) + iV5(7r"^ qab) . (6.13) 



Notice that the hamiltonian (6.13) is a combination of constraints, i.e. it van- 



ishes identically on solutions, which is a standard feature of generally covariant 
systems. 

The symplectic structure is very easily found to be 



{p''\x),qM} = Wix-y) 



(6.14) 



and zero otherwise. 
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Imagine we want to (canonically) quantize this system. Firstly we should 
find a Hilbert space in which to represent the algebra (6.14). Notice that our 
configuration space is the space of all the 3- metric on the space manifold S. We 
shall call this space Met(S). Thus what we are looking for is something like 
L^(Met(E), yu) with respect to a suitable measure /i defined on Met(E). This is 
the first great problem of our program of quantization. It is extremely difficult 
to define such a measure! Indeed this space in not even a vector space, since 
Einstein equations are non-linear. 

Moreover, there are other difficulties. Namely We have to impose the con- 
straints (6.10), (6.11). They're expression in terms of configuration variables 
and momenta is extremely unfit to quantization, they are non-polynomial and 
contain nasty factors like y/q creating creates all sorts of ambiguities. 

By the way, these kind of difficulties are much the same as the ones of Yang- 
Mills QFT. In that case, however, the perturbative approach (namely to 'lin- 
earize' the equations of motion) proves to be successful, i.e. the Yang-Mills 
theory is renormalizable. 



6.2. The tetrad/triad formalism 

We have some more work to do in order to have a suitable form to quantize. In 
this section we shall talk a bit about the introduction and the importance of the 
tetrad/triad formalism. Let's for a moment forget the ADM spacetime splitting 
and work covariantly in four dimensions. The tetrad field, or vierbein field, or 
frame field is a rule that assigns to each point of spacetime an orthonormal local 
inertial frame. Let's do it carefully, since this is often a misleading concept. 
You recall that on each spacetime point it is possible to choose a coordinate 
system in which the Levi-Civita connection is zero (but not its derivatives!). 
This set of coordinates is usually called Riemann normal coordinate system, 
or "free fall" coordinates. Indeed they are just the mathematical expression 
of the equivalence principle: this system embodies the possibilities of (locally) 
"eliminating" the gravitational field by "falling" in it. Let us call , J = 0, 1, 2, 3 
these coordinates. Obviously we can't define in general a single fiat coordinate 
system on the whole manifolcQ if so the manifold would be fiat. Their key 
feature is that they are orthonormal, that is 

9ij = ViJ, (6-15) 

locally, i.e. near a specific point: this expression is not true on all the manifold. 
If we move to a "distant" event, it breaks down, in general. However we can 

■^Or, better, these coordinates will be flat only in a small region of spacetime. 
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always change coordinates to another C,^ system, in which this flatness is true in 
that point. The transformation matrix between these free fall coordinates- field 
and a generic coordinate system, is what we call a tetracQ namely 



= ^(^) • (6-16) 



The important thing is that on each point we can put a tetrad ej(a;) which is 
thus a tetrad field. There is a fundamental relation between the tetrad and the 
metric: 

ds^ = rz/jd^^d^-^ = r//jej(x)e;^(x)dx^dx^ = ^^^.dx^dx^ , (6.17) 
which means 

g^,{x) = el{x)et{x)vij • (6.18) 

Thus one can always "reconstruct" the metric from the tetrad. By the way, 
the tetrad carries some redundancies: indeed the metric has 10 independent 
components, while the tetrad 16. This is due to the fact that we can always 
Lorentz-rotate a free fall system to get another equivalent free fall system: 

ei{x)^el{x)=A'j{x)e'^{x) . (6.19) 

This is true whatever Lorentz transformation we attach at each particular point, 
i.e. it is a Lorentz gauge transformation. Indeed the Lorentz group S0(3,l) is a 
6-parameter group, which matches our counting of independent components. 

This is the intuitive physical picture. Actually not so many words are necessary 
to formally introduce the tetrads. Tetrads are simply a basis of orthonormal 
vector fields. I.e. they are four vector fields Cj such that g{ei, Cj) = r}ij. Every 
vector can be expanded in a coordinate basis thus 

e/ = e^jd^ , (6.20) 

while, in general, the tetrad are not a coordinate basis. 

There is however a high brow (but useful) way of saying these same things, in 
the language of bundles. We are indeed picking a local isomorphism (sometimes 
called a trivialization) 

e : M X M" ^ TM , (6.21) 

which permits us to describe (again, locally) the tangent bundle TAi as the 
trivial bundle x (obviously for us n = 4, but this is a rather general 
construction)!^ 



^Actually these are the components of a co-tetrad, since it is a 1-form; it is however common 

to call it tetrad as well. 
^Notice that the existence of this local isomorphism is part of the definition of a difi^erentiable 

manifold, namely to locally look like M". 
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This trivialization is what high brow people call a frame field. Notice that, 
seen from this perspective, everything sounds quite natural: a section of the 
trivial bundle is just a function on Ai with value in M", and we can take a basis 
of these sections Each section of the trivial bundle can thus be expanded as 
s = s^^j If we take e(^/) we obtain a basis of sections of the tangent bundle, 
which can thus be expanded in a coordinate basis 

e(6) = e>id, . (6.22) 

The trivial bundle carries a metric structure, the Minkowski metric. Thus we 
can rise and lower the fiat indices with f]jj. This is a quite powerful tool, but it 
is useful only if we require 

g{v,w) = 'n{e{v),e{w)) , (6.23) 

namely that the local isomorphism sends rj to g. The condition for this to be 
fulfilled is simply 

^(e(0),e(O))=W (6.24) 

that is, we want orthonormal frame fields in order to work with the fiat metric 
in the trivial bundle and then simply "translate" the results in TM. via e. 



To get the expression of equation (6.18) between the metric g and the trivial 



metric, we need the inverse frame field e ^ (the co-tetrad, which is usually called 
tetrad as well, as argued above). Indeed 

g^, = r]{e-\d^), e-\d,)) = r){e[ii, eiO) = vijeie^ . (6.25) 

Remark. An important remark is to be done now. We have a natural step 
forward to do here: recognize that, along with the trivial vector bundle, we get - 



by virtue of the invariance under (6.19) - a principal S0{3, l)-bundle, sometimes 
called the frame bundle. We can think of it as if we were attaching a copy of 
5*0(3, 1) to each point in spacetime; a section of this bundle (namely a choice of 
a specific tetrad among the equivalent ones) is just a gauge choice. 

Now I want to briefiy show (or hint) that the whole General Relativity, and in 



particular the Einstein-Hilbert action (6.1), can be seen as a theory of a tetrad 



field e^(x) and connection. First of all, we have seen that there are two bundles 
in the play: a vector bundle TAi and a principal bundle with S0{3, l)-fiber. 
We have a well defined connection in TAi, the Levi-Civita connection F, defined 



^Notice that here we are in the trivial bundle, which is not the tangent bundle, they are 
only locally isomorphic (this is the same as saying that flat coordinates can be defined only 
locally). This justifies our use of different kind of indices, capital latin letters, which are 
sometimes called flat indices. 
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as the only torsionless connection compatible with the metric g: Wg = 0, with 
V = 9 + r the covariant derivative associated with the parallel transport defined 
by r. We can define a connection on the frame bundle as well, we shall call it 
uj. It will define a covariant differentiation which will "act on" the flat indices, 
i.e. it deflnes a parallel transport on the frame bundle. When we have objects 
with 'a leg in TAi and one in the frame bundle' we have the following "total" 
covariant derivative. 

V,vi = d,vi - T^y^ + uj^jvi . (6.26) 

Now, the requirement is that this covariant derivative has to be compatible with 
the tetrad (in order to have scalar products invariant under parallel transport 
both in the tangent and in the frame bundle), i.e. 

= . (6.27) 

This requirement links together the Levi-Civita connection and the frame bundle 
connection, namely 

^Ij = eiW.e^j , (6.28) 

thus implying that the frame bundle connection u is not arbitrary and contains 
informations on the Levi-Civita connection as well. Such a connection is called 
the spin connection, and we shall indicate it as w(e) when we want to stress 



that it is tetrad-compatible (namely that equation (6.28) holds). One can intu- 
itively think that the Levi-Civita connection is the gauge potential due to the 
diffeomorphism gauge freedom of general relativity, while the a; connection is 
the gauge potential of 50(3, 1) gauge freedom. 
The curvature is deflned as usual 

= du" + w;^ A w-^^ , (6.29) 



and it's of course a 2-form. If we use the spin connection "deflnition" (6.28) we 
have easily 

Fliiuie)) = eiejR'^l, , (6.30) 
where R is the Riemann tensor constructed out of the metric deflned by e 



(through the usual relation (6.18)). This is a central result for our goal: the 
Riemann tensor is nothing but the curvature of the spin connection. With this 
machinery at hand, it is now a matter of calculations to show that the following 
action 

5[e, u] = \eijKL I e' Ae' A F^^(a;) (6.31) 



has the same equations of motions of the Einstein-Hilbert action (6.1). Notice 



that we have explicitly underlined that the action (6.31) is to be considered 
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with both e and u as independent variables. Many of the readers have certainly 
recognized the first-order or Palatini formulation of GR. The variation with 
respect to the connection simply implies that cu is the spin connection on-shell. 
The vanishing variation with respect to the tetrad gives Einstein equations. We 
could have done the same thing with a second order (standard) formulation, in 



which the action is the same as (6.31 ), but we take the curvature to be the one of 



the spin-connection ab initio: again, the equation of motions for e are Einstein 
field equations. However, it was our goal to stress that Einstein gravity can be 
seen as a theory of 'tetrads and connections as independent variables'. 



Remark. Notice that the action (6.31) is actually an extension of the Einstein- 



Hilbert gravity, in fact (6.31) is well defined also in the degenerate case where 



the determinant of the metric is 0, i.e. when the metric is not invertible. In this 
case the tetrad/triad formulation does not crash and still has something to say. 

Let us now go back to our ADM splitting. We introduce the frame field 
formalism, but we shall need it only on the spatial manifold S. Thus we call it 



a triad field. The analogous of equation (6.18) shall be: 

Qab = ^a^b^ij i 



(6.32) 



where now we use i, j. . .to label the internal 3-space. We introduce also the the 
densitized triads 

1 



^- 2 



and 



One can check that 



K. 



C ^ijk^b^c 1 



(6.33) 
(6.34) 

(6.35) 



Thus we now have, for the action (6.12) 



S[E^,Ki,Na,N,N^] 



where 



dt d X 



xE 



E^K' 

I a 



-N,V\El Ki) - NS{E], Ki) - N^G,{E], Ki)] . 
G,{E^,Ki) = e.,,E^^K'^ . 



(6.36) 
(6.37) 



Some remarks: 
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there are here three constraints the vector constraint imposing spa- 
tial diffeomorphism-invariance, the scalar constraint S imposing time diff- 
invariance and the Gauss constraint Gi, which is due to the redundancy in 
using tetrads instead of metrics: namely we have a S0{3) gauge freedom 
(all orthonormal frames in an euclidean 3d manifold are such up to rota- 
tions), which is the spatial analog of the Lorentz gauge freedom analyzed 
above. The Gi constraint is there in order to impose physical states to be 
S'0(3)-invariant. 

This action is written in terms of densitized triads E and their conjugate 



momenta, the K^s defined in (6.34). The latter are, essentially, a redefini- 
tion fo the extrinsic curvature. Thus here the connection is still in a second 
plane: we shall need another change of variables in order to have it clearly 
in the play. 



6.3. The Ashtekar-Barbero variables 

First of all, let us note that the spin connection u^J (based on S) is a 1-form 
with value in the Lie algebra so(3)[^ Thus, we can expand a; on a basis of the 
algebra, for example on the standard Pauli matrices Tj = at/ 2 

= ^ainr , (6.38) 

which (trivially) allows us to speak of w^, with a single internal index i. 

There is now one key change of variables that can be done now. We define a 
new object as follows 

= u;: + , (6.39) 

where 7 is a real arbitrary parameter, called the Immirzi parameter. It is easy 
to see that A is still a connection. Indeed transforms as a vector under local 
5*0(3) transformations, while u transforms as a connection, namely 

4 g{x)co\x)g-\x) + gdg-\x) , 4 g{x)K'\x)g-\x) , (6.40) 

thus A^ transforms as a connection as well. This is known as the Ashtekar- 
Barbero connectior]^ The remarkable fact about this variable, is that its conju- 
gate momentum is just the densitized triad, i.e. 

{Alix), E'^iy)} = ^5^/jix - y) , {Et{x), E]{y)} = {A^^{x), Ai{y)} = , 

(6.41) 

^This is always the case: a connection on a principal G-bundle is always a 1-form with values 
in g- 

^For more details on the Ashtekare-Barbero formulation of GR you can see |150| . For the 
original Ashtekar's contribution see |149| . 
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The Einstein-Hilbert action becomes 

S[ElA^^,Na,N,N']= I dtd'x 

-N,V\EJ, Ai) - NS{EJ, Ai) - N'G,{E^, AD] . (6.42) 

with constraints 

V,iE], Ai) = EJFat - (1 + l')KG, , (6.43a) 

S{E^, Ai) = [e\F^, - 2(1 + i')Km) , (6.43b) 

G,{E], Ai) = DaEt = daEt + e^.uKE"^ , (6.43c) 

where Fab is the curvature of the Ashtekar-Barbero connection and Da is its 
covariant derivative. 

We have reached our goah a formulation of gravity in terms of connection 
and its conjugate momentum, the triad. It is nothing but a 5*0(3) gauge the- 
ory, plus constraints that implement space diffeomorphism invariance and time 
diffeomorphism invariance. In a compact sentence we can say we have reduced 
Einstein gravity to a background independent 5*0(3) Yang-Mills theory. 



6.4. Smearing of the algebra 

In order to give a geometric interpretation to the actual variables (i?J, A\) we 
shall now smear them appropriately. Let us start from the densitized triad Ej{x). 
This object is a 2-form, thus we shall smear it on a surface as follows 

E,{S) = I (x(a))n„dV , (6.44) 
Js 

where ria = ^abc^^ is the normal to the 5* surface. The object here defined, 
Ei{S) is clearly the flux of the triad field across the surface, and is generally 
referred to simply as 'the flux'. 



The Ashtekar-Barbero connection (6.39) is a one form, thus it is natural to 



smear it along a one-dimensional path. Take a path 7 with parametrization 
x"(s) : [0,1] S and take Aa = A^Ti e SU{2), with generators of 5f/(2). 
The we can take the integral 

A = / Alixis))T,^^ds . (6.45) 
7 Jo ds 
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As one can see the triad field is identified on surfaces and its conjugate momen- 
tum, the connection, on paths. This is rather natural on a three dimensional 
manifold E. 

It is useful for future developments, to define the holonomy of the connection 
along the path 7 

h^[A] =Vexp(^j , (6.46) 
where the V symbol stands for path ordering. 



For future reference it is worth to introduce the so-called Hoist action for GR. 



We have seen that Einstein-Hilbert action (6.1) is equivalent to (6.31), which is 



written in terms of tetrad fields and connections, in a first order formalism. We 



can add to (6.31 ) a purely topological sector, which classically does not give any 



contribution to the dynamics 



S[e,co] = ^euKL j A A F^^(a;) + ^ j a Aej A 



(6.47) 



The presence of the Immirzi parameter 7 is of course arbitrary at this level. I 



shall not go into details, however it can be proved that this action (6.47) is a 4d 



manner to introduce the Ashtekar connection. Namely, introducing that topolog- 



ical sector amounts to shift the connection variable with the prescription (6.39). 
See e.g. |T5nl [I5T1 IT53] . 
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7.1. The program 

The program of Loop Quantum Gravity is organized as follows: 

(i) One starts with the canonical quantization of the fundamental variables 
and choosing a representation reproducing the commutator algebra. The 
program of Loop Quantum Gravity starts with choosing the connection 
as the configuration space (position-like) variable. Next the kinematical 
Hilbert space "Hkm will be defined as the space of square integrable func- 
tions on the space of all the possible connections ipiA), with respect to 
some appropriate measure defined on the space, known as the Ashtekar- 
Lewandowski measure. Notice that this Hilbert space is called kinematical 
since it is not the actual Hilbert space of the theory, because we still have 
to deal with the constraints, which - quite obviously - restrict this space 
to one of its subspaces, the physical Hilbert space ?^phys- 

(ii) the following step is to deal with the constraints. The Gauss and Diffeo- 
morphism constraints have a natural unitary action on "Hkm, thus their 
quantization is straightforward. Quite loosely the subspace thus obtained 
is often indicated as H^n as well. Sometimes, to stress the actual distinc- 
tion, one can find H^i^ C H^i^ C "Hkin where the sup-scripts stand for 
Diff constrained and Gauss constrained. Thus one obtains the space of 
solutions to six of the nine constraint equations. 

(iii) the problem and still an open issue of Loop Quantum Gravity, is of course 
to find the space of solution of the scalar constraint, which actually drive 
the time evolution of the system, i.e. the dynamics. Thus the physical 
Hilbert space "Hphys is still to be found, and also a proper scalar product 
between physical states, i.e. quantum transition amplitudes, which are the 
real clue of a every quantum theory. 

We will now briefly sketch each of these three steps, i.e. the definition of the 
kinematical Hilbert space with the Ashtekar-Lewandowski measure, the imple- 
mentation of the Gauss and diffeomorphism constraints, and we shall discuss a 
bit the problems inherent to the quantization of the scalar constraint. 
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I encourage the interested reader to have a look to \151\ \15'2\ HSBJ for very 
good reviews of the basics of LQG and to |154l[l55j for more modern approaches 
to this subject. 



7.2. The kinematical Hilbert space Tiwm 

Firstly we define the space of cylindrical functions based on a graph 7 denoted 

cyi-Q 

A graph 7 is a collection of simple paths which meet at most at their endpoints. 
We call links I all these paths and denote with L their total number. We say 
that7 is in S when all the links (as paths) belongs to S themselves. Given a 
graph 7 C S and a smooth function / : SU{2)^ — )■ C, we define a cylindrical 
function tp^j G Cyl^ as 

i;,j[A] = f{h,[A],hi,[A],...hi^[A]) . (7.1) 



You have to think these paths as the ones of section 6.4, i.e. as the natural 
one-dimensional objects embedded in the spatial hypersurface S upon which the 
connection A is integrated. 

Next we define the space of cylindrical function on S as 

Cyl = U Cyl, , (7.2) 

where the union means over all the graphs in S. Now for the measure. Let's 
define the Ashtekar-Lewandowski measure [157]: 

/^al(V^7,/) = / n^^' fihi„hi^,...hij , (7.3) 

^ IC-y 

where dh is the normalized Haar measure over SU{2). Recall that the Haar 
measure is normalized to 1, thus ^al{^) = 1- 

With this structure one can define the inner product in the space of Cyl: 

(^7,/ I ^i',g) = f^AL{i^*ji^Y,g) = Yl dhi f*{hii ... V) dihi ---hij . (7.4) 

/C7U7' 

Now we are able to define the kinematical Hilbert space Hkm as the Cauchy- 
completion of the space of cylindrical function in the Ashtekar-Lewandowski 



strongly recommend the reading of |156| for a very clear and suggestive presentation of 
this subject. 
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measure. That this is indeed a Hilbert space is not obvious, and it was proved 
by Ashtekar and Lewandowski [157j . 

Having defined properly the Hilbert space, we can now move to find a basis. 
We shall use the Peter- Weyl theorem, which can be seen as a form of generalized 
Fourier expansion. Namely, every function / G L'^[SU{2), d/iHaar] can be written 
as follows 

fig) = J2 V^UT'Dl^g) , (7.5) 

i 

where the D's are nothing but the Wigner 5'?7(2)-representation matrices, j is 
the spin label and run from to infinity in half integer steps, and the / (Fourier) 
coefficients are given by 

' = v/2jTT / diin^rf{g)Dl^,{g) . (7.6) 

JSU{2) 

The completeness relation reads 

5{gh-') = ^(2j + l)TT{D={gh-')) . (7.7) 
j 

We can easily apply this to any cylindrical function 

i^.AA] = n v^j^ E /r.7r'"-"'^^.n,(/^u^]) • • ■ Di^^nAKm 

«=1 jl-jL 

(7.8) 

where the Fourier coefficient are obtained by taking the inner product (7.4) of the 
ijj function with the tensor product of the Wigner matrices. Thus we have found 
a complete orthonormal basis of "Hkinj namely - calling (pl^^, = a/2j~+TD^^, 
the normalized Wigner matrix - the functions 

L 

i=l 

with the remark of taking all the possible graphs in S and all the values of spins 
associated to the links of the graph. 



7.3. The Gauss constraint. H^j^ 

Now we want to restrict the Hilbert space on its gauge invariant subspace, i.e. 
the subspace in which the solutions of the Gauss constraint live. What is to do, 
is simply to take only the states of Hkin which are SU{2) invariant. For this 
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purpose, let us try to understand how SU (2) gauge transformations act on the 



basis (7.9). 



It is indeed very easy to infer the resuh of a (finite) SU{2) gauge transforma- 
tion on a general cylindrical function from the behavior of the holonomy: 

hi[A] ^ gixmhi[A]g-\xil)) , (7.10) 

where x{s) is the parametrized path. Usually, given a path/link /, one calls s{l) 
the source of the path, i.e. the x(0) point, and t{l) the target x{l), thus writing 
the previous equation as: hi[A] A gs(i)hi[A]g~^jy Then, the action on a general 
graph is 

^7, /[^] = f (hi,,... hi J 4 f{gs(i,)hi,[A]g;^ly . . . gs(i^)hi^[A]g-lj , (7.11) 

and we would like to find a basis of the space of cylindrical functions with such 
a property. Let us denote with U{g) the operator that acts the transformation 



(7.11), i.e. - writing it for the basis (7.9) - 



U{9)<PLAhi) = <PLm'i9sii)hig;^)) , (7.12) 

or, more generally, 

L L 
i=l i=l 

Notice that the action of a gauge transformation is on the nodes only: We can 
focus on the request of invariance on a single node, and then extend trivially the 
results to every node of a graph. In order to make things as clear as possible 
and also useful for future reference, let us take a 4-valent node uq. Let us write 
down in a smart way its generic state 

f[A] = E (0^4^; (^'J^]) • • • <P':.n.'ShiAA])) , (7.14) 

where R is all the rest of the factorization, li to I4 are the four links converging 
at Uq- What we are going to do is to render invariant this small piece of graph, 
that is the one sourrounding this node. The idea is to group average, i.e. to take 
the action of a gauge transformation in uq and to integrate over all the possible 
such transformations: 



^ * JSU{2) ^ ^ 

(7.15) 
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where we have assumed that Uq is the source for all the four links, but it would 
be the same otherwise. It is clear that the result is invariant under the action 
of U on the node uq, since the Haar measure is such! Now, we can simply do 
this on each node of a graph and obtain the desired invariant basis. But let 
us investigate a bit more such a projection on the single node Uq: since the (f) 
functions are essentially Wigner matrices it is obvious that 

cl>i(gh) = D^{g)<l>l(h) , (7.16) 

and then the projected state looks like 



Isi 



{K.m'Sg) . . . D=:,^,{g)) • • • K^m'ShiAA]) . (7.17) 



ISU{2) 

Let us give a name to this operator 

P"° , , = /" dg ( D'' ,(g)... D'' , (g)) . (7.18) 

J SU (2) 

This is indeed a projection operator from the tensor product of the four repre- 
sentation spaces to its SU (2) invariant subspace, where the projection property 
pnopno _ pno jg ^-^^g Haar measure invariance, as it is easy to check. We 

can write it as 

pno . yh ^ yh ^ yh ^ yu _^ l^^[yh ^ yh ^ yjz ^ ■j/i4] ^ (^^.19) 

We can pick an orthonormal basis in Invfy-'^ (8) V^^ V^^ (g) V^'^] - let us call it 
\l) - and decompose P"° as 

P«o = ^|i)(i|. (7.20) 

Obviously the dimension of the invariant subspace depends on the spins that 
contribute to the node. The basis vectors i.e. an orthonormal basis of the 
invariant subspace of a tensor product of vector spaces, are usually known as 
intertwiner operators or simply intertwiners. 



We can do a very simple calculation just to explain how things works. Take a 3-valent 
node with spins Ji, J2j3- we have to deal with V^^ V^^ fS> V-'''-\ which is often written 
in the sloppy way ji j2 <Si js- Then we decompose it into the sum of irreducible 
representations, namely 

ji <X) j2 <8) h = (I ji - j2| © • • • ii + j2) <8) js , (7.21) 
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and so on, which is of course the well-known Clebsh-Gordon condition. Take three 
values, for example ji = 1/2, j2 = 1, js = 3/2, then the decomposition reads 

1 3 13 3 

(g)io=(©)(^=ie2eoeie2e3. (7.22) 

The invariant subspace is, by definition, the representation space. The dimension of 
Inv[ji(5Dj2(SDj3] is given by the multiplicity of the representation in the decomposition. 
It is not hard to see that for a 3-valent node there will be at most a one dimensional 
invariant subspace, and this is the case when the 'selection rule' 

i3G{|ii-i2|...Ji+j2} (7.23) 



holds. This means that for every 3-valent node that satisfy the selection rule (7.231, 
the projector into the invariant subspace is unique up to normalization. Actually it is 
just the (normalized) Wigner 3j-symbol |160) 

Every node with valence more than 3 can be decomposed into 3-valent contractions. 
For instance, a 4-valent intertwiner can be written as 

ai,02,a3,a4 _ ai,02,/3, a3,a4 or^ 

the label i is a spin representation, it ranges over the Clebsh-Gordon- allowed repre- 
sentations joining the four spins. Graphically is everything quite self-explanatory: you 
can see the 4-valent intertwiner as 




J4 

where i is sometimes said to be a virtual link, since it is a spin, but it is spans all the 
allowed spins. Of course one can join the four different links in more than a way. In 
this regard the following equality holds 

. . , V" 

-^(2j + l)A2 ^\ J (7.26) 




J2/ \ J4 



I refer to the appendix of |151| for more details on SU{2) recoupling theory. 
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In general, for a t;-valent node n, one can define with the same procedure of 
group averaging, a projection operator P", given by 

pn _ \ ^ Ov .ay * 27") 

m,\...my,n\...ny / j m\...my n\...ny ' v' '/ 

Olv 

where denotes the elements of the basis. This is clearly a generalization of 



equation (7.20). 



Thus, we have come to a result: an orthonormal basis of l-i^^^ is given by the 
following 

ST,{i},m[A] = (g)0^'(/ii[A]) ■ (g).„ , (7.28) 

/ n 

where / runs over the links and n overs the nodes of the graph F. The ■ stands 
for the contraction between the Wigner matrices and the intertwiner operators. 
These elements are manifestly SU{2) invariant, thanks to the contraction with 
the invariant intertwiner operators, and they are know as spin network state^ 
For future reference, and also to understand better the relation between the 
graph and the intertwiners defined on it, let us analyze a bit the situation of a 
4-valent node. 



7.4. The vector/difFeomorphism constraint. TL^^^^ 

We shall now deal with space diffeomorphisms. We proceed just as for H^^^, 
i.e. we first see how diffeomorphisms act on the vectors of "H^^, then we try 
to "average" over all the possible actions, obtaining a diff-invariant subspace. 
Here, however, one must be more careful, since diff-invariant functions surely 
won't be inside "H^^, because the orbits of the action of diffeomorphisms are 
not compact. It is just as if you want to constrain wave functions defined on a 
cylinder ijj{6,x), ip G L'^{S^ x M), to pgip = and pxijj = 0. This constraints say 
simply that ip cannot depend neither on x nor on 6, it is a constant. But while 
the integration on 6 gives no problem, since 5*^ is compact, the integration over 
X is not bounded, and the constrained states do not belong to L^. However it is 
always possible to choose a suitable dense subset $ C of test functions and 
define the constrained states as distributions on this space. Then one gets the 
Gelfand triple (see, e.g. |163| ) $ C C $*, where $* is the 'extension' of 
we are looking for. 

If is a diffeomorphism of S, then its action on a cylindrical function ^/'^j 
(7.1) is obvious 

Uom-yA^] = ^P^^^,j[A] ■ (7-29) 



^See |158l 1159] for interesting details. 
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Then the vector constraint can be rephrased as 

UD[(l)]i^ = iJ , (7.30) 

where, as we have just said, the solutions can be found only in distributional 
sense, namely Cyl C "Hkm C Cyl*, with ip e Cyl* the space of linear functionals 
on Cyl. Now it is not worth to go on in technicalities: one simply average the 
distributions on Cyl over all possible diffeomorphisms, thus obtaining only diff- 



invariant states. A look to (7.29) suggests that the resulting space is formed by 
the equivalence classes of graphs under spatial diffeomorphisms. These graphs 
are, mathematically speaking, /^nois|^ and the spin network states after this 
procedure are also known as s-knots. They are a basis of the space H^^. 

Summarizing/simplifying: We have first obtained the kinematical Hilbert 
space from a suitable definition of a measure on Cyl. Then we have seen that im- 
posing the Gauss constraint amounts to insert intertwiner operators at each node 
of the basis functions. Finally we have shown that the diff constraints simply 
say that we have to take the equivalence classes of the base graphs (knots). 



7.5. Quantization of the algebra and geometric 
operators 

We have been able to interpret and "solve" the Gauss constraint without actually 



quantize its corresponding operator in equation (6.43). Here we shall face the 



problem of quantizing the triad operator E'^ and the connection operator, or. 



better, their smeared version as introduced in section 6.4 It is rather simple: 



we have defined the Hilbert space "Hkm as (an opportune definition of) L'^{A), 



thus on then spin- network basis (7.28), the connection operator shall act by 



multiplication, i.e. - considering for simplicity the fundamental representation 

h^[A]K[A] = h^[A]K[A] . (7.31) 
The flux shall act by derivation 

MS)K[A] = -ih^ [ dV \ = ±zhjKAAhKAA] . (7.32) 

Js SAl{x{a)) 

Some remarks: the flux should be seen as a surface (namely the S surface) (see 



figure 7.1); the holonomy (i.e. the smeared connection) is represented by the 
path e. If the surface S and the path e does non cross each other then the action 
above is identically zero. If they do cross each other, then the path is "cut" in 



knot is an embedding of a circle in a 3d space, up to isotopies. See, e.g. 1147; for details. 
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Figure 7.1.: The flux surface S 'punctured' by an holonomy e = ei U 62- 



two sub-paths ei and 62 and the flux operator inserts an SU (2) generator among 
them. The sign in front depends on the relative orientation of 5* and e. Notice 
that recovering all the fundamental constants one gets an overall Planck length 
squared for Ei{S) (replacing the h), 



= y ^ ^ 10-33^m ; (7.33) 

this will be important to understand the scales we are working with. 

Now that we have understood how the fundamental variables of our theory 
work, once quantized as operator living in our Hilbert space, we can turn to really 
interesting questions about the structure that emerges from this framework of 
canonical Quantum Gravity. 

Actually it is quite simple to see that the area of a surface 5* embedded in our 
space manifold S can be written in a amazingly simple way in terms of triad 
variables, namely: 

A{S) = I da'da^^E^E^5'^nanb , (7.34) 
^ s 

but in our quantum theory the triad is an operator that acts on "Hfjj^, what shall 
become the area of a surface? 

First of all let us understand how the product of two triads acts on holonomies. 



With (7.32) is easy to see that 

E,{S)EjiS)K[A] = -ll^^K,[A]nTjK,[A] . (7.35) 

When the two fluxes are contracted one obtains the Casimir operator of the 
representation. In this case we have simply holonomies, i.e. fundamental repre- 
sentation, thus = TiT'^ = — |l2- Notice that the Casimir commutes with all 
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S 



Figure 7.2.: A surface S intersected by a generic spin network. You can see that 
refining again the regularization does not change the result, since 
each sub-cell has already at most one puncture. 



the group elements: this fundamental feature allows us to "recompose" the path 
e, namely 

E,{S)E\S)K[A] = -iyC^K[A] . (7.36) 

This means that the holonomy is an eigenstate of the square of fluxes! In a 
generic spin-j representation one gets 

%{S)E\S)D^{K[A]) = I'^^'jU + 1)D^{K[A]) , (7.37) 

since the Casimir reads = + l)l2j+i. 



Now let us turn back to the area (7.34). We have to write it in a way to 



include the smeared version of the triad, the flux, in order to be able to act with 
it on a general state vector. Thus we have to pick up a regularization of the area 



(7.34): We decompose the surface S into N 2-cells, and write the integral as the 



limit of the Riemann sum 



A{S) = lim An{S) , (7.38) 

N^oo 



with 

N 

Ar,{S) = Y, = VMSi)mSi) . (7.39) 
1=1 

Summarizing: This amounts to divide 5* in 2-cells, take the area of each of 
them Ai\i{S), and let go to infinity. It is just a way to regularize the integral 
over S. Next we have to deal with An{S). We shall see what kind of action the 



operator y Ei(S j) E"^ (S i) does on a generic state vector ipr- The trick is to take 
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the regularization sufficiently fine so that each Si intersects once and only once 
the embedded graph T. If so, we already know how this operator acts, namely 



Ei{Si)E^{Si) i/jr = iK^JvUv + 1) , (7.40) 

where p denotes the single link of F that puncture that specific 5/. Refining 
again the decomposition will have no consequences, since at most we have no 
intersection at all for some cells, but that gives a zero contribute to the area (see 



figure 7.2). Thus, at the end. 



A{S) ^Pr= Yl ^^1\/MJp + 1) . (7.41) 
pe5ur 

The great clues of this formula are: 

• the area operator is discrete, 

• its eigenfunctions are just the spin network states, 

• a part from the Immirzi parameter (which is a free parameter of LQG, 
even if there are some proposals to fix its value through the Black Hole 
entropy calculation |16H [IS2j) the area eigenvalues are of the scale of the 
Planck length squared, which is precisely what intuitively we expected. 

We shall see in a moment what a powerful intuitive picture these two key points 
provide. First, let us deal with the volume operator. We have to say that this 
issue contains some technicalities, for instance there are at present two distinct 
mathematically well defined volume operators. The distinctions come out in the 
regularization process. We shall not discuss it here, we merely present the results 
that are in agreement for both the versions, which - incidentally - are the truly 
interesting part of the story. The volume operator has the same two remarkable 
features of the area operator, namely it has discrete spectrum and it is diagonal 
in the spin network basis, and is of order 1^. 

Now we are ready to argue the picture that quantization is suggesting us: space 
geometry is discrete at the Planck scale. Each spin network is a polymer-like 
excitation of space: volume excitation being dual to the nodes of a spin network 
and area excitation dual to the links. The area excitation is proportional to the 
spin quantum number, attached to every link of a spin network, while volume 
ones are determined by the intertwiner space of the node to which it is dual. 
This interpretation has also the power of being fully background independent, 
space geometry is determined in a pure combinatorial way by the spin network 
state, which thus can be seen as a quantum geometry of space. Something which 
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is obvious, but it is worth some words, is that this is not a built-in discretization, 
as it for instance in the lattice formulation of Quantum Gravity: quantization 
is telling us that space has a discrete nature, whose building blocks are of the 
order of the Planck scale. 

Of course one can rise many objections to all the structure we have built so far, 
but it cannot be denied that this picture is extremely suggestive and powerful, 
and this is at the heart of Loop Quantum Gravity itself. 



7.6. The scalar constraint 

One step remains to be analyzed in order to complete the program of (canonical) 



Loop Quantum Gravity, i.e. to deal with the scalar constraint (6.43). Let us 
recall its expression, in a smeared form 

SiN) = ^ d^x A^^2= {^\F^b - 2(1 + l')Km) (7.42) 

= S^{N) - 2(1 + 7^)T(A^) . (7.43) 

The second line is a usual shorthand notation to separate the so-called euclidean 
contribution 

S^'iN) = [ d^x N e'\F^, (7.44) 

from the rest 

The non-linearity if this expression is really awkward, in order to be able to 
quantize it properly. Notice that this is just the same problem one had before 
the loop representation. And it is no mystery that this problem remains the real 
big problem of canonical Quantum Gravity. 

However, the rich structure we have built so far, surely has made possible many 
important steps towards a better comprehension of the quantum scalar constraint 
and (particularly) of its action on the kinematical state space. We shall not deal 
here with the many complicated open issues in this respect |165[ 1166] but we 
want at least to explain what has and what has not been achieved by now. 

One important improvement has been put forward by Thiemann |164j . He 
observed that, if one introduces 

K = [ KiEt , (7.46) 
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then the following identities hold 



{S^{l),V} 



K 



7 



3/2 



(7.47) 
(7 



4 



Vdet(E) 7 



(7.49) 



with V = j \J<lei{E). This is actually very useful, since (7.43) becomes 



(7.50) 



and 



N 



T{E) = J d'x -,e'^'%AA:dS'^{l),V}}{Al{S^{l),V}}{A'^,V} . (7.51) 

Now the trick proceeds by rewriting both the connection and the curvature in 
terms of holonomies and so come to an expression that involves only the volume 
and the holonomies. This expression can be transformed into an operator, since 
we already know how h and V act on a generic state of our Hilbert space. We do 
not list here all the passages, be sufficient to say that thanks to the well known 



(7.52) 



where Ca is a path along the direction - one can express both A and F in 



terms of h. The integral in (7.43) must be regularized as for the area and volume 



operator, thus the space must be decomposed into a cellular decomposition, a 
triangulation, i.e. a collection of tetrahedra bound together. We give here the 



regularized expression for the euclidean part of (7.43) 



S^{N) = lim^Ar/eV6cTr(F„6(A){A„r}) = 

lim T Nie^'" kh^iM] - hj [A])h;/[A]{h-J[A],V}] , (7.53) 



is the volume of a cell, Uah is an infinitesimal loop in the a6-plane, around the 
face of the J-th cell; while e\ is an infinitesimal path along the a-direction, along 



an edge of the J-th cell (see figure 7.3). Notice that the dependence on the cell 



scale disappears in terms of holonomies. 
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Figure 7.3.: One tetrahedron of the regularization cell-decomposition. One face- 
loop and one edge-path are shown. 



H 



(7.55) 



Figure 7.4.: Typical action of the hamiltonian operator on a node of spin net- 
work. In the framework we have sketched s = 1/2, but we shall see 
later that links with generic spin label could be consistent as well, 
and this is one of the (many) ambiguities in the quantization of the 
scalar constraint. 



Now we can formally 'quantize' this operator (7.53) by putting an hat over 



the volume and the holonomies and see how it acts on a spin network: 



(7.54) 



This operator is well-defined and it is very easy to see how it acts on "H^^: it 
has the property of the volume operator to act only on nodes of spin networks. 



while the holonomies in (7.54) modifies the spin network by creating new links 
in the fundamental representation (spin 1/2) (see figure 7.4). 



(We have dealt only with the euclidean part, but a similar analysis can be 
done for the T part of (7.43) as well, and the following results still hold.) This 



way of acting is at the real clue of Loop Quantum Gravity with respect to the 
quantization of the scalar constraint. For example, take all the spin network 
states with no nodes: they are just Wilson loops, i.e. the trace of the holonomy 
around a closed path. All the Wilson loops, irrespectively from the spin, are 
solutions of the quantum scalar constraint equation! Actually this simple but 
important fact is what started the interest in Loop Quantum Gravity: it provides 
a set of explicit solutions of the quantum theory of gravity. Moreover, equation 
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(7.54) tells us that the action on nodes is rather peculiar: it creates special links, 
sometimes called exceptional links. They are special in the sense that the new 
nodes they carry are exactly of zero volume, thus being 'invisible' to a further 
action of the hamiltonian operator. This allow the following picture: a generic 
solution of the quantum scalar constraint is labeled by graphs with nodes of the 
following kind 




a 




+ ... + /? 




+ 




(7.56) 



also known as dressed nodes, i.e. with infinite superpositions of exceptional 
edges. The Q label stands for the collection of weights of the superposition. 

Up till now the good points. Now let us spend some words for the big problems 
still to be solved: 



• the action of the hamiltonian operator that we have described above is 
said to be ultra-local [167J . This feature has raised some concerns whether 
this kind of theory is at all capable of reproducing general relativity in the 
classical limit |167| . 

• there is a large degree of ambiguity in the definition of the quantum scalar 
constraint! One of such ambiguities arises from using the holonomies in 
the fundamental representation: we could have used any representation, 
thus the action of the scalar operator would have been different, creating 
links with arbitrary spin. This doesn't affect the ultra-locality issue, but 
certainly we would end up with an infinite set of dynamical theories, all 
different. Another ambiguity is in the regularization scheme: the scalar 
quantum operator is regularization-dependent! 



there is still a mathematical problem in analyzing thoroughly the limit in 
(7.54) (and the analogous one for the T part of (7.43)), the limit must exist 
and should give a well defined operator acting on the space of s-knots. 



These kind of ambiguities have stimulated both the research into a better 
comprehension of the scalar constraint itself, for example the Master constraint 
approach |168j ^ and new paths of research, the most important of which is the 
spinfoam formalism, to which we dedicate the next sections of this chapter. 

^See also the book by Thiemann |166] for a more self-contained and thorough presentation 
of the entire scalar constraint problem. 
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7.7. Concluding remarks 

I have tried in this chapter to review the very basics of LQG. Of course many 
are the arguments I had to omit not to go beyond the goals I have in mind. Here 
I show a schematic hst of the most important issues I have not treated, but that 
certainly deserve attention, in order to understand the global importance of such 
theory: 

• First of all the scalar constraints definitely deserves more attention. I refer 
to [166j . 

• Black Hole entropy. The idea is to quantize a sector of the theory contain- 
ing an isolated horizon and then to count the number of physical states 
compatible with a given macroscopic area of the horizon. See |169j for 
detailed calculations and |151] for discussions. 

• Loop Quantum Cosmology. This is the cosmological sector of LQG, which 
has rapidly developed since its birth in 2000. There are indeed quite in- 
teresting and stimulating results, particularly concerning the initial sin- 
gularity of GR, i.e. precisely where we expected Quantum Gravity to 
tell something. I refer to the review |174) and to the papers [170]- [T73] . 
Recently also spinfoam calculations of cosmological problems have been 
performed [HSl [HSl [177]. 
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We have up till now described what "happens quantum mechanically" on spatial 
hypersurfaces, and we have argued the difficulties in dealing with the evolution 
of these quantum states in a timelike direction, difficulties connected with the 
hamiltonian operator. However - with in mind the goal of some kind of sum- 
over-histories formulation of the dynamics - we can take a quantum 3-geometry, 
for simplicity a single spin network state (i.e. an eigenstate of area and volume) 
|s') and imagine it evolving in a timelike direction. Imagine its "world-sheet": 
every node would sweep timelike paths, that we call edges, every link would 
sweep faces and there will be some spacetime points in which two or three edges 
meet, that is events in which, for instance, one edge split in two (or whatever) 
or two or more edges converge in one, etc. . . We call these points vertices and 
they are the points in which the original quantum state \s') changes by means 
of the dynamics, i.e. by action of the hamiltonian operatoi^ Time evolution 
has thus produced a 2-complex with vertices, edges and faces. This is a dressed 
2-complex, since the coloring of the spin network will induce a coloring in the 
2-complex as well: irreducible representations to faces and intertwiners to edges. 
This dressed 2-complex is what is called a spinfoam. 

The spinfoam approach to (Loop) Quantum Gravity began with the work of 
Reiseberger and Rovelli |178[ 1179] , which first put out the idea of defining quan- 
tum amplitudes as sum over histories starting from the hamiltonian formulation 
(see next section for more details). Actually, work had already been done in this 
respect (see for example [ISOJ), but not in the LQG framework. The general 
idea of a spinfoam derived from a path-integral discretization procedure (see 
section 8.2), was instead proposed by Baez in the seminal paper |181) (in which 
the name "spin foam" was first used) and in the lectures |182] . which are still a 
very good 'beginner's guide' to spinfoam models, perhaps one of the best in the 
literature. For other good, but a bit dated, reviews see for example |15H I152J: for 
more up-to-date (but a bit more technical) resumes see instead [154[ 11551 il83j . 



^thus we have already an idea about how two quantum states differing from a single vertex 



i.e. a action on a node by the hamiltonian operator - should be like, see section 7.6 
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Figure 8.1.: Example of spinfoam: evolution from a spin network s' to s passing 
through an intermediate spin network through two vertices, i.e. ac- 
tions of the hamiltonian operators. You can see the edges (evolution 
of nodes of s'), faces (evolution of links of s' - one of the faces has 
been colored as an example) and vertices. 

8.1. Sum-over-histories from hamiltonian 
formulation 

We can now make a rather heuristic reasoning trying to "derive" a path inte- 
gral formulation from the hamiltonian formulation, just as one does in standard 
Quantum Mechanics. I want to stress again the heuristic approach of what fol- 
lows, that is mainly taken by [151j (but you can find more detailed discussions 
in |184) ). In the next section we shall see a far more rigorous (and alternative) 
definition of the sum over histories approach to the Quantum Gravity dynamics. 
Rigorous proofs that the two derivations are actually the same thing have been 
found only for the three dimensional case [185J. 

The spirit of the spinfoam approach is to try a Feynman-like procedure in 
a gravitational and background independent context. In standard Quantum 
Mechanics the Feynman idea is somehow summarized in the following expression 

{y,t'\x,t)^ L ^Pge^^M . (8.1) 

On the left hand side we have the scalar product (transition amplitude) between 
two different position eigenstates at two different times. The same amplitude 
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can be written in the Heisenberg picture as {y\U {t, t') \ x), with U {t, t') evolution 
operator form time t to t'. In a gravitational context one should then give sense 
to the following expression 

j Vg,, e*^-^[^l (8.2) 

or, to be more specific, 

A{g,g')r^ [ _ Vg,, e^'-"^^^ , (8.3) 

J9\t=i—9 
9\t=o=9 

which formally represents the transition amplitude from a space with metric g 
to one with metric g'. The specific values of the time label are irrelevant for the 
diffeomorphism invariance holds. There is one big trouble when approaching an 



integral like the one in (8.3): we do not know a non perturbative definition of the 
measure Vg^,^ and perturb at ively we know that the theory is non-renormalizable 
(and moreover in this case we should break the background independence). 



In order to give a concrete significance to the expression (8.3) one usually 
starts from the canonical formulation, which is just what is done for the actual 
definition of path integrals in all the others fields. Let's for a moment resume 
what is the standard way to define the sum- over- paths procedure. The idea is 
to take the time evolution operator e~*^* and to have it acting step by step (in 
each step you have simply the action of if on a state, which you can calculate), 
then taking the limit of this time step going to zero. This gives a mathematically 
precise definition of the Feynman path integral representation of the propagator. 

Unfortunately, we are not able to do the same thing in the Quantum Gravity 
context. The problem is again the hamiltonian operator: we do not know how 



to quantize it in a proper way, as was pointed out in 7^ Thus, we cannot follow 
this route pretending to obtain new rigorous insights. 

However, we can still derive some basic properties that the Quantum Gravity 
path integral should satisfy. We shall see that this will be enough to find at least 
a form for Quantum Gravity transition amplitudes. Then, in the following sec- 
tion, we shall see an alternative way to follow, somewhat dual to the derivation 
from the canonical formulation, which is called "spinfoam" approach. 



Consider the integral (8.3): between which states the amplitude should be 
computed? It should be computed between eigenstates of the three-geometry, 
i.e. with states with a definite 3d metric: but these are nothing but the spin- 
network states. Thus what we would like to calculate is actually something 
like 



{s I s') 



phys 



{s\A\ s') 



kin 



(8.4) 
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Let us try to explain thoroughly this last expression. |s) is the usual spin- network 



state (7.28), or, better, the s-knot states, i.e. with basis graph an equivalence 



class of graphs under spatial diffs (see section 7.4). This is also the reason for 
the subscript "kin", to stress that the spin networks belong to the kinematical 
Hilbert space T^km, but not to the physical Hilbert space. So A (for 'amplitude') 
represents the projector on the kernel of the hamiltonian operator H, i.e. the 
projector onto the physical Hilbert space. If we assume for simplicity that the 
hamiltonian H has a non-negative spectrum, then we could (formally) write 

A = lim e"^* , (8.5) 

t— >oo 

indeed, if \n) is a basis that diagonalizes H (with eigenvalues En), then 

A = lim V |n)e"^"*(n| = V 5o,i?Jn)(n| , (8.6) 



namely A projects onto the lowest-energy subspace, i.e. the kernel of if, if we 
assume a non-negative spectrum. 

Proceeding with these formal manipulations, we can also write 



A = lim TTe-^(^)* = lim e"^^'^ , (8.7) 

X 

hence 



A{s,s') = lim (s|e-^^'" | s') 



kin ■ 



Quite loosely, if we want the propagator to be 4d diff invariant, then the limit 
is irrelevant, so 

A{s, s') = (s I ^0 ^* ^ I s')ki„ . (8.9) 

Now we can split this expression by inserting identities in the form 1 = 
obtaining something like 

A{s, s') = hm I e'^''^ I ^^)kin(^^ I e--^"'^ ^^^^^^ | .iv-i)kin (8.10) 

A/ — ' * 



7V-s>oo 



kin 



We can see then, that the transition amplitude between two spin-network stats, 
can be expressed as a sum, as follows 

A{s,s') = Y,M^) , (8-11) 
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i.e. as a sum over histories a of spin network. A history cr is a discrete sequence 
of spin-network 

a = {s,sn, ■ ■ ■ ,si,s') . (8.12) 

Moreover the amphtude is just a product of the amphtudes between the single 
steps in the history 

A{a) = l[A,{a) , (8.13) 

V 

where we have labeled with v each of the above said single steps (this terminology 
will be clarified in a moment). Now, to take a further step, let us recall that 
the hamiltonian operator acts only on the nodes of the spin network. Thus the 



single amplitude in (8.13) is non- vanishing only between spin- network that 
differ at a node by the action of H. 

A history of spin-networks a is what is called a spinfoam, and it is exactly the 
"world-sheet" of the spin network, as pointed out at the very beginning of this 
section: it is the time evolution of a spin network state. 

As a conclusion of this brief passage I want to recall again two things: 1. the 
heuristic value of the former passages and 2. the fact that a rigorous derivation 
of a sum-over-histories formula for Quantum Gravity transition amplitude can 
be given in the three dimensional case |185j (and moreover it matches the result 
of the following section). 



8.2. Path integral discretization: BF theory 

Now we would like to introduce the spinfoams formalism in a different (and 
somewhat clearer) way. We shall try to "discretize" the path-integral itself. Or, 
better, we shall discretize spacetime with triangulations a la Regge \186\ 1187] , 
and try to read out how the path integral of general relativity can be adapted 
on this triangulation. 

This is an approach somewhat dual to the hamiltonian derivation of the path 
integral: there one take the physical inner product defined in terms of an evo- 
lution operator (the projector into the hamiltonian kernel) and "decompose" it 
by inserting identity resolutions, de facto giving a discrete definition of the path 
integral. The present approach is the other way around, i.e. one defines the am- 
plitude as the path integral of (an appropriate form of) the General Relativity 
action, and discretize this integral trying to give it a rigorous meaning. For a 
good introduction to this approach see [182]. 

Remark. I point out from the very beginning that the discretization we are 
imposing to define the path integral is not at all as the one of, say, relativity on 
the lattice. Indeed, there one uses the lattice just as a regularization procedure. 
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to be removed at the end by an appropriate continuous limit. Here instead, 
quantization tells us that the fundamental theory is discrete! So we are justified 
in our discretization procedure, and no continuous limit is to be done. Another 
problem will be the triangulation dependence of our final results, but this is a 
completely different issue. 

We shall start with BF theory, which is actually a trivial theory, but it is up 
till now the only theory in which the spinfoam approach (=discretization of path 
integral) can be completely gone through. 

To set up a general BF theory we need a principal bundle - let's call it P 
- with base space the spacetime manifold A4, fiber a gauge Lie group G. The 
basics fields in the theory are a connection A on P and an ad(P)-valued {n — 2)- 
form B on A4. Here ad(P) is the associated vector bundle to P via the adjoint 
action of the group G on its Lie algebra. However, we can forget for the time 
being all these technicalities and just see how things work in calculations. The 
lagrangian of BF theory is the defined as follows 

C = Tr(P A F) , (8.14) 

where, as usual, F is the curvature of the connection A. Now it easy to see that 
this theory is (at least classically) trivial, indeed the equations of motion state 
that 

P = ; d^P = . (8.15) 

The first says that the manifold is fiat, there are no local degrees of freedom. 
Indeed the second says that the parallel transport rule is trivial. Now, take the 
simplest path integral you can imagine for this theory, i.e. the partition function 

Z{M) = J VAVBe'^M^'^^"^^ (8.16) 

VA5F , 

where we have integrated out the B field, using its equation of motion P = 0. 
This expression has, by itself, no rigorous meaning. To give a sense to the above 
formal expression the next step - as stated at the beginning of this section - is 
the triangulation of the manifold M. and the definition of a discretized version 



of (8.16). 



Let's recall that a triangulation of a (sufficiently smooth) manifold is obtained 
by a discretization by means of simplices, precisely n-simplices. Recall that in 
the sum-over-histories framework, we argued that spacetime should be seen as 
the dual to what we called a spinfoam, i.e. dual to spin networks world- sheets 
(just as 3- volumes are dual to nodes of a spin network on space-hypersurfaces). 
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This picture will somehow guide us in the discretization of the path integral, 
indeed we shall define the discrete versions of i?'s and A's on the dual of the 
triangulation. 

We call A such a triangulation. Now we take the so called dual 2-skeleton 
of this triangulation. It is built in this way: put a vertex at the center of 
each n-simplex; an edge intersecting each {n — l)-simplex and one (polygonal) 
face intersecting each (n — 2)-simplex. We could go on, defining dual volumes 
intersecting (n — 3)-simplices, and so on till defining a dual n-complex to each 
point (0-simplex) of A, thus building the dual triangulation A*. But we shall 
need only its two dimensional subset, the skeleton. 

A 2-skeleton C A* 

n-simplex vertex v 

(n-l)-simplex edge e 

(n-2)-simplex face / 

For example, if n = 3, then the triangulation is made up by tetrahedra glued 
along faces. For each tetrahedron we have one vertex of the dual skeleton, 4 
edges (one for each triangle) and 6 faces (one for each edge of the tetrahedron). 

A 2-skeleton C A* 
tetrehedra vertex v 

triangle edge e 

edge (of a triangle) face / 

While, for n = 4 - the interesting case - we have 

A 2-skeleton C A* 

4-simplex vertex v 

tetrehedra edge e 

triangle face / 

Now we want to define our BF field theory on this dual skeleton, i.e. the 
points on the manifold will be "replaced" by the vertices, edges and faces of the 
dual skeleton. 

A connection is a one form (precisely a g- valued one form), thus it is natural 
to associate it to a n — 1 dimensional object in spacetime. And this is just what 
we have called edge of the dual skeletoij^ Thus our connection is a prescription 
to associate a group element to each edge of the skeleton (think of it as the 

^In this sense, the discretization defined on the dual of A seems a rather natural choice. 
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holonomy of A along that edge). If the connection has to be flat (as is the case 
in BF theory) then we should want that 

ge^f ■ ■■gej.f = 1 , (8.17) 

where: g stands for a group element; e/ are the edges that surround the face /, 
and there are of them (i.e. we have called the number of edges of the face 
This formula thus states that the holonomy around each face is the identity 
and this is the flatness of the manifold. The B field is a g-valued 2-form on Ai 
and we discretize it as a map assigning to each {n — 2)-simplex - i.e. to each 
face / of the skeleton - a element By . 

Equation ( |8.16[ ) in its discrete version, as just prescribed, becomes 

Z^{M) = [Hdg, HdBfc'^^^^'f^f^ = fUdg, J] ^(^e,/ • • • ^7e./) ; (8.18) 

where with S and J-" are denoted the set of all the edges and of all the faces of 
the 2-skeleton, respectively. I have put a subscript to stress that this discrete 
definition is, in general, triangulation -dependent. In these last equalities we have 
performed the discretized version the integral over B. A remark on the expression 
TT[BfUf]: Uf = ge^f ■ ■ ■ gei^f is the (discrete) holonomy around the / face; using 



the identity [// ~ 1 + where G we get the exponent of (8.18) 



The importance of this last formula (8.18) is that it has a precise and definite 



meaning, and it is no more a purely formal prescription. It is actually the first 
spinfoam path integral we encounter in this thesis. Here it is written in terms 
of group variables. Now we shall see how to evaluate the integrals and pass to 
the spin/intertwiner representation, more suitable for an intuitive grasp and for 
a direct link wit spin networks, but completely equivalent. 

Let us use the following decomposition of the group delta function (Peter- Weyl 
decomposition) 

5{g) = Yl dim(p)Tr(p(^7)) , (8.19) 

p6lrrep(G) 

i.e. a sum over the all irreducible representations of the group of the character 
of the representation weighted by its dimension. Thus we have 

Z^{M) = Yl fU'^ae ndim(p/)Tr(p;((7e,/...^7e,/)) ; (8.20) 

p:J"-5-Irrep(G) e&S f&F 

where the sum is over all the manners to associate to the faces of the 2-skeleton 
irreducible representation of the group (operation usually called "coloring"). We 



■^Notice that the dual faces can have an arbitrary number of edges, it depends on how the 
triangulation is done and on the manifold M.. 
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can go further. It is crucial the following identity; 



LL* 



dim(pi) 



>12 



(8.21) 



if pi ~ p2* and otherwise. Maybe it is useful to write this identity graphically 
like this 



d^ ( .9 



Pi 



P2 




(8.22) 



A generalization of the above formula is the following: 

P2 



dg 



pi 



P3 




(8.23) 



P2 



and so on, with arbitrary number of representations involved. It is important 
to catch the meaning of these last formulae. Actually it is quite simple: the left 
hand side is a group averaging of the product of a certain number of represen- 
tations, thus it belongs to the invariant subspace of the tensor product of those 
representations. The right hand side is obviously a projector onto the invariant 
subspace: the two sides are the same thing. 

Let's go back to our BF theory. Let's make the things as simple as possible 
first: n = 2, our spacetime is a surface. Thus it is triangulated by triangles and 



the dual 2-skeleton is made by polygons. A look to figure [872] will certainly clarify 
what's going on. It is easy to see that each group element is shared by two faces 
(which is the same thing of saying that each edge is common to two and only 
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Figure 8.2.: (Part of the) triangulation of the manifold (dashed hnes) and its 
dual 2-skeleton (thick lines). One face / has been labeled with its 
group elements g/s on the edges. 



two faces). In our path integral formula (8.20) we thus see that each integration 
actually concerns the product of two representations, precisely the two faces 
that share the edge of that group element. Thus we can use formula (8.21) 



which implies that the integral is non-vanishing only if all the representations 
(i.e. on all the faces) are equal. The integral can then be evaluated and the final 
result is 



xiM) 



(8.24) 



pSlrrep(G) 



where = \^\ ~ 1^1 + l-^l is the Euler characteristic of the manifold, and it 

is a topological invariant quantity. For a manifold of genus g it equals 2 — 2g. 
Thus our partition function for a 2d BF theory converges only for surfaces with 
genus greater than 1 (thus it diverges for sphere and torus). See section 8.5 for 
some discussions about divergences in spinfoam models. 

Let's step up to dimension 3. In this case the triangulation of Ai is made up 
by tetrahedra. Again, it is worthwhile to have a look to figures 8^ and 8^ We 
have shown an edge of the dual skeleton, which is shared by three faces. Thus 
the typical integral will be now of the kind of the left hand side of equation 
(8.23), and thus the integrations split the edge in the sum of intertwiners of the 
kind shown in figure 8^ Now, it is not difficult to see that the intertwiners of a 
single tetrahedron combine to form a tetrahedron as well (a dual one) labeled by 
the representations of its 6 edges and by the intertwiners on its vertices. Thus 
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Figure 8.3.: A tetrahedron of a 3d triangulation is shown. The thick hnes are 
the four edges, one for each triangle. One (typical) face has been 
colored, the one corresponding to the dotted line. 




Figure 8.4.: (on the left) A tetrahedron of a triangulation of a 3d manifold. The 
three thick lines represent actually a single edge, the one dual to 
the triangle punctured, which is shared by three faces (one per each 
side of the triangle). On the right figure, the integration over the g 



group element has been evaluated by means of equation (8.23) 
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we have an explicit formula for the partition function, i.e. 



4^^(-M)=EEndini(p/)n 




(8.25) 



The first sum is over p : ^ Irrep(G'), i.e. a labeling of faces with G- 
representations, while the second on labeling of edges with intertwiners. The 
graphic "atom" in this formula is just the contraction of four 3-valent intertwin- 
ers, and depends obviously on 6 spins. If G = SU{2), it is nothing but the 
Wigner 6j-symbol. You can equivalently see this atom as a spin network (with 6 
links and 4 nodes) in a tetrahedron-like pattern, with all the six group elements 
put on the identity, which is often called evaluation of a spin network (recall 



what a spin network is precisely, see equation (7.28)). 

The trick is the same in 4d, mutatis mutandis. In this case one has 4 faces 



that share an edge, thus we shall use the formula (8.23) with 4 representations. 



In 4d the triangulation is made by a 4-simplex, and the integrations over the 
groups gives the contraction of intertwiners in a 4-simplex pattern. The final 
formula is 



4^H-M)=EEndi-(p/)n 




].26) 



The atom in this formula is the contraction of five 4-valent intertwiners in a 4- 
simplex like pattern. Incidentally, notice that the atom is precisely a 4-simplex 
(its projection on 2 dimensions). 

Notice that one can always split a 4-valent intertwiner into two 3-valent one^ 
thus one would end with a trivalent spin-network with 15 spins. If the gauge 
group is SU{2), this is called the '15j-symbol'. 



See section 7.3 and in particular equation (7.261. 
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Remark. Take a 3d model of gravity, i.e. 3d Riemannian General Relativity. 
Ponzano and Regge, in 1968, showed |188] that the SU{2)-BF path integral 
(see equation ( 8.25[ )) has the correct semi-classical limit, in the following sense: 



There is a discretized action (called the Regge action), which is a sum over all 
the tetrahedra with which we have triangulated the 3d manifold of 

SR = J2^e0e, (8.27) 

e 

i.e. the sum over the 6 edges of each tetrahedron of the product of the length 
of that edge and the dihedral angle (the angle between the two normals to the 
faces incident to that edge) - that is an approximation of the integral of the 
Ricci curvature, i.e. an approximation of the Einstein-Hilbert action. Now, the 



tetrahedral spin-network of equation (8.25), with gauge group SU{2), has the 
following asymptotics for large spins 



3nV 



cos 



{Sr + D , (8.28) 



where V is the volume of the tetrahedron and = je + 1/2. 

This is wonderful, since for large spins - which means for scales much grater 
than the Planck length - one exactly recovers something like e*"^ with the right 
General Relativity action This means that, at least in this simple and phys- 
ically un-interesting case, our discretization procedure is consistent, and - re- 
assured by this consistency - we feel more confident with the four dimensional 
case as well. 

The BF theory 3d model is commonly known as the Ponzano-Regge model. 

Up till now, we have dealt only with manifold without boundary, i.e. - to use 
an hamiltonian jargon - we have dealt only with vacuum-to- vacuum transition 



amplitudes. In order to define something like (s | s')^^^^ (see (8.4)), we have to 
consider triangulation of manifold with boundary as well. It is quite natural, 
but under this apparent simplicity many subtleties hide, so one must be careful. 
The simple part is: take a manifold without boundary and cut it in a spacelike 
direction. You get two manifold with spacelike boundary. Thinking in terms of 
the spinfoam, that is in terms of the skeleton of the triangulation, the boundary 
is made of links, results from the cutting of faces, and nodes, cut of edges. The 
reader has surely recognized that we are actually doing the inverse procedure 
with respect to world-sheet sweeping of a spinfoam from a spin-network (as in 



the hamiltonian-to-spinfoam approach 8.1) , i.e. spinfoam sectioning. 



^Actually one gets just the real part of the exponential. This is a technicality, roughly 
speaking the reason is that given the lengths of the edges of a tetrahedron, we still can 
rotate and reflect the tetrahedron. 
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The variables of the cut edges and faces will be variables of the boundary - 
be them group elements or spin/intertwiner (or whatever), depending on the 
representation we want to use - not to be summed/integrated over. All the con- 
struction is exactly the same, with some variables fixed, the boundary variables. 

Now that we have the notion of spinfoam with boundary, we can take the 
following picture: take a spinfoam (with or without boundary) and imagine to 
surround each vertex v with a little 3-sphere. You will get an ensamble of little 
bubbles - sometimes called atoms - each representing a little spinfoam formed 
by a single vertex with a boundary tp^, which is a two dimensional graph with 



links and nodes representing the vertex. Look to formula (8.25) or (8.26), this 
bubble is just the symbol on the far right, the vertex amplitud^ It is now clear 
(if it wasn't already) that a spinfoam amplitude is just a sum over labellings of 
product of amplitudes: the face amplitude (about which we will have much more 
to say) and, most important, the vertex amplitude. 

This is quite often taken as a definition of a spinfoam amplitude. 

The subtleties about boundary states arise in the Quantum Gravity context, 
i.e. when one tries to define quantum transition amplitude for general relativity 
and not for a generic G-BF theory. Indeed in that case, the boundary has to 
represent a state of "Hkin, the kinematical Hilbert space, i.e. SU (2) spin network 
states. This will be a key point. 

Before going on to explain what happens in the (interesting) case of gravity, 
I want to focus on another little but useful issue in BF theory. I have said that 



in formulae (8.25), (8.26), the term on the far right is actually the amplitude of 
a little spinfoam atom. What is like the amplitude of a spinfoam atom (i.e. the 
vertex amplitude) written in the terms of integral over the group (in the sense of 



equation (8.16)? To really understand this (and to really understand the vertex 
amplitude in general), one has to think at what happen when you cut away with 
a, {n — l)-sphere from the dual skeleton of an ra-triangulation. Let us focus on 
n = 3 and = 4. Actually, it is far more easy to think a little about it and 
catch it by oneself, rather than long and cumbersome (and useless) explanations. 
I just give a few hints: take a tetrahedron, and draw the dual. There will be 
one vertex, four edges departing and 6 faces. Cutting this dual will give another 



curved) tetrahedron graph (see figure 8.5), whose links are just the cut of the 



six faces of the spinfoam, and whose nodes are the cut of the edges. In the same 
manner, in a 4-simplex triangulation, when you cut away a 3 — sphere around 
a vertex you get another 4-simplex whose edges and nodes are the cut of faces 
and edges of the spinfoam. 



^Incidentally, this "bubble picture" is what suggested the name "spinfoam" in the early days 
of this theory. 
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64 



Figure 8.5.: (on the left) atom 'bubble' of a 3d triangulation: dotted line for 
the A-tetrahedron, thick lines for the edges of A*, colored lines 
for the "cut" tetrahedron. On the right the same graph for a 4d 
triangulation. 



Thus, one can write the ampltiude simply as 
MU^j) = / n n ^(9^U^,g-') = [U'^gell Sig^Ui^fgO , (8.29) 

i i<j e f 

where Uij = l,...,5;i < j) or Ui are the four(3d)/ten(4d) external vari- 
ables attached to the links of the boundary graph, the notation is quite self- 
explanatory. The integrals are of course over G. If one decompose the delta into 
representations and performs the group integrals - just as we have done before 
for the partition function - one ends precisely with 

(p,, .„) = 6j{pi, 6e) , At\pi, .„) = 15j(p/, ie) , (8.30) 

i.e. the vertex amplitude in the spin/intertwiner representations (i.e. if we take 
as boundary variables irreducibles on links and intertwiners on nodes) is - as it 
had to be - the tetrahedron-like or 4-simplex-like intertwiner contraction. 



8.3. Spinfoam models for Quantum Gravity 

Let us now come to the interesting part of the story: the spinfoam 'machinery' 
applied to General Relativity, rather than to the (trivial) BF theory. This is 
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still a very open issue, thus everything that follows is in some degree work in 
progress. 

The first fundamental step is to recognize that GR can be cast into "BF the- 
ory + constraints". Thus the general guideline will be to treat the path-integral 
discretization d la BF theory, and then impose the constraints. It is not surpris- 
ing that all the big troubles will be in that second step. But let us do everything 
at its proper time. 



8.3.1. The Barrett-Crane model 

The first spinfoam model for quantum gravity is the Barrett-Crane model [189]. 
Take an SO{A)-BF theory in 4d. Its action can be written (see ( 8.14[ )) 



S[B,uj] = I BjjAF'\u) , (8.31) 



where F is the curvature of the connection u. If we replace 

1 



Bij = -euKLc"" A , (8.32) 



we get precisely Einstein-Hilbert action, in the form of equation (6.31). Thus, 
we could formally write the action for General Relativity as 

S[B,u,<P,fi] = j {BjjAF'-'{u) + <PjjklB'' AB'''^ + fie'-'''''<PjjKL) • (8.33) 

Let us check it: is a 0-form while is a 4- form, they both have no dynamics 
and thus, classically, they just give constraints, namely 

e''''''<PijKL = , (8.34) 
i.e. (pijKL = —(pjiKL = —4>iJLK = (pKLij, for the variation with respect to /i. 



while for 0, using (8.34) 



e^'^'P'^BiiBfJ' = e e^^^^ , (8.35) 



where e = ^^euKLB^J^B^J^e^'^'"' . This last expression (8.35) contains 20 equa- 
tions - often called simplicity constraints - one for each independent component 
of 0. They constrain 20 of the 36 independent components of the B field. The 



solutions to (8.35) are of two kinds 



S/j = ±^ewLe^Ae^ ; 5" = ±e^ A ; (8.36) 
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in terms of the 16 degrees of freedom of the tetrad field e^. Only the first type of 
solution gives the GR action (6.31 while the second give a topological (trivial) 



action. This must be remembered when to quantize the action 



For the following, it is useful to rewrite the constraints (8.35) as 



■ B^j^a = , 



B*^,y ■ Bar 



±2e 



(8.37a) 
^8.37b) 
(8.37c) 



with fiuar all different, and with e = f^e^uardx^ A dx'^ A dx"" A dx"^. The first two 
constraints in these form are also known as diagonal and off- diagonal simplicity 
constraints. 

Remark. Notice that we have taken an 5*0(4) connection - or, better, its double 
cover Spin{4) - thus this is an euclidean version of General Relativity. To 
obtain the actual GR one should take an SL{2, C) (the double cover of 5*0(3, 1)) 
BF theory. The problem is the non compactness of SL{2, C) that requires special 
care in handling the (divergent) integrals over the group. It is a general custom 
in Quantum Gravity to investigate the simpler euclidean case first. 



Equation (8.33) is our GR "BF + constraints" action. Now the goal is to 



take the BF spinfoam model and analyze what kind of consequences have the 
imposition of constraints and, firstly, to understand how to impose them. The 
intuitive idea is that the constraints should restrict in some (interesting) sense 



the spinfoam sum in (8.26). In particular the constraints must impose that on 



the boundary of the spinfoam one gets the kinematical Hilbert space of LQG, and 
so they must somehow restrict 5*0(4) to SU{2) on the boundary. However, the 
imposition of the constraints is actually a rather cumbersome issue, to be dealt 
with very carefully. I shall sketch here briefly the key points. The interested 
reader is encouraged to read the review on Barrett-Crane model |191] . However, 
the issue of constraints imposition will be gone through more carefully in the 



case of the Engle, Pereira, Rovelli, Livine (EPRL) model (see section 8.3.2 for 
details and references). 

The idea is to take the following path integral 



j VAVBS[B - (e A e)*] exp (^i j B A F[A] 



Z^^{M) = I VAVB5[B - (e A e)*] exp ( z / B A F[A] ) . (8.38) 
Naively speaking, this means that one must restrict the integral to those config- 



urations satisfying the delta function in (8.38). This turns out to be the case if 



^Up to an overall sign. 
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one restricts the topological spinfoam sum (8.26) - with G = Spin{4) - to the so- 



called Spin{A) simple representations: recall that Spin{A) = SU{2)xSU{2), thus 
you can label Spin{4) irreducible representations with couples of spins (j"^, 
simple representations are the ones with j+ = j' . 
The spinfoam partition function results in 



E ndi-(p/)n 



simple p f 




].39) 



where, as just said, the sum is over simple Spin{4) representations; BC is a fixed 
intertwiner operator between (and depending only on) the 4 representations con- 
verging in it (following the 4-simplex pattern). The vertex amplitude ultimately 
depends on 10 spins, and it is thus referred to as the lOj-symbol. Of course, 
thanks to the 'simplicity' of the representations, we can rewrite this partition 
function as a sum over SU (2) irreducible, using the following property of the BC 
intertwiner 




tl,...,t5 




H J34 t3 




namely 



if V ti,...,t5 




H J34 1-3 




Notice that, according to (8.26), the face amplitude has been set equal to the 



dimension of the face representation. Here this representation is simple, so 
dim(j®j) = (2j + l)2. 
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However, this happens to be a subtle point, since we know (e.g. from QFT) 
that constraints change the integration measure of the path integral. Thus the 
use of the topological face amplitude is not justified at all in dealing with Quan- 
tum Gravity spinfoam models. This is a key point, regarding the choice of the 
face amplitude, and it will be treated more carefully in section |9| 

We have sketched briefly the structure of the Barrett-Crane model. However 
it must be said that this spinfoam model, despite all its success in the few 
years after its first appearance [189J, has been proved to be unfit to Quantum 
Gravity |190| : the big problem is that the boundary states of the BC model are 
only a small subgroup of the spin network states. This is due to the way by 
which constraints are imposed. 

This is the main reason for which much work has been done in order to ame- 
liorate the model. There have been a few proposal starting from 2007. In the 
following section we present the model which is considered the present day best 
candidate to attempt doing calculations in Quantum Gravity. 



8.3.2. The EPRL model 

Here I present the Engle, Pereira, Rovelli, Livine proposal for the vertex ampli- 
tude of spinfoam models. I refer the reader to the original papers |192l 11931 11941 
1195] for a thorough discussion. 

I don't want to give here a rigorous (and cumbersome) derivation and jus- 
tification of the model, I just want to stress the key points. We start with a 
triangulation A with the usual association of group elements to (dual) edges 
etc. . . (see section 8.2), but we use now the General Relativity action in its Hoist 



form (6.47), which can be compactly written as 



S[e,Lo] = j (^(eAe)* + ^(eAe)^ AF(a;) . (8.42) 
with * the Hodge dual operator, namely 

Fh = leuKLF^"- . (8.43) 

This action has the merit of including the Immirzi parameter 7 of Loop Quantum 
Gravity as a pre-factor of the topological sector, which has no consequences on 
the equations of motion (they remain the Einstein field equations) but allows 
the formulation in terms of Ashtekar-Barbero variables. The program is, as for 
the Barrett-Crane model, to write then action as ^BF + constraints' then write 
the BF spinfoam sum and only at the end apply the constraints. The difference 
with respect to BC are in the starting action and (more important) in the way of 
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imposing the constraints. Here we shall work directly in the Lorentz framework, 
thus with gauge group SL{2,C) rather then 5*0(4). 



Hoist action (6.47), (8.42) can be written as 



S[e,uj] = j (^B+ ^B*^ A F{uj) + constraints. (8.44) 

where the constraints must impose, as in BC, B = (eAe)*. The discrete variables 
shall be, as previously, Bf & g associated to triangles/dual faces and associated 
to tetrahedra/dual edges. We call Uf the holonomy around the face /, i.e. the 
product of the group elements of the edges bounding /. 

This time I shall give a more detailed derivation of the EPRL formula. 



We begin by giving a discretized form of the simpHcity constraints (8.37). We can 



easily discretize the constraints in the form (8.37) as 



B*j ■ Bf, 
B*f ■ Bf, 



, 
, 

±121^ 



;.45a) 
.45b) 
;.45c) 



where: V is the volume of the simplex; in the second equations / and /' are faces 
sharing an edge (equivalently: they are associated to triangles living on the same 
tetrahedron) while in the third they are don't (equivalently: they are attached to 
triangles belonging to two distinct tetrahedra). Actually, the correct way to deduce 



the discrete Bf variable from the 2-form B^i, is the followin, 



.46) 



Instead of requiring the last constraint (8.45c), we take the following closure constraint 

Bh+Bf,+Bf.^+Bf,=0, (8.47) 



for the four faces sharing an edge. It is easy to check that diagonal+off diago- 
nal+closure is an equivalent system of constraints. 

Remark. In the BC model the diagonal constraint implies the face representation to 
be simple, while the off diagonal one implies the uniqueness (and the specific form) of 
the BC intertwiner. Both constraints are imposed strongly. 

^Actually the situation is a little more tricky, I refer to the good paper by Engle, Pereira and 
Rovelli [193 for a thorough explanation of this point. 
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We have seen that these constraints admit more solutions than GR, precisely a trivial 
topological sector. Now we exclude the trivial sector by imposing a slightly different 
form of the constraints, i.e. we require that for each tetrahedron exist a vector nj such 
that, for each triangle of the that tetrahedron, holda^' 



ni{B}Y'^ = . (8.48) 



This constraint is intended to replace the off diagonal simplicity constraint ( 8.45b[ ). 
Geometrically the nj represents a vector normal to the tetrahedron/edge. 

Having cleared the meaning of the discrete versions of the simplicity constraints 



( 8.35 1 , we rephrase them in a convenient way. The conjugate momenta to the holonomies 



are 



Which, once inverted 



jf = Bf + ^B) . (8.49) 



= ( J [Jf --Jf] ^ (8.50) 



thus we can reformulate the constraints as 



Cff = J}Jf[^-^)+^JfJf = 0, (8.51) 



Cj = nj[ J* + ) = . (8.52) 



Now we choose a specific nj. A typical choice is nj = 6% that means that all the 
tetrahedra are spacelike. In other words, we are selecting a specific SU{2) subgroup 
of the full SL{2, C). Obviously this is a gauge choice, and must not influence physical 



results. With this choice the constraint (8.52) becomes 



C} = ieVf + ^^?=^} + ^i^} = 0, (8.53) 

where L are the generators of the SU{2) subgroup of SL{2, C) that leaves nj invariant; 
K are the generators of boosts in the nj direction. 

Let us now deal with the quantization of the constraints. In order to do that we must 
identify a Hilbert space in which we define operators. Taking a single vertex bubble, 
recall that the graph on its boundary 7^, naturally defines the following boundary 
Hilbert space 

L\SL{2,C)^) , (8.54) 

where L, N are the number of links and nodes of 7^, (cut from the spinfoam by the 
3-sphere). We impose the constraints as follows. 



^Indeed the off diagonal simplicity constraints imply that the triangles of each tetrahedron 
lie on a common hypersurface. If they are satisfied, there will be a direction orthogonal 
to all the faces. 
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The closure constraint (8.47) imposes 5^(2, C)-invariance on this space, i.e. it im- 
plements the (usual) quotient 

l2(5L(2,C)^/SL(2,C)^) . (8.55) 

The simplicity constraints (8.51), (8.53) are defined on faces, thus they act each on a 
single copy of the group, namely on L'^{SL{2, C)). The diagonal simplicity constraint 
(8.51) on this space reads 



Ci 1 



1 



+ -Ci" 

7 



"0 



with Ci and C2 the Casimir operator of g, i.e. 

Ci = J ■ J = 2{L^ - K'^); C2 = J*-J 



-4L-K 



].56) 



.57) 



The quotations mark means that we have to decide how to impose this constraint. 
However notice that having expressed it in terms of Casimir operators eigenvalues, it 
commutes with all the other operators, thus we can impose it strongly, i.e. requiring 
it to annihilate physical states. 

The constraint (8.53), on the other hand, is more subtle. The technique used in 
|193) . first proposed by Thiemann |196 ] is to pack them into a master constraint 



M 



f 



^(C7i)2 u ^ a 



.58) 



Classically it is of course equivalent to imposing the equal zero separately. The 
power of this approach, is that Mf is now a combination of Casimirs 

1 \ Ci C2 



1 + 



r 



272 27 



. 



Combining (8.56) and (8.59) we get the following (definitive) set of 2 constraints: 

C2 



1 



1 



7 



+ -Ci " = " , 

7 



C2 - 47L2 " = " 



^.59) 

,60) 
,61) 



Having done this, we now simply have to see what consequences these constraints have 
on states of L'^{SL{2,C)). It is easy to see that, labeling with {p,k) {p real and k 
half-integer) the S'L(2, C)-irreducible representations, thus having the decomposition 



L^{SL{2,C)) = Y,'H(p,k) 

{P,k) 



(p,k) 1 



n 



0/ 

j'=k 



- with j' in the sum denoting the SU{2) spin-/ irreducible - of SL{2, 
into S'C/(2)-irreducibles, the constraints impose, for each face, 

Pf = ijf ' ^/ = jf > 



(8.62) 
)-irreducibles 
(8.63) 
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for some half-integer jf, and they restrict the decomposition to the lowest spin, j' = 
k = j. 

You see that the constraints have selected a SU{2) subgroup of SL{2, C)! Precisely 
they tell two things: 1. the permitted 5L(2, C) face labels are only the irreducibles 



of the type (pf = 'yjf,kf = jf) for some half-integer jj^^. 2. given the face label 
(jjfjjf) they select out the spin-j SU{2) irreducible. Let's give a precise definition of 
the embedding L2(SC/(2)) ^ L'^{SL{2,C)): 

Y:j^jC n^^j^j) C L\SL{2, C)) , (8.64) 

is the map that sends each SU (2) spin j irreducible in the lowest spin (i.e. - of course 
- the spin j) irreducible inside (7i, i). Thus, the Y map takes each state in Lp'{SU{2)) 
to a state of L^(5L(2, C)). In terms of Wigner matrices it is really simple 

Di ^ YD^Y^ = D(^J'J) . (8.65) 

This is very nice indeed, since we can define boundary variables to be SU{2) states, 
and this is just what we expect to have, since on the boundary we want to put spin 
networks. 

Remark. I have done all the calculation for SL(2, C) since it is the "reality". How- 
ever, as I had occasion to say, it is sometimes useful to see what happens in the simpler 
euclidean case, i.e. for 50(4). I do not repeat the discussion - which proceeds very 
similarly - the result is the following: the simplicity constraints reduce the 50(4) irre- 
ducibles {jpj'j) to the ones given by ("f+jfj^-jf), jf 'running' over SU{2) irreducibles 
and with [193] 

|l±7l 

7± = • (8.66) 

Thus, in the euclidean 50(4) case we have an embedding of SU{2) in 50(4), quite 
similarly to what happens in the lorentzian case. 

We are almost done. We just have to see the consequences of what I have just said 
in terms of concrete formulae. 



Let us focus on a single vertex. The idea is to take the BF vertex amplitude (8.29) 



(or equivalently the 4-simplex spin network evaluation, i.e the graph on the far right 



of (8.26)) and to restrict the boundary variables to satisfy the simplicity constraints. 



in the form we have just said: 



A^^^^{Ui)= [ ll^9nllP(Ui,gsii)9;^)), (8.67) 

JSL{2,C)« „ J ^ ' 



10 



Notice that this also implies that the effective sum over p, which should be an integral, 
restricts to a sum. 
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where we have no more a delta inside the face product, since now Ui is an SU{2) 
element, but a 'generalized' delta, i.e. 

PiU,g) = Y,{2j + l)TT{D^{U)Y^D^^^'^\g)Y) , (8.68) 
i 

with U G SU{2), g £ SL(2, C). Performing the SU{2) integrals, just as in the standard 



BF theory case (8.2), one can easily find the spin/intertwiner representation of the 



vertex amplitude done in the BF theory case. 



One ends up with 

(8.69) 

where y^^" ^ ^ = {i(p.k) \ Y Putting together all the vertex amplitudes we get 

where the last parenthesis indicates that some simplices are cut by the boundary, 
thus some spins and intertwiners will be fixed (and not summed) to the boundary 
value, given in the left hand side 



11 



Notice again that we have no clue about the face amplitude: should it be 
the SU{2) dimension of the representation attached to that face (as in a SU (2)- 
BF theory)? Should it be the SL{2, C) /SO{A) dimension (as in a SL{2, C) /SO{A)- 
BF theory)? The imposing of constraint changes the measure of the path inte- 
gral, and we have no control about that in this derivation. 



8.4. Spinfoam: a unified view 

After all this rather cumbersome (and quite chronological) model-making I would 
like to give a more coherent and compact view of spinfoam approach and of its 
relation with Loop Quantum Gravity. As pointed out in |154] it is time to provide 
a "top-to-bottom" framework, in which set some properties we want these models 
to satisfy and deduce from them specific models. This has (at least) the merit of 

^^The partition function is obviously a particular case of this formula, i.e. 

Zi^^^ = Y^Y[Aj{jf)\{A'^^"\j^,) . (8.71) 
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clarifying the mess of technicalities that plague - in my opinion - the spinfoam 
model-making approach. 

Remark. What follows is just an attempt to give a unified and more coherent 
framework of spinfoam approach, but it is not a well established and rigorous 
"chapter" of Loop Quantum Gravity. You should take what follows as a reasoning 
ex post on how we can insert all the models into a single framework. For this 
section I mainly follow the paper |154] . 

First of all: what do we want from a spinfoam model? We want a way to 
compute transition amplitudes 

{s\A\s') , (8.72) 
between Loop Quantum Gravity spin network states. A represents the projection 



operator into the kernel of the hamiltonian operator (recall equation (8.4)). In a 
covariant 4d context this is rephrased: take a 4d manifold with boundary state 
ip, we want to calculate 

(AI^P). (8.73) 

The boundary state i/j is typically formed by two spin networks (in terms of 
graphs = U F^j) the "initial" and "final" spin networks, but we can 
be as generic as we want. This last expression is the amplitude associated to 
the specific state ip. The linear functional A (or, if you want, (0|y4, but here 
we enter in the subtle (and fascinating) issue of vacuum in Quantum Gravity, 
maybe we will spend some words about it in the rest of the chapter) is the heart 
of the spinfoam model, and must be thought of as an evolution operator. What 



properties do we ask (8.73) to have? 

• Superposition. Namely we want 

{A\^) = J2M^), (8.74) 

i.e. we want that the amplitude is expressed as a "sum over histories". This 
is what we argued with a rather heuristic touch in section |8} Obviously 
the specific set over which to sum is at this level completely undetermined. 

• Locality. This is a request on the single history-amplitude A{a). We require 

A{a)^l[A,. (8.75) 

This amounts to ask that each history has an amplitude that is the product 
of elementary amplitudes, or vertex amplitudes (the "atoms" we have been 
talking about). 
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• Local Lorentz invariance. Recall that classical general relativity, in tetrad 
formulation, has a local Lorentz invariance, namely a SL{2, C) gauge in- 
variance. However the boundary states know nothing of SL{2, C), they are 
built up over SU{2) gauge invariance. Thus there must be some embed- 
ding from S'?7(2)-gauge invariant states to SL{2, C)-invariant ones. This 
map, usually called simply / is what actually determines the spinfoam 
amplitude. 

Now we have the key ingredients. We have to focus on a single vertex, sur- 
rounded by a kinematical state ipj. We define 

A,{I) = {A,\ijj) = , (8.76) 

with 

/ = PsL(2,c) o Y (8.77) 
the map that embeds the boundary state in L^(5'L(2, C)^) and then projects 
onto the ^L(2, C)-invariant on the nodes, i.e. onto L^{SL{2,C)^ / SL{2,C)^)- 
stands for a state in the boundary Hilbert space of the vertex v, and / 
denotes a set of quantum numbers (think of a spin network state). 
Writing it a bit more explicitly 

(A IV/) = / Y\dg^{l\U{g„)Y\iji) . (8.78) 

ipj is a generic boundary state: it could be a spin network state \ji,in) giving 
Av{ji,Ln) or, in the holonomy('position') representation, \Ui) giving Ay{Ui). (1| 
stands for the SL{2, C) holonomy state where all the links are set to the identity 
(it corresponds to the evaluation of the state on which it acts). U{g) is the action 
of g E SL{2, C) over the state Y\ipj). In a suggestive way one could write 

{A,\= [ T\dgMU{g^)Y , (8.79) 

JSL{2,C)'^ „ 

as a functional over the kinematical Hilbert space of Loop Quantum Gravity. Of 



course this definition gives the EPRL formula (8.69), when {ipj) = \ji,Ln)i and 



gives exactly (8.67) in terms of holonomies. 

Moreover, this framework is completely general, in the sense that it is able to 
reproduce every spinfoam model that matches Loop Quantum Gravity kinemati- 
cal Hilbert space on the boundary. Obviously various models can be reproduced 
by a choice of the / function, i.e. of the Y map, which is, ultimately, what 
determines the specific vertex amplitude, and - indeed - is itself completely de- 
termined by how we impose the constraints on the BF model (recall that there 
are infinite SU{2) irreducibles inside each SL{2,C) irreducible). 
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Once got the vertex amplitude one get the total amplitude by 

(T V 

notice that this sum is two-fold: fix a triangulation and you have to sum over the 
labeling of internal faces and edges; then you have to sum over all the possible 
triangulation in order to really catch all the possible "paths" bounded by ipjl 
This make apparent the problem in dealing with triangulation dependence and 
in how to weight different triangulations (the attentive reader was surely already 
aware of this problem since the very beginning of the discretization stuff). I will 
briefly review the topic (which is a totally open issue) in the following section. 

Moreover here there is no hint about something attached to the faces of each 
spinfoam. We know from the BF theory derivation that something there must 
be, and that it very likely has to do with the dimension of the representations 
attached to the faces (see chapter |9| . 

8.5. Spinfoams as a field theory: Group Field 
Theory 

The present section is a bit beyond the scope of this review of Loop Quantum 
Gravity and spinfoam models, particularly since to understand the research I 
did in this field ^ the argument of this section is not strictly necessary, nor - I 
must admit - closely related. 

However, I spent quite a long time in Marseille studying Group Field Theory, 
and I seriously think it has more than a few chances of becoming a bridge 
between Quantum Gravity in its spinfoam formulation and the prolific and well- 
established framework of Quantum Field Theory. 

In a sense that will be clear - I do hope - from what follows. Group Field 
Theory is precisely a way of putting the spinfoam models in the framework of 
QFT, a very special kind of QFT. 

GFT was first introduced by Boulatov in 1992 |197j . and had, since then, a 
very promising development as a Quantum Gravity framework. The nowadays 
work in this topic is mainly to study the derivation of the spinfoam models 
(particularly the EPRL model) in this framework; to study quantum corrections 
and renormalization issues; control the sum over triangulations; all done with 
QFT-like tools. For a panoramic of the present-day research see the following 
papers |199l [198], while for good reviews of GFT see |200l [2011 1202] . 

The very basic idea is the following: find a field theory whose Feynman dia- 
grams are spinfoams. The natural ancestor of GFT are the well known matrix 
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models: their Feynman diagrams are ribbon-like, which can be seen as dual to 
a triangulated surface [203^ I2U4| . 

A matrix model is a model whose action is something like 

S[M] = yTr(M2) + ^Ti{M') , (8.81) 

where M is a x matrix. I have taken a potential term of order 3, but of 
course it is just an example. 

Going directly to the quantum side of the story, let us try to calculate the 
partition function for this action 

Z = JvM e-^[^^l = jvM e-f Tt^' ^ 1 (^^^ " (TrM^)" . (8.82) 

What are the differences from standard QFT? There are (at least) two: 1. it is 
non local, since our 'fields' Mij are here objects which depend on two 'points'; 
2. it is not a field theory, since the 'space' variables are discrete, i, j = 1, . . . , N . 

We shall drop the non-field stuff in a moment; the true 'news' is the dropping 
of locality, which might sound rather blasphemous, but recall that we are in a 
completely combinatorial framework and these fields do not pretend to represent 
causal propagating particles, but blocks of spacetime, as we shall see . Deriving 
the propagator and interaction vertex of the theory is straightforward 

Gijki = Sij6ki , (8.83) 

k I 



Vijklmn — ^jm^kl^mn . / \ ' (8.84) 




Now imagine to glue interaction vertices through propagators to obtain Feynman 
diagrams, you will get - as hinted above - ribbon-like structures. If we take the 
dual of each of these graphs you obtain precisely a triangulation of a 2-manifold. 
In this sense each propagator represents a side of a triangle, while the vertex 
represent the triangle itself, namely 
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We can rephrase this by saying that each propagator is an edge and each vertex 
a face of the dual 2-skeleton of a triangulation. This gives also a powerful tool 



for calculation explicitly the partition function (8.82), indeed it has been proven 
that it is given by 

Z = J2^9WN^~^' ^ (8-85) 

9 

where g is the genus of the triangulated surface and w is a weight depending on 
the diagram symmetry factor. This is precisely the partition function of a 2d 



BF theory (8.24)! 



You guess the clue: going up in the number of indices - i.e. building tensor 
models \205\ 1206) - we can have the hope of getting skeletons of 3d and 4d 
manifolds, i.e. spinfoams. 

Actually, this is not the case. Tensor models has been proven to fail in this 
respect. Essentially: 1. they are too simple to incorporate even the 3d i?F theory 
case and 2. they produce both manifolds and pseudo-manifolds^^ |205l I206"| . 



The natural way to proceed, is to pass from matrices to fields, defined on an 
appropriate domain space. In GFT - hence the name - this domain space is 
chosen to be an appropriate number of copies of a compact group G 



= \ I d^id^2 </>^(^i,^2) + ^ / d^id^2d^3 0(^1,^2)0(^2,^3)0(^3,^1) , 

(8.86) 

this is the GFT analogous of the matrix model (8.81). Instead of analyzing its 



features, let us move on seeing the key features of a general GFT. 

• is a (typically) real valued field on G", with G compact Lie group and n 
will be the dimension of the 'produced' triangulated manifold. We require 
a gauge invariance under the G-(right)action 



1, 



I 9n/ 



^i9i9,- ■ ■ ,9ng) , V^„^ G G . 



i7) 



• the GFT action (for a real valued field) has the (very) generic structure 

5[0] = T[0] + Ay[0] , (8.88) 

with a "kinetic" part which has no dynamical meaning, since at this stage 
there is nothing as a time variable; and an interaction term which de- 
fines how the various group elements are attached with one another: the 
interaction term is actually the 'atom' of hidden spinfoam model, and it 
determines the kind of n-complexes by which the manifold is triangulated. 

Sort of triangulation of manifold with singularities |207| . 
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The Feynman diagrams are 2-complexes, whose combinatorial structure 
is entirely given by the interaction term of S[(f)]. They are interpreted 
as dual to a n-discretization in terms of n-complexes (if the interaction 
term is particularly simple - actually the only case I will consider - these 
complexes are indeed simplices, and the discretization is a triangulation). 

The value of the Feynman diagram is seen as a spinfoam amplitude, the 
amplitude associated with the discretization of that specific diagram. 

The GFT partition function 

Z= jv<P e-^[^l = jv<P e"^[*]e-^^[<^] = «;(A)A^(A)Aa , (8.89) 

which is the typical perturbative expansion in terms of Feynman diagram: 
A denote equivalently the Feynman diagram or its associated discretiza- 
tion; V^(A) is the number of vertices of the diagram, w{A) is its symmetry 
factor and is the value of the diagram (the spinfoam amplitude). 

Notice that this gives also a very natural way of dealing with triangula- 
tion dependence of spinfoam: each triangulation has a different spinfoam 
amplitude and we simply sum them with the weight given by the coupling 
and the symmetry factor. 

The correlation functions of GFT give the amplitude in going from one 
state to another, i.e. the (most wanted) (s | A | s') 

A[r,^]= J2 «^(A)A^(A)Aa = P^(0)e-^[<^] , (8.90) 

A|aA=r 

where ip is a state on the (boundary) graph F and P is some polynomial 
function of that is able to render the specific graph F labeled by variables 
so to reproduce ■?/ 
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Let us now review some GFT models, stating only the results 

• Ponzano-Regge. The 3d SU (2) GFT model with trivial kinetic term and 
tetrahedron-like interaction term, generates exactly the 3d BF path inte- 
gral ( [8^ for SU{2). The action is 

S[(p] = \ j (^91(^92^93 (t)^ {91,92, 9?) 

. 6 

n ^9i(t>{9i, 92, 93)^(93, 9i, 96)^i9e, 92, 94)^(95, 9a, 9i) ■ (8.91) 



4! 



=1 



^•^Think of n-point functions of QFT and you will get the idea. However, I will not go in much 
more details on this. I refer the reader to 11831 and references therein. 
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Figure 8.6.: Propagator and interaction term of the 3d GFT model (8.91). The 



propagator is simply a delta over the group, while the interaction 
term is defined by a tetrahedron-like contraction. 



The interaction term is easily understood looking to figure 8.6 It is in- 
deed a tetrahedron like graph by itself; moreover, if you take the dual in 
this sense strand — )■ edge (of a triangle), field (i.e. 3-strand) — )■ triangle, 
interaction — )■ tetrahedron. In this sense the interaction term is precisely 
telling us how to glue the 3 edges of 4 triangles to produce a tetrahedron. 

In this precise sense, gluing interaction terms through propagators you 
create Feynman diagrams as well as a 3d triangulation. 

Taking the Fourier representation (i.e. the Peter- Weyl decomposition) of 
the field you get for the vertex the Ponzano-Regge amplitude. 



4d GFT model. Take the following action 



S[4>] 



A 



dgidg2dg3dg4, 4>'^{gi, g2, gs, g^) + -^^ / Yldgi(f){gi, g2, gs, gA)x 



10 



i=l 



0(5-4, 55, 56, 57)0(57, 53, 58, 59)0(59, 56, 52, 5io)0(5io, 58, 55, 5i) • (8.92) 



Here the interaction vertex is a 4-simplex (see figure 8.7), thus Feynman 
diagrams are dual to a 4d triangulation, and the vertex value, once decom- 
posed into irreducible representations of the underlying group, is exactly 



the 4d BF vertex amplitude, as in equation (8.26). 



Barrett- Grane. As a final example I give a sketch of how it is possible to 



deduce the Barrett-Crane spinfoam model |8.3.1| from a GFT model. The 
gauge group is of course G = Spin{4). The trick is to introduce a (slightly) 
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10 629 9387 



Figure 8.7.: Interaction term of the 4d GFT model (8.92) 



more complicated kinetic term in (8.92), namely 

1 



m = - I d'g{P<P{g,,...,g,)y 



with 



P(f>{gi, ■ ■ ■ , gi) = ^g I 

JSL(2.C) Jsi 



i.93) 



dg I d u (j){gighi, g2gh2, g3gh3, g4gh4) . 

'SL{2,C) ^5(7(2)4 

(8.94) 

It is matter of calculations to prove that this GFT gives exactly the BC 



spinfoam amplitude (8.39) 



For a derivation of the EPRL spinfoam model from a GFT, I recommend the 
reading of |198] . 



Bubble divergences A very delicate and important issue of the spinfoam ap- 
proach to Quantum Gravity is the one regarding divergences. I want here only to 
inform the reader of the importance of the issue, and give some useful references. 
Firstly a remark: inside the LQG framework one always talks about infrared di- 
vergences, i.e. divergences arising when summing over the large-distance-scale 
degrees of freedom, i.e. large spins. In LQG there are no ultraviolet divergences, 
since there is an ultimate minimum length-scale, of the order of the Planck length 
(see section 7.5). 



138 



8.5. Spinfoams field theory: Group Field Theory 

Bubble divergences arise in the sum over spins. You can see it as a loop 
integral of QFT. Simple power counting techniques reveal the presence of diver- 
gences depending on the nature of the manifold and of the kind of vertex one is 
considering (how many edges converge to that vertex, how many faces, and so 
on). Many studies have been done in this respect (See |208l I209"| 1210] and refer- 
ences therein). Moreover, it has been noted that the face amplitude is crucial in 
determining whether a bubble is or is not divergent [210j . We shall talk about 
this in chapter [9j 

Here I want to stress that the GFT approach gives powerful tools to handle 
and analyze these divergent bubbles and may provide a unified strategy to fix 
themp^ I refer the reader to the recent work |211) and to |212l I213| . 



This concludes this brief detour into GFT framework. As I have said at the 
beginning of this section, GFT is not the focus of my research papers, so this is 
not the right place to discuss it in more details. What I wanted to give, is the 
idea of the existence of a (in my opinion) promising unifying framework, inserting 
Quantum Gravity into a QFT approach while preserving its non-perturbative 
nature. Indeed recall that the perturbation of the Feynman diagram series in 
GFT is a perturbation in the coupling A, i.e. is a perturbation on the number 



of vertices of the triangulations, that is on the "complexity" of the grapf and 



not a perturbation around a fixed background geometry. I really encourage 
the interested reader to look up the bibliography I have referred to during this 
section, in particular \200\ \201\ I202| for good reviews on the subject. 



Incidentally, notice that the GFT approach provide with an intriguing duality between UV 
and IR divergences: indeed the Quantum Gravity IR divergences (namely the large spin 
divergences) can also be seen as UV divergences on the group. 

Of course we have no clue on the coupling parameter, but this is another problem. 
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amplitude in Quantum 
Gravity spinfoam models 

The present chapter is intended to present the content of the paper JB], in which 
we propose a way for fixing the face amplitude of a general spinfoam model. The 
proposal is motivated by the requiring of a sort of "unitarity" of time evolution of 
space geometries, in a sense that will be made clear in what follows. We found 
that this requirement imposes the face amplitude to be equal to the dimension of 
the SU (2) -projected representation of the S0{A)(SL{2,C)) one attached to the 
face. 

9.1. Introduction and resume of the content of 
the paper ||5]| 



In this section I present the content of the paper [5], which collects the results 
of the work I have done in Marseille, in collaboration with Eugenio Bianchi and 
Carlo Rovelli. 

I have repeatedly focused on the fact that, while for BF theory the face ampli- 
tude of the spinfoam sum is well determined by the path integral discretization 
procedure and it is given by the dimension of the representation labeling the 



face (8. 24), (8. 25), (8. 26), for Quantum Gravity the situation is much more sub- 
tle. One starts with a certain BF theory, obtains a spinfoam sum formula, and 
then imposes constraints in the vertex amplitude. What happens to the face 
amplitude? The BC model, for instance, simply let it be the BF face ampli- 
tude, i.e. the square of the SU{2) irreducible dimension. In S'0(4)-models with 
the Immirzi parameter (such as S'0(4)-EPRL) one ha^ 

Ap^ = (2j+ + l)(2j_ + 1) = (27+ + l)(27_j; + 1). (9.1) 



^In [S] we have worked in the eucHdean (5*0(4)) case. However, everything can be done in 



the lorentzian (5i(2, C)) as well, mutatis mutandis. See section 8.3.2 for discussions about 
this point. 
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However this is not well motivated, indeed doubts can be raised against this 
argument. For instance, Alexandrov |i221j has stressed the fact that the imple- 
mentation of second class constraints into a Feynman path integral in general 
requires a modification of the measure, and here the face amplitude plays pre- 
cisely the role of such measure, since we have an independent 
way of fixing the face amplitude? 

Let me recall that all the spinfoam models share the same form of partition 
function, namely the sum 

^^ = E n^P/ U'^M^^e) , (9.2) 



where, as usual (recall section 8.2), that sum is intended to be a sum over all 
the possible labellings of the faces (edges) of A* with irreducible representations 
(intertwiners) of an appropriate group G. 

In pj we argued that the face amplitude is uniquely determined for any spin- 



foam sum of the form (9.2) by three inputs: 1. the choice of the boundary 
Hilbert space, 2. the requirement that the composition law holds when gluing 
two-complexes; and 3. a particular "locality" requirement, or, more precisely, a 
requirement on the local composition of group elements. 

We argued that these requirements are implemented if the partition function 
Z is given by the expression 

Z^= f dUJ l[A,{UJ) n mT-UT) , (9.3) 

V f 

where UJ G G, Vi...Vk are the vertices surrounding the face /, and A^{UJ) is the 
vertex amplitude A^{jf, i^) expressed in the group element basis [222J. Then we 



showed that this expression leads directly to (9.2), with arbitrary vertex ampli- 
tude, but a fixed choice of face amplitude, which turns out to be the dimension 
of the representation j of the group G, 

Aj = dim(j) . (9.4) 



In particular, for Quantum Gravity this implies that the BF face amplitude (9.1 ) 
is ruled out, and should be replaced (both in the Euclidean and in the Lorentzian 
case) by the SU (2) dimension 

Aj = dim(j) = 2j + 1 . (9.5) 



Equation (9.3) is the key expression of the whole paper. 



I organize this chapter as the paper from which is taken, i.e.: in section 9.2 I 



show that 5*0(4) BF theory (the prototypical spinfoam model) can be expressed 
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Figure 9.1.: Schematic definition of tlie group elements h^e, UJ and Ue associated 
to a portion of a face / of the two-complex. 



in the form (9.3). Then I discuss the three requirements above and I show that 
(9.3) implements these requirements. (Section 9.3[ ). Finally I show that (9.3) 
gives (9.2) with the face amplitude (9.4) (Section IV). 

The problem of fixing the face amplitude has been discussed also by Bojowald 
and Perez in [220j. Bojowald and Perez demand that the amplitude be invariant 
under suitable refinements of the two-complex. This request is strictly related 
to the composition law that we considered in [5], and the results we obtained 
are consistent with those of |22U] . 



9.2. BF theory 

Take the general expression of the BF partition function in terms of group 



elements (see equation (8.16)) 



(9.6) 



where Ue are group elements associated to the oriented edges of a, and (ei, Cat) 
are the edges that surround the face /. Let us introduce group elements h^e, 
labeled by a vertex v and an adjacent edge e, such that 



(9.7) 



where v and v' are the source and the target of the edge e (see figure 9.1 ). Then 



we can trivially rewrite (9.6) as 



(9. 
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Now define the group elements 



(9.9) 



associated to a single vertex v and two edges e and e' that emerge from v and 



bound the face / (see figure [93] ) . Using these, we can rewrite (9.6) as 



where the first product is over faces f" that belong to the vertex v, and then a 
product over all the vertices of the 2-complex. 

Notice that this expression has precisely the form (9.3), where the vertex 
amplitude is 

''dK,l[6{uj,K,h-^]) , (9.10) 
/" 

which is the well-known expression of the 15 j Wigner symbol (the vertex am- 
plitude of BF in the spin network basis) in the basis of the group elements 
(cfr. (18^291). 

We have shown that the BF theory spinfoam amplitude can be put in the 
form (9.3). We shall now argue that (9.3) is the general form of a local spinfoam 
model that obeys the composition law. 



9.3. Three inputs 

(a) Hubert space structure. Equation ( |9.2[ ) is a coded expression to define the 
amplitudes 

AAiJhLn) =^YlAj^YlA^{jf,Le;jl,Ln) , (9.11) 

defined for a triangulation A with boundary, where the boundary graph F is 
formed by links / and nodes n. The spins ji are associated to the links /, as well 
as to the faces / that are bounded by /; the intertwiners i„ are associated to the 
nodes n, as well as to the edges e that are bounded by n. The amplitude of the 
vertices that are adjacent to these boundary faces and edges depend also on the 
external (thus fixed) variables (ji,6„). 

In a quantum theory, the amplitude A{ji, t„) must be interpreted as a (covari- 
ant) vector in a space Hr of quantum states]^ We assume that this space has a 

^If r has two disconnected components interpreted as "in" and an "out" spaces, then Hj^ 
can be identified as the tensor product of the "in" and an "out" spaces of non-relativistic 
quantum mechanics. In the general case, Hr is the boundary quantum state in the sense 
of the boundary formulation of quantum theory |151L 1223] . 
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Hilbert space structure, which we know. In particular, we assume that 

nr = L2[G\dUi] , (9.12) 
where L is the number of hnks in F and dUj is the Haar measure. Thus we can 



interpret (9.11) as 



AAiJhtn) = {jl,Ln\A) 



where \ji, in) is the spin network function (cfr. equation (7.28)) 



(9.13) 
(9.14) 



Using the scalar product defined by (9.12), we have 



dUi ^,,„XUi)^f,yAUi) 
JJdim(jO 5j^f^ YlS,„,>^ . 



(9.15) 



where dim(j) is the dimension of the representation j. Therefore the spin- 
network functions ipji,i,„{Ui) are not normalized. (These dim(j) normalization 
factors are due to the convention chosen: they have nothing to do with the 



dimension of the representation that appears in (9.4).) The resolution of the 
identity in this basis is 



1 = E (ndim(jO) \ji,tn){j 



I, in 



(9.16) 



(b) Composition law. In non relativistic quantum mechanics, if U{ti, to) is the 
evolution operator from time to to time ti, the composition law reads 

U{t2,to) = Uit2,t,)U{t^,to) . (9.17) 

That is, if \n) is an orthonormal basis, 

(/ \U{t2, to) \t)=J2{f \U{t2, t^) I n) (n \U{t,, to) \ i) . 

n 

Let us write an analogous condition of the spinfoam sum. Consider for simplicity 
a two-complex a = Ui U (T2 without boundary, obtained by gluing two two- 
complexes a I and cr2 along their common boundary F. Then we require that W 
satisfies the composition law 



(9.18) 
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- where now I label Z directly with the 2-skeleton cr C A* - as discussed by 
Atiyah in |224] . Notice that to formulate this condition we need the Hilbert 
space structure in the space of the boundary states. 

(c) Locality. As a vector in H^, the amplitude A{ji,Ln) can be expressed on 
the group element basis 

A{Ui) = {Ui\A) = J2 ili^Mji)) ^,,.„(f/0^(jMn) . (9.19) 

Similarly, the vertex amplitude can be expanded in the group element basis 

A,{UJ) = {UJ\A,) (9.20) 
= E (ndim(j7))^,^,,(f/J)^(j7,C)- 

Notice that here the group element UJ and the spin jj are associated to a vertex 
V and a face / adjacent to v. Similarly, the intertwiner is associated to a 
vertex v and a node n adjacent to v. Consider a boundary link / that bounds 



a face / (see figure 9.2). Let Vi...Vk be the vertices that are adjacent to this 
face. We say that the model is local if the relation between the boundary group 
element Ui and the vertices group elements UJ is given by 

Ui = Uf ... Uf . (9.21) 

In other words: if the boundary group element is simply the product of the group 
elements around the face. 



Notice that a spinfoam model defined by (9.3) is local and satisfies composition 
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law in the sense above. In fact, (9.3) generalizes immediately to 



AM) = I duj hamj) n m?-UT) 

V internal / 

X n mr-U^Ur'). (9.22) 

external / 

Here the first product over / is over the ("internal") faces that do not have an 
external boundary; while the second is over the ("external") faces / that are 
also bounded by the vertices Vi,...,Vk and by the the link /. It is immediate 
to see that locality is implemented, since the second delta enforces the locality 



condition (9.21 ). 

Furthermore, when gluing two amplitudes along a common boundary we have 
immediately that 

j dUi Ajui) A^.iUi) = Z^iua, , (9.23) 

because the two delta functions containing Ui collapse into a single delta function 
associated to the face /, which becomes internal. 



Thus, (9.3) is a general form of the amplitude where these conditions hold. 

In |220j . Bojowald and Perez have considered the possibility of fixing the face 
amplitude by requiring the amplitude of a given spin/intertwiner configuration 
to be equal to the amplitude of the same spin/intertwiner configuration on a 
finer two-simplex where additional faces carry the trivial representation. This 
requirement imply essentially that the amplitude does not change by splitting a 



face into two faces. It is easy to see that (9.3) satisfies this condition. Therefore 



(9.3) satisfies also the Bojowald-Perez condition. 



9.4. Face amplitude 



sufficient to insert (9.20) into (9.3). This gives 



Finally, let us show that (9.3) implies (9.2) and (9.4). To this purpose, it is 



= duj]lJ2iU^^^iff))H^duj)Aifpii) 



X 



n mr-uT) ■ 



(9.24) 
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Expand then the delta function in a sum over characters 



Z. = [ dUJ nE(ndini(j7)) i^,^.,^M)MjhO (9.25) 

X H 5^dim(j^) TT{D^f{U}^)---D»{U}^)) . 

f if 

We can now perform the group integrals. Each UJ appears precisely twice in 
the integral: once in the sum over jj and the other in the sum over jj. Each 
integration over the group, gives a delta function (recall equation 



.21)) 



1 



(9.26) 



which can be used to kill the sum over jj dropping the v subscript. Notice that 
the dimensional factor involved here, exactly cancels - once done the product- 



the first dimensional factor in (9.25), i.e. the one coming from the spin network 
normalization. Following the contraction path of the indices, it is easy to see 
that these contract the two intertwiners at the opposite side of each edge. Since 
intertwiners are orthonormal, this gives a delta function 6,v ,v' which reduces the 
sums over intertwiners to a single sum over '■= = i-n ■ Bringing everything 
together, we have 



(9.27) 



This is precisely equation (9.2), with the face amplitude given by (9.4). 



Notice that the face amplitude is well defined, in the sense that it cannot be 
absorbed into the vertex amplitude (as any edge amplitude can). The reason 
is that any factor in the vertex amplitude depending on the spin of the face 
contributes to the total amplitude at a power k, where k is the number of sides 
of the face. The face amplitude, instead, is a contribution to the total amplitude 
that does not depend on k. This is also the reason why the normalization 
chosen for the spinfoam basis does not affect the present discussion: it affects 
the expression for the vertex amplitude, not that for the face amplitude. 

By an analogous calculation one can show that the same result holds for the 



amplitudes W: equation (9.11) follows from (9.22) expanded on a spin network 
basis. 



In conclusion, we have shown that the general form (9.3) of the partition func 



tion, which implements locality and the composition law, implies that the face 
amplitude of the spinfoam model is given by the dimension of the representation 



of the group G which appears in the boundary scalar product (9.12). 
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In general relativity, in both the Euclidean and the Lorentzian cases, the 
boundary space is 

Hr = L2[SU{2)'',dUi] , (9.28) 

therefore the face amplitude is dj = dimsu{2){j) = 2j + 1, and not the 5*0(4) 
dimension (9.1), as previously supposed. 

Notice that such dj = 2j + l amplitude defines a theory that is far less divergent 
than the theory defined by (9.1). In fact, the potential divergence of a bubble 
is suppressed by a power of j with respect to (9.1). In [210], it has been shown 
that the dj = 2j + 1 face amplitude yields a finite main radiative correction to 
a five-valent vertex if all external legs set to zero. 
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Conclusion 



Let me now briefly review in a bird-flight manner the results which have been 
achieved. 

We have worked on phantom dark energy models in classical Cosmology. Let 
me recall that phantom models (i.e. models with some kind of phantom energy) 
are not ruled out by observations, even if the standard model of Cosmology, 
ACDM, scientifically still "wins", by virtue of Occam's razor. Crossing of the 
phantom divide line is not excluded as well. We have analyzed what seem to 
us one of the best candidates to reproduce this crossing: two-field cosmological 
models with one scalar and one phantom field. 

We have discovered that, starting from an expression for the Hubble variable 
h{t), there is an infinite number of models (i.e. of potentials) that have, among 
their solutions, precisely that h{t). This freedom, namely to have infinite dif- 
ferent models which are compatible with the same dynamics for the universe, is 
the most important result presented in [1]. 

One possible approach to discriminate between such models, is to couple the 
cosmological fields with other, observable, fields, such as cosmic magnetic fields. 
The calculations and the numerical simulations have indeed shown that such a 
coupling gives in principle observable results, capable of select between different 
models [2j. 

However, we are aware that phantom fields have been heavily criticized, and 
actually they do deserve the label of 'exotic'. The negative sign of the kinetic 
term is indeed something that disturbs even just to write it. Thus, we tried to 
find out some mechanism so that the negative sign was only effective, while the 
"fundamental" field is a standard one, with no stability problem. Accepting to 
work in the PT symmetric Quantum Theory framework, we succeeded in this 
respect, and find a possible way to have an effective phantom field, stable to 
quantum fluctuations |3]. 

These were the contributions in the phantom/crossing-of-the-phantom-line 
branch of Cosmology. Aside, we analyzed one-fleld models giving evolutions 
for the universe with an interesting kind of singularities: sort of Big Bang/Big 
Crunch with a flnite non-zero radius [4j. Notice that this is in line with all 
Quantum Gravity theories, namely to have a minimum length scale. Thus, we 
thought that the study of non-zero radius singularities could be useful also in 
perspective of future results from the Quantum Gravity community. 
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Conclusion 

Quantum Gravity is essential to properly understand what is space and what 
is time. Maybe the total absence of experimental results give to these kind of 
studies an halo of philosophical vagueness. I could agree with this, but still I 
firmly think that this kind of fundamental research must be carried on. Indeed, 
many results of LQG and spinfoam theory are very suggestive and stimulating, 
even if on purely theoretical basis. The attempt to find a 'sum-over-histories' 
formulation of Quantum Gravity, based on the LQG approach to quantization, is 
a recent and promising field of research, where results crowd from many different 
and far areas of Theoretical Physics. Specifically we have proposed a way of fixing 
the face amplitude of a general spinfoam model. Up till now, the face amplitude 
has been taken to be the one of the BF theory underlying the quantization 
procedure. But we argued that this is not the case, if one want to assure the 
proper "gluing" of spinfoam amplitudes [5]. 
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